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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  MEMORANDUM  1337 


ANALYTICAL  STUDY  OF  SHIMMY  OF  AIRPLANE  WHEELS* 
By  Christian  Bourcier  de  Carbon 


SUMMARY 


The  problem  of  shimmy  of  airplane  wheels  is  particularly  important 
in  the  case  of  the  tricycle  nose  wheel. 

There  is  no  rational  theory  of  this  dangerous  phenomenon,  although 
some  tests  have  been  made  in  the  U.  S.  A.  and  in  Germany. 

The  present  report  deals  with  a  rather  simple  theory  that  agrees 
with  the  known  tests  and  enables  dimensions  of  the  land-gear  to  be  com- 
puted so  as  to  avoid  shimmy  without  resorting  to  dampers.   Tests  with  a 
full-scale  landing  gear  are  described. 


I.  THE  PROBLEM  OF  SHIMMY 


When  a  wheel  fitted  with  a  tire  is  designed  to  pivot  freely  about 
a  vertical  pivot  AA *  (fig.  1)  and  this  pivot  is  given  a  horizontal 
forward  motion  while  the  wheel  is  made  to  roll  on  the  ground,  it  happens 
that  it  spontaneously  assumes  a  self-sustained  oscillating  motion  about 
the  pivot  AA'.   This  phenomenon  popularly  termed  shimmy  or  wiggle  is 
frequently  observed  on  the  tail  wheel  of  a  conventional  airplane  and 
on  the  nose  wheel  of  tricycle  landing  gears.   In  the  latter  case  it  can 
become  extremely  violent  and  consequently  very  annoying,  and  occasionally 
induce  failure  of  the  landing  gear.   In  fact,  it  has  been  called  the 
"bete  noire"  (bugbear)  of  the  tricycle  landing  gear:   "the  problem  of 
nose  wheel  shimmy,  long  the  bugaboo  of  tricycle  landing  gear  .  .  .  ." 
(Aerodigest,  March  15,  1°M,  p.  13^) . 

The  Americans,  promoters  of  the  tricycle  landing  gear,  have  already 
made  large-scale  theoretical  and  experimental  studies  of  this  phenomenon, 
but  as  yet,  there  seems  to  be  no  complete  and  acceptable  theory.   The 
writer  proposes  to  supply  this  theory  in  the  present  report.   The  sug- 
gested explanation  involves  only  the  elementary  mechanical  properties 
of  the  tires  which  are  all  well  defined  and  easily  measurable. 

"Etude  Theorique  du  Shimmy  des  Roues  d  "Avion",  Office  National 
d 'Etudes  et  de  Recherches  Ae'ronautiques ,  Publication  No.  7 ,    19^-8. 
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Before  any  particular  hypothesis,  it  is,  in  fact,  quite  evident 
that  shimmy  can  have  no  cause  other  than  the  reactions  exerted  hy  the 
ground  on  the  vheel-reactions  transmitted  by  the  tire.   These  reactions 
can  only  be  contact  reactions,  due  to  the  adherence  of  the  tire;  they 
are  necessarily  reduced  to  a  force  and  to  a  torque.   To  analyze  these 
reactions  more  correctly  it  is  necessary,  first  of  all,  to  form  the 
theory  of  the  elementary  mechanical  properties  of  the  tire. 


II.  THE  MECHANICS  OF  THE  TIRE 

A  tire  is  not  merely  an  elastic  element  possessing  vertical  elas- 
ticity which  is  essentially  its  reason  for  existence,  that  is,  a  simple 
vertical  spring  placed  between  ground  and  wheel;  it  is  a  complex 
mechanical  unit  having  a  certain  number  of  other  elementary  mechanical 
properties.   For  the  ensuing  theory  the  four  following  fundamental  fac- 
tors are  considered: 

1.  Lateral  Elasticity 

2.  Torsional  Elasticity 

3.  Drifting 
k.  Turning 

1.  Lateral  Elasticity 

At  rest,  under  the  action  of  a  relatively  moderate  force  F,  applied 
parallel  to  the  wheel,  that  is,  perpendicular  to  its  plane,  the  wheel  is 
displaced  with  respect  to  the  ground,  and  through  the  simple  effect  of 
the  deformation  of  the  tire,  by  an  amount  A  proportional  and  parallel 
to  the  force  F.   This  displacement  is  due,  as  shown  in  figure  2,  to  the 
simple  deformation  of  the  tire  near  the  surface  of  contact  with  the 
ground,  without  any  modification  of  that  contact  surface  as  long  as  the 
force  F  does  not  exceed  the  limit  of  adherence. 

In  the  formula 

A  =  TF 

T      is   the   coefficient  of  lateral  elasticity  of  the   tire. 

2.  Torsional  Elasticity 

At  rest,  under  the  action  of  a  relatively  moderate  moment  torque  M 
and  vertical  axis  applied  to  the  wheel,  the  latter  is  subjected  to  a 
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rotation  of  angle  a  proportional  to  the  moment  M.   As  in  the  preceding 
case,  the  mechanism  of  this  rotation  is  easy  to  understand.   It  results 
from  the  deformation  of  the  tire  near  the  surface  of  contact  with  the 
ground,  this  surface  itself  remaining  unaltered  as  a  result  of  the 
adherence. 

In  the  formula 

a  =  SM 

S   is  called  the  coefficient  of  torsional  elasticity  of  the  tire. 

3.  Drifting 

Under  the  action  of  the  side  force  F  a  wheel  mounted  with  a  tire 
is,  according  to  the  foregoing,  subjected,  first,  to  the  lateral  dis- 
placement A  =  TF,  that  is,  its  center  is  shifted  from  0  to  0 '  and  its 
plane  from  P  to  P '  (fig.  3).   In  addition,  if  the  wheel  is  made  to  roll, 
it  is  seen  that  its  track  is  no  longer  contained  in  the  plane  of  the 
wheel,  it  moves  without  sliding  along  a  straight  line  O'M,  forming  an 
angle  6  proportional  to  the  force  F. 

In  the  formula 

5  =  DF 

D   is  termed  the  drift  coefficient  of  the  tire.   This  phenomenon  is  in 
effect,  a  drift  comparable  to  the  drift  of  a  ship  in  side  wind.   It  is 
particularly  evident  on  a  bicycle,  or  better  yet,  on  a  motorcycle  with 
an  insufficiently  inflated  front  wheel.   It  is  this  deflection  of  the 
rubber  tires  that  enables  mothers  to  steer  with  ease  baby  carriages  with 
non-swiveling  wheels  without  being  forced  to  make  the  wheels  skid  on 
the  ground  or  to  raise  two  of  the  wheels  in  order  to  turn  the  carriage. 

A  remarkable  thing,  this  drift  is  produced  without  the  least  sliding 
of  the  tire  over  the  ground.   It,  as  well  as  the  other  properties  of 
the  tire,  can  be  explained  by  a  very  simple  diagram.   Supposing  the 
tire  is  equivalent  to  an  infinite  number  of  small  coiled  springs  mounted 
radially  on  the  periphery  of  the  wheel  and  endowed  with  elasticity  in 
the  radial  and  the  lateral  directions  (fig.  h) .      This  evidently  repre- 
sents the  vertical  and  lateral  elasticity  of  the  tire,  as  well  as  the 
torsional  elasticity.   It  also  explains  the  drift.   In  effect: 

Suppose  that  the  wheel  is  at  rest;  a  certain  number  of  springs  n 
are  in  contact  with  the  ground,  as  outlined  in  figure  5.   Now,  when  a 
side  force  F  is  applied  to  the  wheel,  it  first  is  subject  to  a  lateral 
displacement  A  =  TF  as  a  result  of  the  simple  lateral  deformation  of 
the   n  springs  in  contact  with  the  ground  without  the  least  sliding 
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of  their  contacts  with  the  ground.   As  a  consequence  each  of  these 
springs  assumes  a  certain  lateral  tension,  the  sum  of  these  n  tensions 
balancing  the  force  F.   Now  the  wheel  is  made  to  roll.   Obviously,  on 
account  of  the  displacement  A  the  contact  of  the  spring  n  +  1  with 
the  ground  is  no  longer  in  the  alignment  of  the  first  n's,  but  dis- 
placed laterally  by  an  amount  equal  to  A.   On  the  other  hand,  at  the 
moment  that  the  spring  n  +  1  makes  contact  with  the  ground,  the 
spring  1  loses  this  contact;  but  this  spring  has  a  certain  lateral 
tension  which  suddenly  stops  as  soon  as  the  spring  leaves  the  ground, 
while  at  the  instant  of  making  contact  with  the  ground  the  spring  n  +  1 
has  as  yet  no  tension.   The  sum  of  the  reactions  of  the  ground  is  there- 
fore diminished  by  the  tension  of  spring  1  and  it  no  longer  balances 
the  force  F.   To  recover  this  reduction  the  wheel  is  again  subjected  to 
a  slight  lateral  displacement  so  that  the  point  n  +  2   is  again  dis- 
placed or  forced  off  laterally  with  respect  to  the  point  n  +  1.   The 
process  is  repeated  every  time  one  spring  leaves  the  ground  and  another 
makes  contact.   It  is  seen  that  the  wheel  rapidly  assumes  the  state  of 
a  lateral  displacement  proportional  to  the  path  traveled  in  the  direc- 
tion of  rolling  of  the  wheel.   The  normal  state  can  be  considered  as 
reached  as  soon  as  the  point  n  has  left  the  ground,  that  is,  as  soon 
as  the  wheel  has  covered  a  distance  equal  to  the  length  I     of  the 
surface  of  contact.   This  length  being  slight,  it  all  happens  as  if 
that  state  was  instantaneous.   Since  each  one  of  the  elementary  dis- 
placements involved  in  this  analysis  is  proportional  to  F,  the  same 
applies  to  their  sum,  hence  the  formula 

drift  =  DFx 

x  =  distance  traveled  in  direction  of  rolling,  and  consequently 

angle  of  drift  5  =  DF 

Thus  the  drift  is  explained  by  the  natural  elasticity  of  the  tire 
without  involving  the  least  effect  of  skidding  or  sliding  of  the  tire 
over  the  ground.   This  phenomenon  is  almost  the  same  as  the  longitudinal 
creep  of  transmission  belts;  it  is  fully  comparable  with  pseudo-slipping 
in  which  the  drive  wheels  of  a  vehicle  make  a  number  of  revolutions 
greater  than  that  resulting  from  the  distance  traveled  by  the  vehicle 
and  greater  than  that  of  nondriving  wheels  of  the  same  diameter  carrying 
the  same  load.   Pseudo-slipping  which  has  been  the  subject  of  many 
experimental  studies  for  locomotives,  is  explainable  in  the  same  manner 
without  involving  real  slipping  between  wheel  and  rail;  it  results  from 
elastic  deformations  of  the  steel  near  the  surface  of  contact  and  its 
simplified  theory  is  the  same  as  that  of  drift.   Pseudo-slipping  is  a 
longitudinal  drift  and  the  drift  a  lateral  pseudo-slipping.   Likewise, 
braking  also  produces  negative  pseudo-slipping. 
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Within  certain  limits  the  drifting  of  a  wheel  can  be  compared  to 
the  deformation  of  a  solid  body.   The  displacement  A  plays  a  part 
similar  to  that  of  the  reversible  elastic  deformation,  and  drift  DFx 
similar  to  that  of  the  irreversible  permanent  deformation  arising  from 
a  viscous  creep,  the  path  covered  x  in  the  drifting  playing  the  part 
of  the  time  in  creep. 

The  torque  accompanying  the  drift.-  The  preceding  analysis  indi- 
cates that  the  drift  is  necessarily  accompanied  by  a  very  important 
mechanical  phenomenon,  namely,  the  displacement  of  the  lateral  reaction 
of  the  ground  toward  the  rear  of  the  wheel.   Return  to  figure  h   and 
consider  the  wheel  at  the  moment  the  spring  k  +  n  comes  in  contact 
with  the  ground.   On  account  of  the  drift  the  points  k  +  1,  k  +  2,  .  .  . 
k  +  n,  are  alined,  conformably  as  in  figure  6,  along  a  straight  line 
forming  an  angle  5  with  the  plane  OP  of  the  wheel.   On  the  other  hand, 
at  the  instant  point  k  +  n  makes  contact  with  the  ground,  the  plane 
of  the  wheel  obviously  passes  through  the  point  k  +  n,  which  is  termed 
the  head  point  of  contact  of  the  wheel  with  the  ground.   Since  the  ten- 
sion is  proportional  to  the  distance  between  ground  contact  and  wheel 
plane  the  result  is  a  lateral  tension  of  some  of  the  n  springs  in 
contact  with  the  ground  proportional  to  the  distance  between  ground 
contact  and  head  point.   This  tension  varies  therefore  from  zero  for 
head  point  k  +  n  to  a  maximum  value  for  the  tail  point  k  +  1.   The 
reaction  F'  of  the  ground,  the  resultant  of  these  different  tensions, 
therefore  passes  through  the  point  G  located  at  two-thirds  of  the 
line  of  contact  from  the  head  point.   The  ground  reaction  F'   therefore 
does  not  pass  through  the  center  M  of  the  line  of  contact.   While 
parallel,  equal  and  of  opposite  sign  of  F,  the  reaction  F'   is  there- 
fore not  directly  opposite  to  the  force  F  which,  passing  through  the 
projection  0  of  the  center  of  the  wheel,  passes  also  through  M.   The 
result  is  a  torque  C  about  the  vertical  axis  of  the  wheel. 

With  e   denoting  the  distance  MG  and  noting  that  the  angle  of 
drift  S   is  practically  always  small,  this  torque  has  the  value 

C  =  eF 

In  the  case  of  the  wheel  shown  in  figure  h   it  is  readily  apparent 

that   e  =  MG  =  -:,   where  I   =  length  of  contact,  that  is,  distance  between 
0 

head  point  and  tail  point.   But  it  is  fitting  to  note  at  the  same  time 
that,  while  the  scheme  is  helpful  for  understanding  the  drift  and  the 
phenomenon  of  the  accompanying  torque,  the  tire  is  nevertheless  a  more 
complicated  mechanism.   The  simplified  scheme  supposes  that  the  small 
radial  springs,  equivalent  to  a  tire,  are  independent.   But  owing  to  the 
continuity  of  the  pneumatic  tire,  these  springs  must  be  considered 
strongly  tied  elastically  to  one  another  in  such  a  way  that  the  deforma- 
tion of  one  of  them  also  involves  almost  as  much  deformation  of  the  two 
adjacent  ones. 
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However,  it  can  be  shown  that  between   e  and   Z  the  important 
relation 

l<e<i 
6  2 

exists. 

For  the  present,  without  stopping  to  demonstrate  this  formula,  the 
following  are  assumed  experimental  facts: 

(1)  Every  relatively  moderate  lateral  force  F  applied  to  a  tire 
in  motion  produces  a  drift  proportional  to  the  angle  5  =  BF. 

(2)  This  drift  is  accompanied  by  a  couple  C  =  eF  tending  to 
orient  the  wheel  in  the  true  direction  of  displacement,  that  is  to 
reduce  the  angle  of  drift  spontaneously. 

The  coefficient   e  for  a  given  load  supported  by  the  wheel  is  a 
characteristic  constant  of  the  tire  which  is  termed  length  of  displace- 
ment; and  torque  Fe  the  torque  accompanying  the  drift. 

k.    Turning 

Under  the  action  of  a  couple  of  relatively  moderate  torque  M 
about  the  vertical  axis,  applied  to  the  center  of  the  wheel,  the  latter 
is  subjected,  first,  to  a  static  rotation  of  angle  cc,  as  shown  in  the 
foregoing.   But,  when  the  wheel  is  made  to  roll,  it  no  longer  moves 
along  a  straight  track,  but  along  a  circular  trajectory  of  radius  p, 
according  to  figure  f.      The  arc  of  circle  AB  is  such  that  its  curva- 
ture is  proportional  to  the  moment  M.   Hence 

I  =  RM 
P 

that  is,  for  a  trajectory  AB  of  length  s,  the  wheel  has  turned  through 
an  angle 

3  =  RMs 

R  is  termed  the  coefficient  of  turning  of  the  tire. 

To  a  certain  degree  the  rotation  of  the  wheel  can  be  compared  with 
the  deformation  of  a  solid  body.   The  angle  a  then  plays  a  part  simi- 
lar to  that  of  the  reversible  elastic  deformation,  and  angle  3  similar 
to  that  of  the  irreversible  permanent  deformation. 
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It  is  seen  that  the  moment  M  necessary  to  make  the  wheel  turn 
to  a  permanent  angle  3   is  inversely  proportional  to  the  traveled 
distance  s.   It  is  this  property  of  turning  that  explains  the  ease  of 
steering  the  wheels  of  an  automobile  in  operation  while  turning  the 
wheels  of  an  automobile  when  standing  still  is  much 'more  difficult 
because  it  can  only  be  achieved  by  an  entirely  different  mechanism, 
namely,  by  exceeding  the  limit  of  adhesion  of  the  tire  with  the  ground. 
If  this  adhesion  were  perfect,  that  is,  if  no  sliding  were  possible, 
the  permanent  turning  of  the  wheel  at  rest  would  be  impossible,  whereas 
the  sliding  of  the  wheel  plays  no  part  in  steering  the  wheel  when 
rolling,  as  will  be  proved  elsewhere  in  a  more  rigorous  analysis  of 
this  phenomenon.   As  for  the  drift,  turning  can  be  explained  on  the 
basis  of  figure  k;    the  process  is  exactly  the  same;  simply  replace  the 
lateral  displacements  by  rotations.   Turning  is  to  the  torsional 
elasticity  what  deviation  is  to  the  lateral  elasticity.   The  reality 
of  turning  is  assumed  as  an  experimentally  established  fact.   As 
regards  the  proportionality  of  rotation  3  to  moment  M   (like  the 
proportionality  of  drift  to  side  force)  from  the  moment  the  existence 
of  such  rotation  is  assumed  the  proportionality  must  be  assumed  proved 
from  the  following  very  general  reasoning:   The  tire  is  a  complex 
elastic  system,  whatever  the  effects  of  the  forces,  they  must  be  pro- 
portional, at  least  as  long  as  the  forces  do  not  exceed  certain  limits. 

The  properties  of  turning  and  drift  in  question  indicate  that  a 
wheel  fitted  with  a  pneumatic  tire,  or  more  general,  any  wheel  fitted 
with  any  tire  has  the  important  property  of  being  able  to  roll  while 
making  a  certain  angle  with  its  trajectory.   This  is  the  essential 
phenomenon  which,  as  will  be  shown,  underlies  and  basically  explains 
the  shimmy  of  pneumatic  tires. 

The  five  characteristic  tire  factors  produced  by  the  foregoing 
analysis  are 

(l)  Coefficient  of  lateral  elasticity  T  defined  by  the  equation 

A  =  TF 
,,    (2)  Coefficient  of  torsional  elasticity  S  defined  by  equation 

a  =  SM 
(3)  Coefficient  of  drift  D  defined  by  equation 

6  =  DF 
(k)   Coefficient  of  turn  R  defined  by  equation 

3  =  RMs 
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(5)  Length  of  displacement  e   defined  by  equation 

C  =  eF 

III.  MATHEMATICAL  THEORY  OF  ELEMENTARY  SHIMMY 


In  order  to  obtain  a  clear  picture  of  the  mechanism  of  this 
phenomenon,  the  number  of  parameters  are  reduced  by  starting  the  study 
with  what  may  be  called  elementary  shimmy  or  shimmy  with  one  degree  of 
freedom,  that  is,  the  shimmy  obtained  when  the  spindle  axis  of  the 
wheel  is  assumed  to  be  impelled  only  by  a  straight,  uniform  horizontal 
motion.   This  is  equivalent  to  assuming  the  pivot  absolutely  rigid  and 
the  mass  of  the  airplane  very  great  with  respect  to  that  of  the  wheel, 
so  that  the  effects  of  wheel  reactions  on  the  trajectory  of  the  caster 
axis  can  be  disregarded.   The  position  of  the  wheel  with  respect  to 
the  airplane  is  thus  defined  by  one  parameter,  the  angle  9     of  its 
plane  with  the  axis  of  the  airplane  (fig.  8). 

The  study  of  elementary  shimmy  is  divided  in  two  parts. 

In  the  first  part  is  assumed  that  the  reactions  of  torsion  and 
turn  are  negligible  compared  with  the  reactions  of  lateral  elasticity 
and  drift,  that  is  that  the  effects  of  torque  M  are  negligible  com- 
pared with  those  of  force  F,  which  is  the  same  as  assuming  that  the 
coefficient  R  is  sufficiently  great  with  respect  to  coefficients  T 
and  D.   Later  on  it  shall  be  shown  that  this  is  practically  the  case 
for  ordinary  tires.   It  produces  a  simplified  theory  of  elementary 
shimmy.   The  second  part  contains  th£  complete  theory  of  elementary 
shimmy. 


A.  Simplified  Theory  of  Elementary  Shimmy 

It  is  assumed  that  the  pivot  is  vertical  and  the  pivoting  free, 
that  is  without  restoring  torque  and  without  braking. 

Let  P  denote  the  track  on  the  ground  of  the  pivoting  axis  at 
instant  t,   x  the  distance  traveled  by  this  point,   0  the  projection 
of  the  wheel  center,  9     the  angle  of  the  line  PO  with  the  trajectory  PX 
of  the  point  P,  that  is  with  the  axis  of  the  airplane,  AMGB  the  axis 
of  symmetry  of  the  contact  surface  of  the  tire  with  the  ground,  A  being 
the  head  point,   B  the  tail  point,  M  the  middle  of  AB  and  G  the 
center  of  the  ground  reactions,  that  is  the  point  so  that  MG  =  6 ,   y 
and  z  the  distances  from  PX  of  points  M  and  0  and  finally 
that  \|/   is  the  angle  of  the  line  of  contact  AB  with  PX. 
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Considering  now  the  place  of  point  M  on  the  ground  whose  tangent 
is  the  straight  line  AB,  since  the  point  M  passes  successively 
throught  points  B  and  A  due  to  the  rolling  of  the  wheel.   Hence 
the  formula  (fig.  8) 

dx   y 


It  is  readily  seen  that  the  angle  of  drift  5   is  equal  to  9  +  ty,   and 


since  9  =  — , 


hence 


dy  +  z 
dx   a 


Now,  if  F   is  the  resultant  of  the  ground  reactions,  the  pre- 
liminary study  made  on  drift  makes  it  possible  to  write 

5  =  DF 


hence  the  first  equation  reads 


^  +  £  =  DF 
dx   a 


On  the  other  hand,  the  action  of  the  lateral  elasticity  yields  the 
second  equation 

z  -  y  =  TF 


Lastly,  the  force  F  being  applied  at  G,  that  is  a  distance  a  +  e 
from  point  P,  produces  a  torque  equal  to  -F (a  +  s)  with  respect  to 
this  point.   In  other  words,  to  the  primary  torque  -Fa  the  accompanying 
torque  -Fe  must  be  added.   So,  if  I  signifies  the  inertia  of  the 
whole  oscillating  system  with  respect  to  the  pivot 
tion  reads 


P,  the  third  equa- 


dfe 
dt2 


-F(a  +  e) 


Noting  that  in  the  first  equation 


dy  =  dy  dt 
dx   dt  dx 


1  dy 
v  dt 
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where  v  =  speed  of  airplane,  and  that  in  the  third  equation 

9   =  £ 
a 

the  system 

±  *y  +  *  =  df  (i) 

v  dt   a 
z  -  y  =  TF  (2) 


Id£|=_F(a+e)  (3) 

a  dt2 


is  obtained.   The  elimination  of  F  from  equations  (l)  and  (3)  leaves 


1  dy   z  ,     ID    d2 


+  T  +   ,     .  ^  =  o  M 


v  dt   a   a(a  +  g)  dt2 

On  the  other  hand,  the  elimination  of  F  from  equations  (2)  and  (3) 
leaves 

'-•♦-r35-^  (5) 

a(a  +  e)  dt2 

*. 

Substitution  of  this  value  of  y  into  the  preceding  equation  finally 
gives  the  differential  equation  of  the  third  order 

■  +  ^^lfz  +  a  dz  +  z  =  Q  (6) 


v(a  +  e)  dt3   a  +  e  dt2   v  dt 


which  governs  the  motion  of  the  wheel  center.   The  general  solution  of 
this  equation  is 


s^t      S2t      sot 
z  =  Cxe    +  C2e    +  036 


* 
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where  Ci,  C2,  and  C3  are  arbitrary  parameters  and  s-j_,  sp,  and  So  are 
the  roots  of  the  characteristic  equation 

£2 s3  +  -2^-   s2  +  2:  3  +  1  =  0-  (7) 

v(a  +  e )      a  +  6      v 

Now,  while  the  mathematical  representation  of  the  wheel  center  can  be 
expressed  by  a  unique  formula,  it  is  well  to  bear  in  mind  that  the 
physical  appearances  of  this  motion  are  entirely  different  depending 
upon  the  positive  or  negative,  real  or  imaginary  values  of  the  charac- 
teristic equation.   Hence  the  necessity  for  discussing  the  values  of 
s1j  s2>   an^  s3   i-n  terms  of  the  coefficients  of  the  characteristic 
equation,  and  as  a  consequence,  as  function  of  six  parameters  T,  D, 
e,  a,  I,  v  which  characterize  the  tire,  the  wheel  and  the  airplane 
speed. 

It  should  be  noted  that  the  characteristic  equation  has  no  posi- 
tive coefficients,  hence  can  have  no  positive  root;  so  if  this  equation 
has  real  roots,  they  must  naturally  be  negative.   On  the  other  hand, 
the  equation  has  always  at  least  one  real  root  since  its  first  member 
varies  continuously  from  -00  to  +°°. 

Suppose  s-^  =  -a  is  this  root. 

From  the  point  of  view  of  real  roots  there  can  be  only  two  possible 
cases : 

(l)  Three  negative  real  roots. 

The  general  solution  is  then  of  the  form 

-at      -pt      -7t 
z  =  C]_e    +  C2e    +  Coe 


a,   3,  and  7  being  positive  numbers.   Whatever  the  initial  conditions 
that  cause  the  variations  of  C-^,  Cg,  Co,  the  motion  of  the  wheel  is, 

in  this  case,  an  aperiodic  and  convergent  oscillation. 

(2)  One  negative  real  root  S]_  =  -a  and  two  conjugate  imaginary  roots. 
It  is  known  that,  when  s2  and  So  are  conjugate  imaginary,  of 

Sot  Sot 

the  form  X,  ±  coi,  the  two  terms  Coe     and  Coe  ->   combine  to  give 
a  term  of  the  form 
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Be^tsin(aDt  +  cp) 

cd  being  a  number  essentially  positive  and  X.  a  positive  or  negative 
real  number. 

The  general  solution  can  then  be  written  as 

z  =  Ae_<rt;  +  Be^tsin(cut  +  cp) 

the  three  constants  of  integration  which  depend  upon  the  initial  condi- 
tions of  motion,  are  the  parameters  A,  B,  cp. 

From  the  point  of  view  of  the  physical  appearances  of  the  motion 
this  case  is  split  in  two: 

(a)  If  X  <  0  the  oscillations  are  convergent  periodic 

(b)  If  X  >   0  the  oscillations  are  divergent  periodic 

From  among  the  possible  forms  of  the  solution  z,  the  only  case 
in  which  there  is  a  divergence  of  oscillation  is  therefore  X  >   0. 

Now  an  attempt  will  be  made  to  calculate  the  angular  frequency  cx> 
and  the  divergence  X     of  the  oscillations  in  terms  of  the  coefficients 
of  the  equation  of  motion  of  the  wheel,  that  is  in  terms  of  the  charac- 
teristic coefficients  of  the  tire  and  the  wheel.   The  condition  of 
divergence  of  the  oscillations,  that  is  the  condition  of  instability 
is  derived  by  writing  A.  >  0.   In  all  other  cases  there  is  convergence, 
that  is  stability. 

Calculation  of  divergence  X.  and  frequency  cd 

As  X     and  cd  are  independent  of  A,  B,  cp,  it  is  assumed  that 
A  =  0,   B  =  1,   cp  =  0.   Hence 

z   =   e^sin  cut 


dz   =  ext 


A,  sin  cut  +  cd  cos  cut 
dt  *-  -1 


d2 


dt2 


z  =  e\t 


|(\2   _  cD2)sin  CDt  +   2  Xca  cos  CDt 


dJz  xt 


|x,(X.2   -    3cD2)sin  CDt  +  cd(3X2   -  cd2)cos   cut 
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which,  written  into  the  equation  of  motion  of  the  wheel  gives  an 
equation  of  the  form 


»*•*  jg(X.,cjo)sin  cut  +  h(X,co)cos  cut]  =  0 


This  equation,  which  must  be  checked  whatever  t  may  be,  gives  the 
system 


g(X-,Co)  =  0 

h(X,co)  =  0 


of  two  equations  that  define  X     and  cu.   Hence  the  system 


IT 


v(a  +  e) 


\(XT  -  3o>2)  + 


JL_(X2  _  co2)  +  £  X   +  1  =  0 
a  +  e  v 


IT 


v(a  +  e) 


(3X2  -  cd2)  + 


^.  +  *  =  o 

a  +  e     v 


which  can  also  be  written 


2T\(\2  +  cu2)  +  vD(\2  +  co2)  _  v  (a  +  e)  =  0 

I 


T(3X2  -  co2)  +  2vDX  +  a(a  +  €)  =  0 

I 


(8) 


(9) 


The  second  of  these  equations  gives 


0,2  =  3X2  +  2vm  +  a(a  +  e) 
T        IT 


(10) 


which,   entered   in  equation    (8),   and     co    eliminated,    gives   the   equation 
8T2IX3  +   8DTIvX2  +  2[D2Iv2  +  a(a  +  e)TJx  +    (aD  -  T) (a  +   6 )v  =  0   (ll) 
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and  defines  the  values  of  X   in  terms  of  the  velocity  v  and  the 
characteristic  coefficients  of  the  tire  and  the  wheel.   With  regard  to 
the  abscissa  v  and  ordinate  X,  this  equation  represents  a  cubic 
passing  through  the  origin,  symmetrical  with  respect  to  the  origin  and 
asymptotic  to  the  axes  of  the  abscisses.   It  is  easy  to  plot  this  cubic 
and  to  identify  the  characteristic  elements. 

Plotting  of  curve  X(v).-  The  preceding  equation  is  of  the  second 
degree  in  v  and  can  be  written 

2D2IXv2  +  [SDTIX2  +  (aD  -  T)(a  +  e)j  v  +  8T2IX3  +  2a(a  +  e)TX  =  0  (12) 

Hence  the  curve  of  the  values  of  v  is  easily  plotted  against  X.   It 
is  readily  apparent  that  for  each  given  real  value  of  X  there  always 
are  two  real  values  of  v,  the  only  condition  being  that  the  discriminant 
be  positive,  that  is 

l6DT2IX2  <  (a  +  e)(aD  -  T)2 
or 


|aD_^|  /a_^=  - 
14-T    V  DI      ' 


Maximum  or  minimum  A, .  -  When 

2   (aD  -  T)2  a  +  e 


X< 


l6T2     DI 


the  curve  X(v)   passes  therefore  through  a  maximum  or  minimum  M  and 
the  value  vm  of  the  corresponding  speed  is  such  that 

p   &T2IX2  +  2a(a  +  e)T 

vm  =  

2D2I 

we  get,  by  replacing  X  by  its  value 

aD 


vm 


mfH  (12a) 

2D   V  DI 
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Tangent  at  the  origin. -  This  is  the  straight  line  of  the  equation 

2aTX  +  (aD  -  T)v  =  0 

hence  its  slope  is 

T  -  aD 
2aT 

Representative  curve.-  In  consequence  of  this  value  of  the  slope 
at  the  origin,  the  representative  curve  can  assume  different  shapes, 
depending  upon  whether  aD  >  T  or  aD  <  T. 

First  case:   aD  >  T 

The  curve  X(v)  has  the  shape  represented  in  figure  9«   The 
velocity  v  being  necessarily  positive,  it  is  seen  that  X   is  always 
negative. 

Second  case:   aD  <  T 

The  curve  X.(v)  has  the  form  represented  by  figure  10. 

The  speed  v  being  necessarily  positive,  X      is,  in  this  case, 
always  positive.   The  only  condition  for  divergence  of  oscillations  of 
the  wheel  is,  as  was  shovm  before,  that  X   >  0,  and  the  only  condition 
of  convergence  is  X.  <  0.   In  consequence  the  condition  of  divergence 
or  instability  is  reduced  to 

aD  <  T 

and  that  of  convergence  or  stability  to 

aD  >  T 

This  result  could  have  been  reached  more  rapidly  by  application 

of  Routh's  general  rules  (or  the  equivalent  method  by  Hurwitz)  of  which 

the  particular  application  to  an  equation  of  the  third  order  with  posi- 
tive coefficients  of  the  form 


d3z  ,  _  d2z 


dz 


A  ^-f  +  B  2-2-     +  C  —  +  Dz  =  0 
dt3     dt2      dt 


reduces  to  the  following  rule:   the  necessary  and  sufficient  condition 
of  convergence  is  DC  >  AD. 
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In  fact,  in  the  case  of  equation  (6)  it  is  seen  that  this  rule 
yields  immediately  aD  >  T  as  the  condition  of  convergence.   But  this 
method  involves  the  use  of  the  original  works  by  Routh  ("On  the  Stability 
of  a  Given  State  of  Motion,"  Adams  Prize  Essay,  1877 )  and  his  treatise 
("Advanced  Rigid  Dynamics")  or  the  works  of  Hurwitz. 

Moreover,  it  supplies  no  accurate  data  about  the  variations  of 
X      in  terms  of  the  airplane  speed  v. 

Resuming  the  preceding  discussion  it  is  seen  that  when  aD  <  T, 
the  oscillations  are  unstable,  that  is  the  wheel,  while  rolling, 
assumes  a  motion  of  divergent  oscillations. 

In  practice  the  aforementioned  conclusions  are  not  vigorously 
checked  except  for  the  start  of  the  phenomenon.   This  is  due  to  the 
fact  that  our  equations  are  applicable  only  to  small  deformations  and 
cease  to  be  rigorous  when  the  oscillations  attain  sizeable  amplitude. 
So,  in  the  initial  phase  there  is  a  phenomenon  of  unstable  oscillations, 
exhibited  by  a  divergence  which  can  become  very  accentuated  when  the 
velocity  is  near  that  of  the  previously  defined  speed  vm.   But  the 
motion  is  never  infinitely  divergent.   As  soon  as  the  oscillations 
have  reached  sufficient  amplitude,  a  state  of  continuous  oscillatory 
motion  is  established. 

The  process  was  frequently  believed  to  be  a  phenomenon  of  resonance, 
arising  from  a  lack  of  symmetry  of  the  tire  or  the  wheel  or  similar 
causes,  and  interpreted  as  such  in  most  theories  of  front  wheel  shimmy 
in  automobiles.   This  explanation  cannot  be  sustained  in  the  face  of 
the  severe  cases  of  shimmy  over  the  smoothest  of  grounds  with  perfectly 
balanced  wheels.   Furthermore,  the  violence  of  the  phenomena  in  some 
of  its  manifestations  and  the  extent  of  the  speed  range  in  which  it  can 
manifest  itself,  are  enough  to  prove  that  no  resonance  phenomenon  is 
involved  but  rather  a  phenomenon  of  vibratory  instability.   The  present 
theory  shows  that  shimmy  can  be  explained  in  the  most  elementary  way 
by  the  inherent  mechanical  properties  of  the  tire. 

The  simplicity  of  the  condition  of  convergence  aD  >  T   is  sur- 
prising.  In  particular,  it  is  extremely  unusual  to  find  that  neither 
the  wheel  inertia  I,  nor  the  offset  e   or  the  velocity  v  figure  in 
the  condition  of  convergence.   This  condition  gives  a  rule  of  extreme 
simplicity;  it  is  sufficient  to  follow  this  rule  in  the  construction  of 

T 
the  undercarriage  a  >  =r  in  order  to  eliminate  elementary  shimmy. 

Third  intermediary  case:   aD  =  T 

In  this  case,  equation  (ll)  gives  immediately  X  =  0. 
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On  the  other,  equation  (10)  gives 


a^  = 


a(a  +  e  )   a  + 


IT       DI 
Therefore  the  general  solution  takes  the  form 

z  =  Ae~at  +  B  sin(a>t  +  9) 
The  term  Ae_a   rapidly  approaches  zero,  leaving  only  the  second  term 


/a  +  £ 
which  represents  a  sinusoidal  motion  whose  frequency  od  =  1/ is 

not  dependent  on  the  airplane  speed  v. 

Frequency  u) 

The  quantity  X.  can  be  eliminated  in  the  same  manner  as  gd  from 
equations  (8)  and  (9)-   Multiplying  equation  (8)  by  -3  and  equation  (9) 
by  2X.  and  then  adding  up,  gives  the  equation 

vIDX2  +   2  {a (a  +  e)  -  UlTa£]  X.  +  3v  \{a.  +   e  )  -  IDoj2J  =  0 
On  the  other  hand,  equation  (9)  can  also  be  written 

3ITX.2  +  2vIDX  +  [a  (a  +  e)  -  iTuFj   =  0 

thus  giving  two  equations  of  the  second  degree  in  X..   It  is  known  that 
the  result  of  the  elimination  of  x  from  the  two  equations 

p 
ax  +  bx  +  c  =  0 

a  *x2  +  b  'x  +  c  '  =  0 
is  the  relation 

(acf  -  ca')2  -  (ab'  -  ba ' ) (be  '  -  cb ' )  =  0 

obtained  by  equating  the  result  of  these  two  equations  to  zero. 
Applying  this  formula  to  the  two  preceding  equations  in  A.,  gives  the 
e  quat  ion 
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v2I  JD[a(a  +   e)    -    ITcd2]    _   c/r  [[a  +   e )    -   IDcd2]  >      =  1+  JV2ID2   -   3T[a(a  +  €  ) 
4lTu£]>   <[a(a  +   e)    -   UlTa£|  [a (a  +   e)    -   ITcd2]    -   3v2ID  [(a  +  e)   -   IDcd2]  f 

(13) 

which  defines  the  values  of  cd  in  terms  of  v  and  the  characteristic 
coefficients  of  the  tire  and  wheel.   This  equation  of  the  sixth  degree 
in  v  and  cd  is  reduced  to  one  of  the  third  degree  by  a  change  of 
the  coordinates   v2  =  x  and  to2  =  y.   Unfortunately  the  direct  study 
of  this  equation  proved  to  be  rather  difficult,  and  was  therefore 
abandoned  in  favor  of  a  more  simple  method  of  defining  the  variations 
of  cd  in  terms  of  v.   In  fact  it  can  be  stated  that  the  writer  has 
studied  the  variations  of  \     in  terms  of  v,  and  was  thus  able  to 
study  the  variations  of  co  in  terms  of  \,   and  the  problem  of  the 
variations  of  the  frequency  cd  with  respect  to  v  can  be  solved  the 
same  way. 

The  velocity  v  is  easily  eliminated  from  equations  (8)  and  (9); 
hence  the  equation 

{X2  +  ^f  +   a_W3  _  a\  2  _  a^/1   a\  2  +  a(a  +  Q2  =    {lk) 
I  \d      i)  I  Vd   T/      z2dt 

which  defines  the  values  of  to  in  terms  of  \     and  the  characteristic 
coefficients  of  the  tire  and  wheel.   This  equation  of  the  fourth  degree 
in  X     and  cd  is  reduced  to  one  of  the  second  degree  by  changing  the 
coordinates  X.  =  X  and  cd2  =  Y.   The  result  is 

(x  +  y)2  f  !L±±(1  _  ±)x  _  SL+-«/i  +  aV  +  *(*  +  O2  =  0       (15) 

I      \B       T/  I      \P        T/  Z2DT 


The  equation  of  the  group  of  asymptotic  directions 

(X  +  Y)2  =  0 

indicates  that  the  conic  involved  is  a  parabola  whose  axis  is  parallel 
to  the  second  bisectrix  of  the  coordinates:  Y  =  -X.  Only  the  portion 
corresponding  to  the  conditions  X  >  0,  Y  >  0,  that  is  located  in  the 
first  quadrant,  will  be  considered. 
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This  arc  of  the  parabola  is  easily  plotted.   For  X  =  0,  that  is 
X  =   0,  the  two  values 

P   a  (a  +  6 ) 
Xa.  =  ^A  =  ~~ represented  by  point  A 


p   a  +  g 
Yg  =  cog  =  represented  by  point  B 


are  obtained. 

By  referring  to  equation  (13)  or  simply  to  equation  (10)  or  equa- 
tions (8)  and  (9)  it  is  easy  to  identify  that  the  first  of  these  values 
corresponds  to  zero  velocities  and  the  second  to  infinite  velocities. 

On  the  other  hand, 

YB  -  YA  =  ^1(T  -  aD) 

hence,  when  aD  <  T,  that  is  in  the  case  of  divergence  of  the  oscilla- 
tions, point  B  is  above  point  A  (fig.  11);  while,  when  T  <  aD,  that  is 
convergence,  point  B  is  below  point  A  (fig.  13). 

Arranging  equation  (15)  according  to  the  decreasing  powers  of  Y 
under  the  form  AY2  +  BY  +  C  =  0  and  posting  the  condition  B2  -  1+AC  >  0 
to  be  fulfilled  so  that  the  values  of  Y  are  real,  results  in 

Y      2  .  (aD  -  T)2  a  +  e 
Am  -   Ain  -^ 5 ~ — 

16T2    di 

p 
a  value  already  obtained  as  maximum  of  \        in  the  study  of  A,   in 

terms  of  v. 

The  corresponding  value  of  Ym  =  t%2  is  readily  obtained  by  con- 
sidering equation  (15)  developed  under  the  form 

AY2  +  BY  +  C  =  0,   hence  Ym  =  -  — 
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which  gives 


Ym  -  ^gT^VP  +  t)  "  Xm 


hence 


Ym  =  <%  = — (7T2  +  lOaDT  -  a2D2) 

16DIT2 


Returning  to  the  study  of  variations  of  Y,  that  is  of  uj2  when 
v  varies  from  zero  to  infinity,  it  is  apparent  that  the  representative 
point  of  the  variations  of  cd^  and  \2   (figs.  11  and  13)  passes  from 
A  to  B  along  the  parabolic  arc  defined  by  equation  (15)  and  that  B 
is  above  A  when  aD  <  T  and  below  in  the  opposite  case.   In  order  to 
define  the  shape  of  the  representative  curve  of  the  variations  of  the 
frequency  a>  in  terms  of  velocity  v,  it  is  necessary  to  establish 
whether  the  variations  of  a^  are  increasing  or  decreasing  as  in  fig- 
ures 11  and  13  or  can  pass  through  a  minimum  L  as  in  figure  15.   As 
to  the  possibility  of  a  maximum  it  is  apparent  that,  geometrically,  it 
is  excluded  by  reason  of  the  general  position  of  the  parabola. 

For  it  to  have  a  minimum  L  it  is  necessary  and  sufficient  that 
the  slope  of  the  tangents  to  A  and  B  have  the  same  sign.   So,  when 
equation  (15)  is  differentiated,  the  value  of  the  slope  p  at  any  one 
point  is 


dy 

p  =  a     = 
dx 


I  \p       T7 


+  2(X  +  Y) 


I  \D       T/ 


2(X  +  Y) 


Since  the  coordinates  X  and  Y  of  points  A  and  B  are  known,  the 
preceding  formula  gives  the  slopes  p^  and  p^;  forming  the  product 
gives 


PA  x  pb 


(aD  +  3T)(aD  -  5T) 
(aD  -  T)2 


which  for  p.p-r,  >  0  reduces  to  aD  >  5T. 
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This  condition  involves  a  fortiori  aD  >  T;  hence  there  can  be  no 
minimum  except  in  the  case  of  convergence. 

The  discussion  will  he  concluded  with  an  examination  of  the  possi- 
bility of  aperiodic  solutions.   For  this  possibility  to  exist  it  is 
necessary  and  sufficient  that  Y  =  ay?  can  assume  negative  values,  that 
is  that  the  parabolic  arc  AB  intersect  axis  OX.   This  geometrical 
condition  is  expressed  by  the  double  analytical  condition  that  the 
equation  in  X  obtained  by  making  Y  =  0  in  equation  (15)  can  have 
at  least  one  real  and  positive  root. 

The  two  roots  of  this  equation  are  both  real  or  imaginary;  to  be 
real  it  is  necessary  and  sufficient  that  the  discriminant  be  positive, 
which  gives  the  condition 

9T2  -  lOaDT  +  a2D2  >  0 

This  condition  requires  that 

aD  <  T 
or 

aD  >  srr 

afa  +  e  )2 
As,  on  the  other  hand,  the  product  of  the  roots    v —  is  positive, 

I2DT 

the  roots  are  either  both  positive  or  negative.   To  be  positive,  their 

3   a 

D   t 


3   a 
sum  must  be  positive,  which  gives  the  condition  ^-  -  —  <   0  that  is 


aD  >  3T. 

To  satisfy  the  required  double  analytical  condition  it  is  therefore 
sufficient  to  know  the  unique  relation  aD  >  9T;  this  is  the  condition 
necessary  and  sufficient  for  the  existence  of  aperiodic  solutions  for 
suitably  chosen  v.   These  velocities  of  aperiodic  conditions  are, 
obviously,  those  comprised  between  the  velocities  v^,  and  V2,  positive 

solutions  of  the  biquadratic  equation  obtained  by  making  to  =  0  in 
equation  (13).   The  corresponding  convergences  are  obtained  the  same 
way  by  making  a>  =  0  in  equation  (l^).  N 

Now  the  curves  oj(v)   representative  of  the  variations  of  the 
frequency  oo  or,  which  is  the  same,  the  curves  oj2   in  terms  of  veloc- 
ity v  can  be  accurately  plotted.   It  is  accomplished  by  associating 
to  each  curve  aY^v)   the  corresponding  curve  (x^(\^). 
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The  foregoing  discussion  distinguished  four  cases: 

(1)  aD  <  T  divergent  periodic  oscillation  (figs.  11  and  12) 

(2)  T  <  aD  <  5T  convergent  periodic  oscillation,  with  decreasing 
frequency  (figs.  13  and  lk) 

(3)  5T  <  aD  <  9T  convergent  periodic  oscillation  with  maximum 
frequency  (figs.  15  and  l6) 

(k)      9T  <  aD  convergent  aperiodic  oscillation  when  v^_  ^  v -^  Vp 
convergent  periodic  in  the  other  cases.   (figs.  17  and  18) 

B.  Complete  Theory  of  Elementary  Shimmy 

The  equations  (l),  (2),  and  (3)  which  form  the  basis  of  the  sim- 
plified elementary  shimmy  were  obtained  by  assuming  the  moment  torque  M 
arising  from  the  reactions  of  torsion  and  turning,  a  torque  exerted 
directly  on  the  tread  by  the  ground  to  be  negligible.   To  take  this 
moment  M  into  account,  equations  (l)  and  (3)  must  be  suitably  modified, 
while  equation  (2)  remains  the  same. 

First  of  all,  the  plus  sign  is  placed  before  M  when  a  restoring 
torque  due  to  a  positive  displacement  z   is  involved,  just  like  the 
force  F  was  given  the  plus  sign  when  a  restoring  torque  due  to  a 
positive  displacement  z  was  involved. 

Next,  consider  equation  (3);  its  first  term  represents  the  angular 
acceleration  of  the  wheel  about  the  caster  axis.   The  quantity  -F(a  +  e) 
was  the  generating  moment  of  this  acceleration.   In  the  complete  theory, 
the  moment  -M  simply  added  to  this  moment  in  such  a  way  that  this 
equation  becomes 

I  d  z    „,     v 
=  -F(a  +  e)  -  M 

a  dt^ 

As  to  equation  (l),  it  expressed  the  relation  between  the  direction 
dy/dx  of  the  contact  surface  and  force  F.   It  had  been  assumed  by 
writing  this  equation  that  the  drift,  that  is  the  orientation  of  the 
contact  surface  immediately  followed  the  variations  of  force  F;  it 
implicitly  implied  that  the  turn  was  instantaneous,  or  in  other  words, 
implicitly  attributed  an  infinite  value  to  the  turn  coefficient  R.   To 
understand  the  line  of  reasoning  that  leads  to  the  equation  desired, 
we  shall  start  with  the  analysis  of  a  particularly  simple  phenomenon: 
drift. 
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Complete  Study  of  Drift 

Consider  a  pneumatic  tire  at  rest  (fig.  19)-   A  force  F  normal 
to  its  plane  is  applied  and  it  is  made  to  roll  while  maintaining  its 
plane  parallel  to  a  fixed  direction  OP.   Up  to  now  it  had  been  assumed 
that  the  contact  surface  shifts  along  a  straight  line  CM  forming  an 
angle  of  drift  5  =  DF  with  the  plane  of  the  wheel,  that  is  instanta- 
neous drift  was  assumed.   But  experience  indicates  that  the  trajectory 
of  the  contact  surface  or  rather  its  center  is  actually  a  curve  CAB 
tangential  to  OP   in  0  and  asymptotic  to  O'M'  parallel  to  OM. 
This  contact  area  can  turn  only  according  to  the  law  of  turning 


3  =  d 


dy 
dx 


RM  dx 


that  is 


dfy 
dx2 


RM 


But  this  couple  M,  the  generator  of  turning,  is  itself  the  dif- 
ference of  two  effects:   the  accompanying  couple  Fe  and  the  torsional 
couple  of  the  tread  equal  to 

1  dy_ 
S  dx 

in  such  a  way  that 


M   =  Fe   -   I  &L 
S   dx 


hence 


dx2 


=  R 


Fe   -I& 
S   dx 


or 


<1  y  dy 

S  +  R  —  =  RSeF 


dx 


2  dx 
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This  then  is  the  differential  equation  of  the  trajectory  OAB  of  the 
center  of  the  contact  surface  in  drift.   This  equation  is  easily 
integrated 


SeFx  + 


an  equation  that  closely  represents  a  curve  having  the  form  of  that 
plotted  in  figure  19.   Moreover,  this  equation  makes  it  possible  to 
present  an  unusually  interesting  result:   consider  the  drift  dy/dx; 
when  x  increases  indefinitely,   dy/dx  tends  toward  SeF.   Comparison 
with  the  formula 

^=DF 
dx 

it  is  seen  that 


Se  =  D 

This  relation  is  retained  and  applied  repeatedly  in  the  subsequent 
study.   Another  significant  feature  of  this  relation  is  that  it  reduces 
the  number  of  characteristic  parameters  of  the  ordinary  tire  by  one. 
Take  for  example,  T,  D,  e,  R  or  else  T,  D,  S,  R.   The  equation  of 
drift  then  reads 

y  =  DFx  +  —  F\e  S   -  1/ 
R 

In  the  light  of  this  digression  on  drift  we  now  return  to  the  most 
general  motion  problem.   The  path  curvature  of  the  center  of  the  contact 
surface  is  always 

d2y 

— §  =  -RM 

dx^ 


The  generating  moment  M  of  turning  of  the  contact  surface  being  still 

the  difference  of  two  effects,  the  torsion  couple  of  the  tread  equal  to 

a/S,  a  being  the  angle  between  the  axis  of  contact  surface  and  wheel 

dy   z 
plane,  that  is  a  =  —  +  —,   and  the  accompanying  couple  F. 

QX     9, 
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Hence 


M   i  *L  +  -5.  -  F, 
S  dx   Sa 


d2y  1  d2y 

But  can  still  be  written  — -  — -  if  v  is  constant.   Hence  the 


system  of  four  equations 


v2  ,-2 


dx£ 


1  d2y 
_£.  —L  =  _rm 

v2  dt2 


(16) 


TF 


(17) 


I  d£z 
a  dt2 


=  -F(a  +  e) 


M 


(18) 


M  =_1^+  ± 
~       Sv  dt   Sa 


F€ 


(19) 


Entering  the  values  of  F  amd  M   in  equations  (17)  and  (19)  and 
then  in  equations  (l6)  and  (18)  gives 


_^_  d2y    i_   dy   ,     n     _  ,.   =  ( 

Rv2  dt2   Sv  dt   T     Ua   T/ 


Sv  dt   T     a  dt2   VSa   T/ 


(20) 


(21) 


The  elimination  of  y  from  the  two  preceding  equations  then  leaves  the 
differential  equation  in  z  defining  the  motion  of  the  wheel.   Thus 

d3z    .    /ie       1_  a    \d2z 

,  Ta        ci3«ir2 


d^z  I 


Rav2  itk       Sav  dt3 


SRav 


TRv2/  dt 


a  +  e   dz       a  +  e 
+ +  z 

TSv     dt  TSa 


0 


(22) 
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Note  1. -  Setting  R  =  oo  in  the  preceding  equation  and  bearing  in 
mind  the  relation  Se  =  D,  yields  the  equation  (6)  of  the  simplified 
theory. 

Note  2.-  The  elimination  of  y  from  equations  (2Q)  and  (21) 
requires  algebraic  equations  which  can  be  simplified  by  the  following 
considerations;  the  elimination  must  result  in  a  linear  differential 
equation  in  z.   If  n  is  the  degree  of  this  equation  it  will  have 
n  particular  solutions  of  the  form  z  =  Ae   ,  the  coefficients  a 
being  real  or  imaginary.   Conversely,  the  elimination  of  z   leaVes  a 
unique  equation  in  y.   Equation  (20)  indicates  that,  if  y  is  of  the 
form  y  =  A^eat,   z   is  also  of  the  form  z  =  Aeat.   Equation  (21) 
indicates  that,  when  z   is  of  the  form  z  =  Aea  ,  there  also  is  a 
solution  y  of  the  form  A^e0^.   Therefore,  every  solution  of  the 

differential  equation  in  z   is  also  a  solution  of  the  differential 
equation  in  y  and  conversely.   Thus  these  two  equations  are  deduced 
by  simply  changing  y  for  z  and  vice  versa.   The  elimination  of  z 
is  much  easier;  simply  enter  its  value  in  equation  (20)  and  then  sub- 
stitute in  equation  (21).   This  method  makes  it  possible  to  write  equa- 
tion (22)  very  quickly. 

Stability  of  oscillation.-  The  preceding  equation  is  a  linear  differential 
equation  of  the  fourth  order  of  the  form 

&qZ       +  a^z111  +  a2Z   +  a^z1  +  a^z  =  0 

and  all  its  coefficients  are  positive. 

Referring  to  Routh's  "Advanced  Rigid  Dynamics"  it  is  seen  that 
the  conditions  necessary  for  such  an  oscillation  to  be  stable  are 

ala2  ^  a0a3 


aiEpai  >  aga^  +  a-i^ar 

Applied  to  equation  (22)  the  first  of  these  conditions  gives, 
after  reductions, 

ie     T   ^  e 
+  > 


a   SRav2   Rv2 
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which,  with  allowance  for  the  fundamental  relation  Se  =  D,  can  be 
written 

2   aD^  (   j 

IRD 

The  second  of  these  conditions  yields,  after  reductions 

Ta   SRav2   TRv2   Sa2 
or 

(LRv2  -  a) (Sea  -  T)  >  0 

that  is  with  the  fundamental  relation  previously  advanced, 

(IRv2  -  a)(aD  -  T)  >  0  (23a) 

IV.  COMPARISON  OF  THE  PRESENT  THEORY  WITH  THE  EXPERIMENTAL  RESULTS 

Let  us  interrupt  the  examination  of  these  two  stability  conditions, 
and  begin  by  studying  the  second.  To  interpret  it,  assume  that  speed  u 
is  defined  by  the  equality 

2    a 
u  =  — 
IR 

so  the  discussion  can  be  summed  up  as  follows: 

First  case:   aD  >  T 

The  condition  (23a)  is  confirmed  if  v  >  u. 

Second  case:   aD  <  T 

The  condition  (23a)  is  confirmed  if  v  <  u. 

It  is  readily  apparent  that  in  both  cases  equation  (23)  is  confirmed 
when  the  condition  (23a)  is  confirmed.   Therefore  the  conditions  of  the 
stability  of  oscillation  are  as  follows: 

First  case:   aD  >  T  with  v  >  u 

Second  case:   aD  <  T  with  v  <  u 
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This   u   is  termed  the  inversion  speed  of  stability.   The  determina- 
tion of  the  order  of  magnitude  of  this  inversion  speed  u   for  the  usual 
pneumatic  tires  seems  to  be  essential.   But  unfortunately  no  exact 
determination  of  the  value  R  is  available  so  that  this  value  could  be 
introduced  in  the  formula 


u^  =  — 
IE 


Nevertheless,  it  shall  be  shown  that  the  apparently  very  different 
experimental  results  can  be  used  to  define  the  order  of  magnitude  of 
the  inversion  speed  for  the  ordinary  pneumatic  tires. 

Referring  to  the  earlier  study  of  drift  and  to  figure  19,  and 
assuming  that  the  equation  for  the  curve   OAB   is  given,  it  is  readily 

seen  that  ^—   =  e      at  any  point  C   of  this  curve.   Now,  for  the 

DD' 
usual  pneumatic  tires  the  curve  OAB  rapidly  tends  toward  its  asymp- 
totic O'D'M',  and  according  to  the  data  given  by  various  experimenters 
it  seems  permissible  to  assume  that  after  a  travel  of  50  cm 

CD'  <  1 
CD*  *  h 


which,  with  the  centimeter  as  the  unit  of  length,  gives 


1 50 
es 


>  k 


hence 


50  £  >  1.1*0,   that  is   R  >  0.028S 
S 

So,  if  the  values  of  S  for  the  usual  tires  are  known,  the  problem  can 
be  solved.  Unfortunately,  no  such  determinations  are  as  yet  available. 
But  it  was  seen  that  Se  =  D.   Thus  the  preceding  relation  becomes 

R  >  0.028  - 
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On  the  other  hand,  it  has  been  indicated  that  if  Z   denotes  the  length 
of  surface  contact 

1  s       s   l 
6<   6  <2 

So  a  length  of  3  cm  for  e   is  an  altogether  admissible  value  and  prob- 
ably very  close  to  the  real  value  in  a  great  number  of  cases.   The 
preceding  inequality  consequently  becomes 

R  >  °-028  D  or  R  >  0.0093D 

3 

Experiments  made  on  automobile  tires  have  shown  that  under  a  side  force 
of  100  kg  the  drift  5  =  DF  reaches  a  value  approaching  3°,  that  is 

- —  =  0.0521+  radians 
180 

Using  the  C.G.S.    system  of  units   the   coefficient     D     is 

0.052i4-  radians  1A      , 

D  =  =  5-35   lO-10      (C.G.S.    system) 

100.000  X  981  dynes 
hence 

R   >5.10-12      (C.G.S.    system) 

Supposing  now  that  the  caster  length  a  is  5  cm  and  the  inertia  I 
is  10  kg  at  1  m  from  the  axis,  that  is   10°   (C.G.S.  system).   In  that 
case 

U2  =  A  =  5 =  10h 

IR   108  x  5  x  10- 12 


u  =  100  cm/sec  =  1  m/sec  or  3-6  km/hr 
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With  a  10  cm  caster  length  it  would  be 

u  =  1. k   m/sec  =  5.1  km/hr 

With  a  20  cm  caster  length 

u  =  2  m/sec  =7-2  km/hr 

The  essential  part  to  remember  is  that  the  inversion  speed  u  is 
small.   This  is,  as  shall  be  shown,  en  extremely  important  result  and 
it  can  be  regarded  as  incontestable,  although  R  has  been  derived  by- 
indirect  means.   The  form 

2    a 
uc  =  — 

IE 


of  the  equation  used  to  define  u  has,  in  fact,  shown  that  a  given 
variation  of  R   involves  only  a  very  small  variation  of  u,  so  that, 
for  R  one- fourth  as  large,  for  example,  the  real  value  of  u  would 
only  be  doubled.   Furthermore,  the  results  of  this  analysis,  and  in 
particular  the  1  m/sec  obtained  for  the  inversion,  coincide  exactly 
with  the  experimental  curves  published  by  John  Wylie  of  the  Douglas 
Aircraft  Company  (ref.  1).   That  writer  gives  the  exact  coefficient  of 
friction  K  of  the  hydraulic  damper  necessary  to  avoid  shimmy  on  the 
landing  gear  of  the  OA-^A  airplane.   This  curve  is  reproduced  in  fig- 
ure 20.   There  is  no  need  for  a  damper  at  velocities  below  1  m/sec, 
but  there  is  for  all  speeds  above  this  value.   The  case  exhibited  by 
this  airplane  was  the  case  a  <  — .   Therefore  Wylie 's  curve  evidences 

an  inversion  speed  with  stability  of  oscillations  below  this  speed  and 
instability  above  it.   Moreover,  the  inversion  speed  resulting  from 
Wylie's  curve  is  exactly  coincident  with  that  evaluated  in  the  present 
numerical  example. 

The  same  writer  also  gives  a  similar  curve  for  the  nose  wheel  of 
the  Douglas  DC-*;.  The  inversion  speed  of  the  stability  for  this  air- 
plane is  around  1.80  m/sec,  hence  is  still  close  to  that  evaluated  in 
the  present  article. 

The  Douglas  BC-k   in  19^-0  was  the  largest  airplane  which  had  been 
equipped  with  a  nose  wheel  landing  gear.   It  was  a  four-engine  all 
metal,  low  wing  monoplane.   Each  engine  developed  1165  hp. ,  its  weight 
was  29.6  tons  and  its  one  nose  wheel  weighed  300  kg.   It  was  designed 
for  U7  passengers,  had  a  382  km/hr  top  speed  and  119  km/hr  landing  speed. 

A  similar  curve  was  obtained  by  E.  Maier  and  M.  Renz  in  Germany 
during  the  war,  on  the  nose  wheel  of  the  Douglas  DB-7  "Boston."  The 
inversion  speed  seems  to  have  been  near  k   m/sec. 
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As  a  result  of  these  experimental  data  the  existence  of  a  low 
inversion  speed   u  can  be  regarded  as  proved  by  experiment.   Therefore 
referring  to  the  stability  condition  (23)  it  can  be  seen  that,  in  order 
to  eliminate  shimmy  at  all  speeds  above  u,  it  is  sufficient  to  realize 
the  condition  aD  >  T.   At  speeds  below  u,  which  are  low  speeds,  it  is 
evident  that  shimmy  should  not  be  very  annoying  and  can,  perhaps  pass 
completely  unnoticed  in  certain  cases,  because,  as  will  be  shown  later, 
simple  relatively  moderate  friction  in  the  hinge  is  enough  to  damp  it 
out  instantly.   The  consequence  of  this  analysis,  that  is,  the  condition 

aD  >  T 

can  be  practically  regarded  as  the  fundamental  stability  condition  of 
the  wheel. 

Can  this  condition  be  easily  achieved  for  the  permissible  values 
of  caster  lengths  a?  This  is  the  next  problem  to  be  treated. 

If  the  exact  values  of  the  coefficients  T  and  D  for  the  usual 
tires  are  known,  the  answer  is  immediate,  but  unfortunately  there  are 
no  such  data.   So  far  it  has  been  possible  to  evaluate  the  order  of 
magnitude  of  the  coefficients  D,  S,  R,  and  e  on  the  basis  of  earlier 
experiments,  but  still  we  know  nothing  of  coefficient  T.   However, 
the  difficulty  can  be  overcome  by  suitable  interpretation  of  seemingly 
very  dissimilar  experiments. 

We  refer  to  the  report  by  B.  v.  Schlippe  and  R.  Dietrich,  entitled 
The  Mechanics  of  Pneumatic  Tires  (ref.  2)  which  is  a  detailed  study  of 
the  behavior  of  the  pneumatic  tire.   Unfortunately,  the  authors  failed 
to  give  the  simple  properties  that  control  the  mechanics  of  the  tire, 
as  they  did  not  take  up  the  study  of  the  most  complex  phenomenon, 
shimmy.   Nevertheless,  their  report  supplies  some  significant  informa- 
tion because  it  contains  a  certain  amount  of  numerical  data.   They 
made  experiments  and  measurements  on  a  260  X  85  tail  wheel  fitted  with 
a  Continental  balloon  tire  with  small  longitudinal  grooves  carved 
around  the  circumference.   The  tire  pressure  was  2.5  atm  and  the  load 
180  kg.   The  wheel  was  connected  to  different  dynamometers  and  kept 
stationary,  and  rolling  was  accomplished  by  a  rotating  drum  90  cm  in 
diameter,  covered  with  emery  paper  to  assure  good  adhesion. 

The  first  experiment  by  these  writers  to  which  attention  is  called, 
is  the  following:   when  a  tire  is  made  to  roll  by  keeping  the  plane  of 
the  wheel  fixed  and  imposing  on  it  a  rectilinear  trajectory  making  a 
certain  angle  0     with  the  wheel  plane  (fig.  22),  the  center  M  of  the 
contact  surface  describes  a  curvilinear  path  MM'  having  as  asymptote 
the  straight  line  (D)  located  a  distance  m  from  the  trajectory  00'  of 
the  wheel  center. 
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The  authors  established,  by  experiment,  the  proportionality 

m  =  K0 

K  being  equal  to  10  if  6     is  expressed  in  radians  and  m  in  cm. 

When  considering  the  side  force  F  necessary  to  keep  the  wheel 
on  the  trajectory  00',  the  angle  B      is,  obviously,  the  angle  of  drift, 
hence  9   =  DF.   On  the  other  hand,   m  is  evidently  the  displacement  of 
the  lateral  elasticity  and  consequently  m  =  TF,  hence 

m  _  t 

9  ~  D 

hence 

m  =  -  e 

D 

Their  coefficient  K  is  none  other  than  the  quotient  T/D,  and 
the  experiment  in  figure  22  gives  this  quotient  directly.   The  quotient 
for  this  wheel  is,  thus,  equal  to  10  cm,  and  the  caster  length  necessary 
to  eliminate  shimmy  is  a  ^  10  cm. 

This  particular  example  indicates  that  the. caster  lengths  to  which 
the  present  theory  leads  are  of  a  reasonably  approximate  magnitude,  at 
least  in  certain  cases.   Hence  nose-wheel  shimmy  can  be  probably 
eliminated  by  a  simple  modification  of  the  position  of  the  pivoting 
axis  of  the  wheel,  thus  making  the  use  of  shimmy  dampers  unnecessary. 

The  experiments  by  Schlippe  and  Dietrich  are  of  further  interest 
for  another  reason:   they  enable  the  characteristic  coefficients  T,  D, 
S,  R,  and  e  of  the  tire  to  be  found  indirectly,  which,  in  the  absence 
of  more  adequate  measurement,  helps  in  defining  the  approximate  magni- 
tude of  these  fundamental  factors.   Their  direct  measurement  of  the 
lateral  elasticity  of  the  wheel  indicated  that  a  side  force  F  of 
2  X  32.5  kg,  or  65  kg,  was  necessary  to  produce  a  1  cm  lateral  displace- 
ment A  of  the  tread  with  respect  to  the  wheel.   Now  if  A  =  TF,  we 
set  in  C.G.S.  units 


65OOO  X  981 
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or 

T  =  157  X  10-10   (C.G.S.  system) 

Since  -=  10,  we  get  D  =  15-7  X  10"10   (C.G.S.  system) 
D 

These  authors  also  indicated  that,  to  maintain  an  angle   0=1° 
or  0.0175  radians  in  the  experiment  described  above  (fig.  22),   a  side 
force  of  11. k   kg  and  a  torque  having  a  moment  C  of  hty.J   cm  X  kg  must 
be  exerted.   So,  since  9  =   DF ,   we  should  have 

0.0175  =  D  X  11 A  X  981.OOO 
hence 

D  =  15.6  X  10" 10   (C.G.S.  system) 

The  agreement  between  this  figure  and  that  obtained  independently  some 
lines  back  is  perfect. 

The  accompanying  torque  measured  by  these  authors  was  C  =  ^9> 7  cm.  kg. 
It  is  known  that 

C  =  eF 

K9.  7 
hence   e  =  — - —  or   €  =  k.  35  cm. 

11.1+ 

This  figure  is  compared  with  the  length  2  of  the  tread  or  ground  con- 
tact surface.   The  authors  indicate  that  this  length  was  about  9  cm. 
Hence 

1   -    -2-  =2.1 


e 


1+.35 


This  result  is  *in  close  agreement  with  the  formula 

2<  |<  6 

given  at  the  beginning  of  the  present  article  while  analyzing  the 
accompanying  torque. 


3^ 
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The  torsional  elasticity  was  measured  by  the  difference  in  lateral 
displacement  Z  of  the  lead  point  and  Z'   of  the  tail  point  of  the 
tread  contact  in  such  a  way  that  if  I      is  the  length  of  the  tread,  the 
torsion  angle  a  is  such  that 


sin  a  = 


therefore,  for  the  small  angles 


a  = 


For  these  conditions  the  authors  set  up  the  formula 

M  =  d(Z  -  Z«) 

where  M   is  the  moment  necessary  to  produce  the  torsion  and  d   is  a 
constant  equal  to  317  kg,  where  Z  and  Z'  are  expressed  in  cm  and 
M  in  cm.  kg,  that  is   d  =  317  X  98IOOO  (C.G.S.).   Now  Z  -  Z'  =  la, 
hence  the  formula  giving  M  becomes 

M  =  dZa 


However,  since  a  =  SM  and   1=9  cm, 


S  =  ^ 


that  is,   S  =  3-57  X  10"10   (C.G.S.  system) 


dl       9  x  317  X  98IOOO 
On  the  other  hand,  D  =  Se   and  6  =  k. 35  cm.   Therefore, 
D  =  ^-35  X  3.57  X  10" 10  that  is  D  =  15 . 5  X  10"10   (C.G.S.  system) 


Two  independent  methods  have  already  given 


D  =  15.7  X  10" 10  and  D  =  15-6  X  10"10 


The  agreement  between  these  three  figures  is  really  remarkable,  it  even 
seems  to  exceed  the  probable  accuracy  of  the  measurement.   At  any  rate, 
these  results  give  a  completely  satisfactory  experimental  check  of  the 
present  theory. 

The  one  characteristic  left  to  define  is  the  turn  coefficient  R. 
Its  determination  will  be  based  upon  the  following  experiment:   by 
causing  the  tire  to  roll  against  or  on  a  circular  wheel  of  1  m  radius, 
the  authors  claimed  that  the  motion  produced  an  axial  force  F  of 
21.2  kg,  tending  to  keep  the  tire  away  from  center  of  the  periphery. 
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If,  as  seems  likely,  according  to  the  conditions  of  the  experiment,  it 
is  assumed  that  the  generating  moment  of  M  of  the  rotation  arises 
solely  from  the  accompanying  torque 

Fe  =  21.2  x  981000  x  k.35  =   90.5  X  106   (C.G.S'.  system) 
the  coefficient  R  can  be  deduced  by  the  formula 


d2y 
dx2 


KM 


If  p   signifies  the  curvature  radius  of  the  track  the  formula 
must  be  used,  hence 


d2y 
dx2 


1 
P 


R  =  -T7  that  is 
pM 


R  = 


100  x  90.5  X  10c 


or  R  =  1.1  X  10"10  (C.G.S.  system) 


It  will  be  noted  that  this  value  of  R  is  in  good  agreement  with 

the  inequality  R  >  5  X  10" -^   indicated  previously.   Still,  this  deter- 
mination does  not  offer  the  same  degree  of  certainty  as  the  preceding 
ones:   in  fact  since  it  was  necessary  to  assume  that  the  wheel  was 
exactly  perpendicular  to  the  radius  of  the  track,  any  error  in  this 
special  circumstance  will  modify  the  centrifugal  force  F  by  the 
superposition  of  a  drift  effect.   To  eliminate  this  potential  source 
of  error  completely  it  is  necessary  to  assure  the  equality  of  the 
centrifugal  force  F  by  changing  the  direction  of  rotation.   In  the 
case  of  a  minor  discrepancy,  the  average  should  be  taken. 

U.S.  Reports  on  Shimmy 

The  problem  of  nose-wheel  shimmy  has  already  formed  the  object  of 
numerous  theoretical  and  experimental  studies,  expecially  in  the 
United  States  where  the  so-called  tricycle  landing  gear  was  born.   It 
was  studied  in  great  detail  by  the  Douglas  Aircraft  Company,  then  a 
little  later,  by  the  Lockheed  Aircraft  Corporation.   On  the  occasion 
of  these  studies  a  mathematical  theory  of  shimmy  was  suggested  by  Wylie 
(ref.  l)  and  by  Arthur  Kantrowitz  of  the  Langley  Memorial  Aeronautical 
Laboratory  of  the  National  Advisory  Committee  for  Aeronautics  (NACA 
Rep.  686). 


Because  it  is  the  only  theory  giving  a  numerical  account  of  the 
phenomena  of  shimmy,  it  is  deemed  practical  to  reproduce  a  literal 
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exposition  of  Kantrowitz's  article.  This  theory  is  proposed  as  being 
"based  on  the  discovery  of  a  new  phenomenon  called  kinematic  shimmy." 
Here  is  the  exposition  of  Kantrowitz. 


1.  Kinematic  Shimmy 

Some  preliminary  experimental  results  on  shimmy  were  obtained  by 
the  N.A.C.A.  with  the  aid  of  the  belt-machine  apparatus  shown  in  fig- 
ure 23-   This  machine  consists  of  a  continuous  fabric  belt  mounted  on 
two  rotating  drums  and  driven  by  a  variable-speed  electric  motor. 
Provision  is  made  for  rolling  a  caster ing  wheel  up  to  about  6  inches 
in  diameter  on  the  belt  in  such  a  way  that  it  is  free  to  move  vertically 
but  not  horizontally. 

On  this  belt  machine,  the  following  phenomenon  (see  fig.  23)  was 
discovered  while  pushing  the  belt  very  slowly  by  hand.   With  the  wheel 
set  at  an  angle  with  the  belt  as  in  (a)  and  the  belt  pushed  slowly,  the 
bottom  of  the  tire  would  deflect  laterally  as  is  shown  in  (b).   When 
the  belt  was  pushed  farther,  the  wheel  straightened  out  gradually  as 
is  shown  in  (c).   The  bottom  of  the  tire  would  then  still  be  deflected, 
however,  and  the  wheel  would  continue  to  turn  as  in  (d).   The  wheel 
would  thus  finally  overshoot,  as  shown  in  (e)  and  (f).   The  process 
would  then  be  repeated  in  the  opposite  direction. 

Figure  2k   is  a  photostatic  record  of  the  track  left  by  the  bottom 
of  the  tire  on  a  piece  of  smoked  metal.   Two  things  will  be  noticed: 
First,  that  the  bottom  of  the  tire  did  not  skid;  and,  second,  that  the 
places  where  the  wheel  angle  is  zero  (indicated  by  zeros  on  the  track) 
correspond  roughly  to  the  places  where  the  lateral  deflection  of  the 
tire  is  a  maximum.   Thus  the  wheel  angle  lags  the  tire  deflection  by 
one -quarter  cycle. 

It  was  noticed  that  the  oscillation  could  be  interrupted  at  any 
point  in  the  cycle  by  interrupting  the  motion  of  the  belt  without 
appreciably  altering  the  phenomenon.   From  this  observation  it  was 
deduced  that  dynamic  forces  play  no  appreciable  part  in  this  oscillation. 

The  distance  along  the  belt  required  for  one  cycle  was  also  found 
not  to  vary  much  with  caster  angle  or  caster  length.   (See  fig.  25.) 
Caster  length  was  therefore  considered  not  to  be  of  fundamental  importance 
in  this  type  of  oscillation. 

It  should  be  pointed  out  that,  in  order  to  observe  the  kinematic 
shimmy,  lateral  restraint  of  the  spindle  is  necessary  to  prevent  the 
spindle  from  moving  laterally  when  the  bottom  of  the  tire  is  deflected 
and  thus  neutralizing  the  tire  deflection.   This  restraint  is  supplied 
by  the  dynamic  reaction  of  the  airplane  when  the  airplane  is  moving 
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forward  rapidly  but  is  not  ordinarily  present  when  the  airplane  is 
moving  forward  slowly.   It  has  been  observed,  however,  on  airplanes 
towed  slowly  with  two  towropes  so  arranged  as  to  provide  some  lateral 
restraint. 

Figure  23(a)  shows  that,  when  the  center  line  of  the  wheel  is  at 
an  angle  8      (see  fig.  25)  with  the  direction  of  motion,  the  bottom  of 
the  tire  deflects.   This  situation  is  represented  schematically  in  fig- 
ure 26.   It  is  seen  that  a  typical  point  on  tne  peripheral  center  line 
must  have  a  component  of  motion  perpendicular  to  the  wheel  center  line 
if  the  tire  is  not  to  skid.   Thus 

dX  =  -sin  9   ds 

(The  minus  sign  follows  from  the  conventions  used  as  shown  in  fig.  25.) 
Since  only  small  oscillations  are  to  be  considered,  the  approximation 

*L=   -6  (1) 

ds 


may  be  substituted. 

The  effect  of  tire  deflection  on  6     will  now  be  considered.   For 
the  purposes  of  rough  calculation,  it  will  be  assumed  that,  as  illus- 
trated in  figure  27,  the  projection  of  the  peripheral  center  line  on 
the  ground  is  a  circular  arc  intersecting  the  wheel  central  plane  at 
the  extremities  of  the  projection  of  the  tire  diameter.   (See  fig.  23(d).) 
Thus,  in  figure  27,   r  is  the  tire  radius.   (it  will  be  assumed  for 
the  time  being  that  the  caster  length  and  the  caster  angle  are  zero. ) 
Now  if  the  tire  is  deflected  in  the  form  of  a  circular  arc,  then  the 
condition  that  the  torque  about  the  spindle  axis  be  zero  is  that  the 
strain  be  symmetrical  about  the  projection  of  the  wheel  axle  on  the 
ground.   Clearly,  if  the  wheel  is  displaced,  it  will  be  turned  about 
the  spindle  axis  by  the  asymmetrical  elastic  forces  until,  if  it  is 
allowed  time  to  reach  equilibrium,  the  symmetrical  strain  condition  is 
reached.   Thus,  if  the  tire  is  deflected  an  amount  X  as  in  figure  27(a) 
and  if  the  wheel  rolls  forward  a  distance  ds   to  the  condition  shown 
in  figure  27(b),  in  order  for  the  strain  to  remain  symmetrical  the  wheel 
must  turn  about  the  spindle  axis  an  amount  d0.   From  figure  27, 
Rd9  =  ds.   The  value  of  R  may  be  readily  obtained  from  geometry  in 
terms  of  r  and  X.   Thus   R2  =  r2  +  (R  -  x)2 f    from  which,  if  X2  «  r2 f 
it  is  seen  that  R  =  r2/2X.   Then  substituting  for  R, 

^=^X  (2) 

ds    r2  K    ' 
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If  the  caster  length  is  finite,  the  strain  will  not  be  symmetrical 
about  the  axle,  as  was  assumed  here,  but  will  be  symmetrical  with  respect 
to  some  line  parallel  to  the  axle  but  a  certain  fixed  distance  ahead. 
Hence,  the  essential  elements  of  the  geometry  are  unchanged  and  all  the 
reasoning  that  led  to  equation  (2)  is  still  valid  for  this  case. 

Since  the  phenomena  represented  by  equations  (l)  and  (2)  occur 
simultaneously,  they  must  be  combined  to  get  the  total  effect.   Thus 

^--^B  (3) 

ds2     r2 


This  differential  equation  corresponds  to  a  free  simple  harmonic 
oscillation  occurring  every  time  the  wheel  moves  a  distance 

S  =  -%==   «£. 
fr2 


Kb- 

Measurements  of  the  space  interval  S  of  kinematic  shimmy  have 
been  made  for  three  similar  tires  of  the  type  illustrated  in  figure  23. 
These  tires  all  had  radii  of  approximately  2  inches  so  that  the  theo- 
retical interval  was  about  0.7^  foot.   Their  experimental  intervals 
were  0.65  foot,  0.7^  foot,  and  0.79  foot.   This  agreement  is  closer 
than  might  have  been  anticipated  in  view  of  the  roughness  of  the  assump- 
tion.  It  will  be  seen  from  equation  (l)  that  X      is  one-fourth  cycle 
out  of  phase  with  9,  which  is  in  agreement  with  the  information 
obtained  from  figure  2k. 

In  order  to  take  account  of  tires  for  which  the  assumption  made 
concerning  the  projection  of  the  peripheral  center  line  is  not  quantita- 
tively valid,  an  empirical  constant  K  will  be  used  in  place  of  2/r2 
in  equations  (2)  and  (3)>  thus  obtaining 


and 


ds 


d£e 

ds2 


=  -K0  (5) 


The  constant  K  can  be  measured  by  observing  the  space  interval  of 
kinematic  shimmy.   Where  experimental  values  of  K  are  available,  they 
will  be  used  rather  than  the  rough  theoretical  value  2/r2. 
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2.  Dynamic  Shimmy- 
In  the  foregoing  derivation  for  the  oscillation  called  kinematic 
shimmy,  it  was  assumed  that  the  strain  of  the  tire  was  always  symmetrical, 
that  is,  the  wheel  was  moving  so  slowly  that  any  torque  arising  from 
dynamic  effects  involved  in  the  oscillation  would  be  negligible.   For 
this  case,  from  equation  (k) 


(x   -  i  &) 

\         K  ds/ 


If,  now,  the  wheel  is  assumed  to  be  moving  at  a  velocity  such  that 
the  effect  of  the  moment  of  inertia  about  the  spindle  axis  is  signifi- 
cant, then  the  strain  can  no  longer  be  symmetrical  and,  for  small 
asymmetries,  the  torque  exerted  by  the  tire  on  the  spindle  will  be 
proportional  to  the  amount  of  the  asymmetry.   Thus,  the  value  in  paren- 
theses will  no  longer  be  zero  but  it  can  be  assumed  that  it  will  be 
proportional  to  the  dynamic  torque;  hence 

dt2    \   K  ds 

where  C-^  is  an  appropriate  constant  of  proportionality  and  includes 
the  moment  of  inertia  about  the  spindle  axis   I.   If  the  forward  veloc- 
ity 'V  of  the  wheel  is  constant 

d29  _  v2  df£ 
dt2      ds2 

and 


V2  d!i  =  C   (x   _   1  &)  (6) 

ds2    \    K  d3/ 

The  constant  C^  may  be  determined  by  deflecting  the  bottom  of 
the  tire  a  known  amount,  moving  the  wheel  forward,  and  balancing  the 
torque  M  exerted  by  the  tire  on  the  spindle  so  that  9     stays  constant. 
The  method  of  deflecting  the  bottom  of  the  tire  a  known  amount  will  be 
described  later.   In  this  case  (6)  becomes 


1+0 
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M 


(6a) 


from  which  C]_  may  be  found.   It  has  been  observed  that  C-j_  increases 
with  increasing  caster  angle.   Thus,  for  a  tire  like  the  one  in  fig- 
ure 23,   C]_  was  71*000  radians  per  second^  per  foot  for  a  caster 
length  L  of  0.17  inch  (caster  angle,  5°;  no  fork  offset)  and  was 
10^,000  radians  per  second  per  foot  for  caster  length  0.68  inch 
(caster  angle,  20°). 


In  the  study  of  kinematic  shimmy,  it  was  also  seen  that  the  only 
change  in  X  was  due  to  the  fact  that  a  component  of  the  forward  motion 
was  perpendicular  to  the  central  plane  of  the  wheel.   This  circumstance 
is  expressed  by  equation  (l).   It  is,  however,  found  that,  if  the  spindle 
is  clamped  at  9=0°     and  the  bottom  of  the  tire  is  deflected,  the 
deflection  will  gradually  neutralize  itself;  that  is,  the  bottom  of  the 


tire  will  roll  under  the  wheel. 


that  exists  in  this  case  contributes  to 
d0/ds   zero  is  illustrated  in  figure  28. 
to  be  proportional  to  the  cause,  there  is 


Thus  the  asymmetrical 
dX/ds .   The 


£") 


strain 


case  of  9     and 
If  the  effect  is  again  supposed 
obtained 


dX 
ds 


^  =  -C2X 


(7) 


In  this  equation,  the  constant  C2   is  a  geometrical  constant  of 

the  tire  that  can  be  obtained  from  static  measurements.   The  order  of 

magnitude  of  Co  can  be  obtained  by  assuming  that  the  periphery  of  the 

tire  intersects  the  extremity  of  the  extended  central  plane  of  the  wheel. 

In  that  case   Co  =  — . 

r 

If  9      is  not  zero,  there  will  be  a  component  of  the  forward  motion 
contributing  to  dX/ds.   As  in  equation  (l),  this  component  will  be  -9. 
Adding  this  component  to  the  part  of  dX/ds   due  to  asymmetry,  then 
(d0/ds  still  assumed  zero) 


dX 
ds 


—  =  -0  _  CoX 


(8) 


This  equation  expresses  that,  for  —  =  0,  the  contribution  of  the 

ds 

asymmetrical  strain  to   dX/ds   was   -CoX.   Also  for  the  symmetrical 

strain,  in  which  case   X  -  — L-^-   =  0   (kinematic  shimmy),  there  was, 

K 

of  course,  no  contribution  due  to  asymmetry.   Assume  now  a  linear  inter- 
polation between  these  two  limiting  cases.   Thus,  finally, 
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dX. 
ds 


K  ds/ 


(9) 


A  method  of  determining  C2   is  provided  by  equation  (7)  which, 
when  integrated,  gives 

log^  —  =  -  C2(s  -  sq) 

If,  with  the  spindle  clamped  at  9  =   0°,  the  tire  is  deflected  a 
known  amount  Xq     and  rolled  ahead  a  known  distance  and  the  new  X 
measured,  C2  may  be  computed.   It  was  seen  earlier  that  C2  was  of 
the  order  of  magnitude  of  l/r;  that  is,  it  would  be  of  the  order  of  6 
for  a  2-inch-radius  tire.   The  constant  C2  was  determined  for  two 
model  tires  under  different  loads  and  found  to  be  6.2  and  3.^-   Con- 
siderable variations  of  this  constant  with  tire  pressure  and  load  have 
been  found. 

A  method  of  obtaining  the  constant  known  X     necessary  for  the 
measurement  of  C]_  is  provided  through  equation  (8).   Here  it  is  seen 
that,  if  the  wheel  is  pushed  along  at  a  constant  angle  0,      X     will 

increase  (negatively)  until  8  =   -C2X.,  in  which  case  ■= —  =  0  and 
equilibrium  is  reached. 

When  the  wheel  is  moving  ahead  at  a  finite  constant  velocity,  the 
phenomena  represented  by  equations  (6)  and  (9)  occur  simultaneously. 
Therefore,  to  get  the  total  effect,  combine  the  two  equations,  thus 
obtaining 

„9  ,3„   /-.   r~  _\  ,P~ 

=  0  (10) 


The  solution  for  the  natural  modes  of  motion  represented  by  equa- 
tion (10)  is 


a-rs     ctps     aos 
8   =  Ae  L     +   Be  d     +   Ce  3  (ll) 

where   the     a's     are   the   three   solutions   of  the  so-called  auxiliary 
equation 

y£  a3  +  (k  +  £2  v2V2  + 1  =  0  (12) 

Cl  U      ci     / 
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One  of  these  ex's   is  real  and  negative  and  corresponds  to  a  non- 
oscillatory  convergence.   The  other  roots  are  conjugate  complex  numbers 
and  correspond  to  the  shimmy  under  consideration.   The  roots  will  be  of 
the  form  a  ±  coi.   If  the  divergence  a  is  positive,  the  oscillation 
will  steadily  increase  in  amplitude  (while  its  amplitude  is  not  large 
enough  for  skidding  to  occur);  and,  if  a   is  negative,  the  oscillation 
will  steadily  decrease  in  amplitude  and  eventually  disappear.   The 
meaning  of  the  quantity  "divergence"  may  be  illustrated  by  saying  that 
it  is  approximately  equal  to  the  natural  logarithm  of  the  ratio  of 
successive  maximum  amplitudes  to  the  distance  between  them.   The 
quantity  ix>     is  equal  to  2jt  times  the  number  of  oscillations  per  foot. 
The  frequency  therefore  is  cjdV/2jt.   The  phase  angle  is  obtained  by  sub- 
stituting for  9      in  equation  (6)  the  value  obtained  from  the  foregoing 
procedure  and  solving  for  X. 

The  divergence,  the  frequency,  and  the  phase  relations  thus  derived 
for  typical  model  tire  constants  are  plotted  in  figure  29.   For  small 
velocities  (0  to  6  ft  per  sec)  the  frequency  corresponds  to  kinematic 
shimmy;  it  is  proportional  to  velocity.   The  divergence  increases 
rapidly,  however,  because  the  spindle  angle  lags  on  account  of  the 
moment  of  inertia  about  the  spindle  axis,  thus  allowing  more  lateral 
deflection  than  would  occur  in  a  kinematic  shimmy.   On  the  next  half 
cycle,  a  larger  spindle  angle  is  reached  and  the  process  repeats.   As 
the  velocity  is  further  increased,  the  lag,  and  hence  the  asymmetry  of 
the  strain,  further  increase  until  the  strain  becomes  almost  entirely 

asymmetrical.   For  this  condition,  $>§■  —  «  X.      Then  the  restoring  torque 

on  the  spindle  is  approximately  proportional  to  X.   (See  equation  (6).) 
The  tire  deflection  X      (measured  negatively)  will,  however,  still  lag 
somewhat  behind  9     because,  after  the  wheel  is  turned  through  a  given 
angle,  a  certain  forward  distance  is  required  for  equilibrium  tire 
deflection  to  be  reached.   Thus,  the  restoring  force  will  again  lag 
the  displacement.   As  the  velocity  increases  in  the  high-velocity  range, 
the  frequency  stays  nearly  constant  (see  fig.  29)  and  the  distance 
corresponding  to  a  single  oscillation  increases.   Hence  this  constant 
lag  becomes  a  smaller  part  of  the  cycle  and  the  divergence_  decreases 
at  high  velocities. 

It  will  be  appreciated  that  the  foregoing  theory  considers  only 
the  fundamental  phenomena  taking  place  in  shimmy.   Other  phenomena 
occurring  simultaneously  have  been  neglected.   Some  of  the  more  important 
of  the  neglected  phenomena  are: 

1.  Miscellaneous  strains  (other  than  lateral  tire  deflection) 
occurring  in  the  tire.   A  rubber  tire  being  an  elastic  body  will  distort 
in  many  complicated  ways  while  shimmying.   In  particular,  there  will 
be  a  twist  in  the  tire  due  to  the  transmission  of  torque  from  the 
ground  to  the  wheel. 
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2.  Two  effects  will  cause  the  stiffness  constants  of  the  tire  to 
change  with  speed.   First,  centrifugal  force  on  the  rubber  will  make 

the  tire  effectively  stiffer  at  high  speeds.   Second,  much  of  the  energy- 
used  to  deflect  the  tire  will  go  into  compressing  the  air.   The  com- 
pressibility of  the  air  will  change  with  the  speed  of  compression  owing 
to  the  different  amounts  of  heat  being  transferred  from  it. 

3.  There  will  be  a  gyrostatic  torque  about  the  spindle  axis  caused 
by  the  interaction  of  the  rotation  of  the  wheel  on  the  axle  and  the 
effective  rotation  of  part  of  the  tire  about  a  longitudinal  axis  on 
account  of  the  lateral  tire  deflection.   This  torque  will  later  be 
shown  to  have  a  noticeable  effect  on  the  results. 

The  inclusion  of  items  1  and  2  in  the  theory  would  obviously  be 
very  difficult.   It  is  therefore  necessary  to  resort  to  experiment  to 
determine  whether  the  present  theory  gives  an  adequate  description  of 
the  phenomena.   If  so,  the  omission  of  these  and  any  other  items  will 
be  justified. 

An  experimental  check  on  the  theory  was  obtained  by  measuring  the 
divergence  and  the  frequency  of  the  shimmy  on  the  belt  machine  at  two 
caster  angles  and  at  a  series  of  velocities.   These  measurements  were 
made  by  placing  a  lighted  flashlight  bulb  on  a  6- inch  sting  ahead  of  a 
model  castering  wheel  with  a  ball-bearing  spindle  and  then  taking  high- 
speed moving  pictures  of  the  flashlight  bulb  with  the  wheel  free.   The 
photographs  were  made  with  time  recordings  on  the  film,  and  the  belt 
carried  an  object  that  interrupted  light  from  a  fixed  flashlight  bulb 
and  thus  recorded  the  belt  speed  on  the  film. 

The  divergence  and  the  frequency  of  the  shimmy  were  obtained  by 
measuring  the  displacements  and  the  times  corresponding  to  successive 
maximum  amplitudes  (while  the  amplitude  was  still  small  enough  to  make 
all  the  assumptions  valid).   The  results  are  plotted  in  figure  30. 

In  order  to  compare  these  results  with  the  theory,  the  constants  C]_, 
C2,  and  K  were  determined  on  the  same  tire  at  the  two  caster  angles  by 
the  previously  described  methods.   It  was  found  that  C2  =  6.  2  feet"^-; 
that  K  =  62.5  feet"2;  and  that,  for  5°  caster  angle,   C"i  =  71,100  feet"2- 
second-2  and,  for  20°  caster  angle,  Ci  =  10^,000  feet-1-second~2.   The 

roots  of  equation  (12)  were  then  found  and  the  divergence  and  the  fre- 
quency of  the  shimmy  were  computed  for  a  series  of  velocities  and  at 
caster  angles  of  5°  and  20°.   These  results  are  plotted  in  figure  31  and, 
for  purposes  of  comparison,  the  experimental  curves  are  also  reproduced. 

The  agreement  between  theory  and  experiment  is  considered  satis- 
factory as  regards  qualitative  results.  It  will  be  noticed,  however, 
that  the  theoretical  values  of  the  divergence  are  decidedly  too  large 
at  high  velocities,  say  25  feet  per  second. 
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Kantrowitz 's  theory  and  his  experimental  results  will  now  be 
compared  with  the  new  theory  proposed  in  the  present  article: 

First  of  all,  the  equation  of  motion  (22)  of  the  wheel  has  for  the 
general  solution  the  function 

z  =  C^e8!*  +  C2ea&  +  C3eS3t  +  Cke3^ 


where  C-^,  C2,  ^3,  ^h      indicate  four  arbitrary  parameters  and  S]_,  S2, 
so,   sjj.  the  four  roots  of  the  characteristic  equation 

I    k    I   o   /ie     1      a  \  o   a  +  e     a  +  e 

34  +  SJ  +  I  +  +  \Sd    + 


Rav2      Sav      VTa   SRav2   TRv2/      TSv       TSa 


The  coefficients  of  this  equation  being  always  positive,  the  four 
roots  are  real,  negative  or  conjugate  negative  in  groups  of  two.   The 
negative  roots,  if  existing,,  correspond  to  convergent  terms  and  there- 
fore cannot  be  generators  of  shimmy;  as  a  result  shimmy  must  arise 
from  the  imaginary  roots.   Assuming  now  that  S]_  and  39  are  two  con- 
jugate complex  numbers  of  the  form  X  ±  toi.   The  corresponding  expo- 
nentials are  then  imaginary  and  it  is  advisable  to  modify  the  expres- 
sion.  Transformed  in  trigonometrical  terms  by  means  of  the  Euler 
formulas,  these  two  exponentials  are  combined  to  produce  a  term  of  the 

form  Ae^tsin(cut  +  <p),  that  is  an  exponential  sinusoidal  oscillation, 
that  is  stable  when  X   is  negative  and  unstable  when  X   is  positive. 

The  theoretical  study  of  shimmy  can  now  be  successfully  completed 
with  the  calculation  of  the  variations  of  X  and  a)  in  terms  of  the 
coefficients  of  equation  (22).   The  same  direct  method  used  in  the 
simplified  theory  for  defining  the  curves  of  X  and  <x>  in  the  terms 
of  velocity  v  could  be  applied,  but  the  task  would  be  drawn  out  and 
difficult.   In  fact  it  will  be  shown  that  it  is  not  absolutely  necessary 
to  interpret  the  phenomenon  of  shimmy  quantitatively  and  qualitatively 
in  its  minute  details.   Use  will  be  made  of  Kantrowitz's  theory  and 
experimental  data. 

We  shall  first  examine  what  the  kinematic  shimmy  can  be.   The 
record  of  the  track  of  the  tire  on  the  ground,  the  place  of  the  center  I 
of  the  control  surface,  of  elongation  y,  was  expressed  by  an  equation 
where  y  had  the  same  form  as  equation  (22).   Replacing  z  by  y  in 
equation  (22),  and  taking  the  path  distance  x  instead  of  time  t  as 
the  independent  variable,  yields 
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dy  __  dy  dx  _  ^  dy 

*2y_ 

v2  iz 

dt    dx  dt     dx 

dt2 

dx2 

d3y  _  v3  d3y 

dt3      dx3 

dV 
dt^  = 

dx^ 

and  equation  (22)  becomes 


Iv2  dV  +  Iv£  d3y  +  /lev2  +  _I_  +  a\d2y  +  a  +  e  dy  +  a  +  e     0   (2m 
Ra  dxl+   Sa  dx3   \^Ta    SRa   TR/dx2    TS   dx    TSa 


When  the  velocity  v  is  very  small  and  approaches  zero,  all  the  terms 
containing  v  disappear  and  the  preceding  equation  becomes  equivalent 
to 


,2, 


(T  +  Sa2)d-i  +  Ra(a  +  e)—  +  R(a  +  e)y  =  0  (25) 

dx2  dx 

Such  an  equation  represents  a  damped  oscillating  motion.   The  tire 
track  will  therefore  be  a  damped  sinusoid  of  equation 

y  =  Ae-)"xsin(aDX  +  cp) 

with 

P   k(T  +   Sa2)R(a  +  e)  -  R2a2(a  +  e)2  ,    „ 

of   = (26) 

MT  +  Sa2)2 


and 


Ra(a  +  e)  . 

X   =  — : '—  (27) 

2(T  +  Sa2) 
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The  wave  length  on  the  ground,  that  is  the  period  W,  is  given  by 

w  =  2£ 


CD 


At  zero  caster  length  a 


CD2  =  51       W 

T 


So,  if  a  =  0,  the  damping  is  zero,  that  is  the  tire  track  on  the 

f¥ 

ground   is   a  sinusoid  of  wave    length      W  =   2jt  l/ — .      Thus    the    first   point 
made   by  Kantrowitz   and  the   phenomenon  of  kinematic   shimmy   is    explained. 

An  attempt  will  now  be  made  to  compare  the  computed  wave  length 

with  that  obtained  experimentally  by  the  American  writer.   To  this  end 

the  characteristics  obtained  for  the  tire  studied  by  Schlippe  and 
Dietrich,  that  is 

T  =   157  X   10-10  C.G.S.    system 

S  =   3-57  X   lO"10  C.G.S.    system 

R  =   1.1  X   10" 10  C.G.S.    system 
e   =  h. 35  cm 


are   used. 

The   result    is 


W  =   2*1/ 121 =    36   cm 

lr  l.l  x  1+.35 

This  wave  length  is  greater  than  the  19-8,  22.5  and  2^.1  cm 
obtained  by  Kantrowitz,  but  that  is  quite  natural;  the  tire  of  the 
German  writers  was  26  cm  in  diameter,  as  against  10.2  cm  for  the 
U.  S.  tire.   It  will  new  be  attempted  to  find  the  effect  of  the  caster 
length  a  on  this  wave  length  from  equation  (26).   We  get 
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for 


a  = 

0 

W  =  36  cm 

a  = 

2  cm 

W  =  3h   cm 

a  = 

5  cm 

W  =  36  cm 

a  = 

6  cm 

W  =  39  cm 

a  = 

10  cm 

W  =  7U  cm 

a  = 

13  cm 

W  =  00 

a  < 

13  cm 

W  is  imaginary 

It  is  seen  that  the  wave  length  is  practically  constant  for  caster 
lengths  up  to  about  6  cm,  because  this  wave  length  passes  through  a 
minimum  near  a  =  2  cm.   This  explains  Kantrowitz's  second  remark: 
"the  distance  along  the  belt  required  for  one  cycle  was  also  found 
not  to  vary  much  with  caster  angle  or  caster  length." 

Unfortunately  the  author  believed  he  could  conclude  immediately 
that  the  caster  length  a  could  be  regarded  as  being  devoid  of  impor- 
tance, which  was  a  little  hasty  and  not  entirely  logical.   The  essential 
conclusion  of  the  theory  presented  here  is  exactly  to  the  contrary. 
The  effect  of  the  caster  length  a  on  the  phenomena  of  shimmy  is  of 
primary  importance. 

Furthermore,  the  preceding  tabulation  indicates  that  wave  length  W 
is  far  from  being  constant  as  Kantrowitz's  theory  assumes;  it  passes 
through  a  minimum  only  at  a  small  value  of  caster  length  a,  and  as 
soon  as  a  exceeds  6  cm,   W  is  seen  to  increase  in  such  a  way  that 
W  =  7^  cm  for  a  =  10  cm.   This  increase  is  speeded  up  more  and  more 
and  W  approaches  infinity  when  a  reaches  13  cm.   The  wave  length  for 
a  values  above  13  cm  is  imaginary,  that  is  the  track  of  the  tire  on 
the  ground  is  an  aperiodic  curve.   The  phenomenon  of  kinematic  shimmy 
has  then  completely  disappeared. 

This  fact  of  the  aperiodic ity  of  kinematic  shimmy  when  a  exceeds 
a  certain  value  is  absolutely  sure.   It  can  be  easily  checked  with  a 
caster  fitted  with  an  elastic  tire. 

Kantrowitz  's  theory  can  therefore  be  only  a  rather  rough  theory 
since  its  starting  point  is  based  on  the  assumption  of  constant  wave 
length  W,  which  definitely  becomes  infinite  when  a  reaches  values 
of  the  order  of  these  (a  >  10  cm).   As  shown  previously,  this  should 
eliminate  the  instability  of  the  oscillations,  that  is  shimmy.   In  any 
case,  the  American  author  was  absolutely  unable  to  reach  our  significant 
conclusions  regarding  the  effect  of  caster  length,  because  his  theory 
starts  by  refusing  a  priori  to  take  this  parameter  into  consideration. 

To  sum  up  the  theory  proposed  in  this  report,  fits  into  a  theory 
on  kinematic  shimmy  represented  by  equation  (25)  in  place  of  the  equation 
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i-i  +  K0  =  0 
ds2 

obtained  by  making  v  =  0  in  Kantrowitz's  equation  (lQ),  and  equa- 
tion (25)  seems  to  correspond  more  to  reality,  because  it  makes  allow- 
ance for  the  variations  of  wave  length  W  in  terms  of  caster  length  a. 

Lastly,  it  should  be  added  that  equations  (25)  and  (26)  indicate 
that  kinematic  shimmy  is  usually  a  damped  oscillation,  except  when 
a  =  0,  in  which  case  the  damping  is  zero.   The  equation 

^  +  K9   =  0 
ds2 

on  the  other  hand,  cannot  allow  for  a  damping,  whose  presence  is  easy 
to  confirm  by  experiment. 

It  will  next  be  attempted  to  compare  certain  conclusions  of 
Kantrowitz's  theory  with  the  corresponding  conclusions  of  the  suggested 
theory.   For  this  we  shall  try  to  express  the  coefficients  K,  C]_,  C2, 
beginning  with  the  characteristic  coefficients  T,  D,  S,  R,  e. 

Coefficient     K 

Kantrowitz's   differential  equation    (5)   assumes      9  =  A  sin  s yK.      The 
wave   length  of  the   kinematic   shimmy   is   then 

T 

[¥ 

But  the  present  theory  indicates  that  W  =  2n  1/ — .   Therefore 

VRe 

K  =  51  (28) 


Coefficient     C-± 

M 
Kantrowitz's   equation    (6a)   gives     C-i    =  -=-,    that   is  with  our  notations 

M  Ia 

Ct    =  -=-    and  since 
1       IA 


M  =  F(a  +  e )      and     A  =  TF 
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we  get 


Cl  =  ^  (29) 


Coefficient  C< 


Kantrowitz 's  equation  (7)  expressed  in  our  notation  reads 

^  =  -CrA     so  that  ^  =  5  =  DF  =  -D  £ 
dx    -^  dx  T 

hence 

C2  =  ° 

This  coefficient  C2  is  essentially  tied  to  the  drift  which  the 
American  author  implicitly  assumes  to  be  instantaneous. 

Numerical  results 

Kantrowitz  obtained  W  =  2l*.l  cm,  or  K  =  2-^-,  hence  K  =  0.068 

W2 
(C.G.S.  system).   Formula  (29)  indicates  that  C]_  must  be  dependent 
in  a  large  measure  on  caster  length  a.   Kantrowitz  in  fact  noted  it 
and  found 

for  a  =  0.43  cm  C]_  =  2.330  (C.G.S.  system) 

a  =  1.73  cm         Ci  =  3.1*10  (C.G.S.  system) 

Moreover,  formula  (29)  makes  it  possible  to  write 

0.1*3  +  e   2.330   I0.1*3 


(30) 


1.73  +  e   3.1*10   I 


1.73 


This  equation  permits  the  computation  of  the  length  e  by  assuming 
Io.  1*3  =  Il.73>  which  is  approximately  true  on  account  of  the  smallness 
of  a. 

Hence 

€  =  2.1*  cm 
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For  C2  Kantrowitz  obtained  C2  =  6.2  feet,  or  C2  =  0.2 
(C.G.S.  system). 

Calculation  of  T,  D,  S,  R,  and  € 

Coefficient   e   has  already  been  defined.   The  other  characteristic 
coefficients  of  the  tire  studied  by  Kantrowitz  can  be  defined  by  the 
formulas  (28),  (29),  and  (30),  in  which  the  inertia  I  is  involved. 
Kantrowitz  gives   I  =  1.06  X  10-^  in  slugs  and  feet.   The  slug  is  a 
mass  unit  of  Ik. 6   kg  and  1  foot  =  30. k8   cm,  so  that 

I  =  l.l+lj-0  C.G.S.  units 
Formula  (29)  gives  then 

a  +  g    0.^3  +  2.k 


T  = 


I-1    I.Wj-O  X  2330 


or 

T  =   8k  X  10-8  in  C.G.S.    units 
Formula   (28)   gives 

KT       0.068  x  Ok  x  IP-8 

G  2.k 

or 

R  =   2.38  X   10-8   in  C.G.S.    units 

likewise   formula    (30)    gives 


^.2 

or 


D  =  C0T   =   0.2  X   8k  X   10' 


•  8 


D  =  16.8  X  10_0  in  C.G.S.  units 
Lastly, 

S  =  7  x  10~8  in  C.G.S.  units 

by  reason  of  the  formula  Se  =  D. 
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The  kinematic  shimmy  of  Kantrowitz 

By  the  use  of  the  derived  characteristic  coefficients  together 
with  formula  (26)  the  determination  of  the  wave  lengths   W  of  kinematic 
shimmy  made  on  the  Schlippe-Dietrich  tire  can  be  applied  to  Kantrowitz's 
tire. 


For 


a  =  0 

W  =  2k   cm 

a  =  2  cm 

W  =  21  cm 

a  =  h   cm 

W  =  27  cm 

a  =  6  cm 

W  =  38  cm 

a  =  11  cm 

W  =  00 

a  >  11  cm 

W  is  imag 

This  tabulation  proves  that  the  conclusions  valid  for  the  German  tire 
are  even  more  valid  for  the  American  tire. 

Inversion  velocity  u 

Having  defined  the  characteristic  coefficients  the  inversion 
velocity  u  defined  previously  by  the  equation 


2    a 
u  =  — rr 
IE 


can  now  be  computed. 

Since  a  high  caster  angle  is  likely  to  modify  this  u,  we  shall  not 
make  the  calculation  for  the  test  series  for  which  this  high  caster  angle 
was  20  but  only  for  the  series  where  it  was  5°,  which  is  small  enough 
to  be  regarded  as  zero.   Then 


u2  = 


0.^3 


i.Mk)  x  2.38  x  10"8 

hence 

u  =  112  cm/sec 
that  is 

u  =  3.70  ft /sec 
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But  on  considering  figure  30  and  the  curve  A  of  the  convergence 
(for  a  caster  angle  of  5°  and  a  =  0.  i+3  cm)  it  is  plain  that  the  exten- 
sion of  the  experimental  curve  does  not  pass  through  the  origin.   This 
circumstance  is  also  relevant  for  the  curves  of  figures  20  and  21 
established  by  Wylie.   The  extension  of  the  experimental  curve  A 
plotted  by  Kantrowitz  exactly  intersects  the  axis  of  the  velocities  in 
the  axis  of  the  abscissa  v  =  3-7  ft/sec;  as  stipulated  by  the  theory. 

In  contrast  to  Kantrowitz's  theory,  this  is  an  additional  accom- 
plishment of  the  present  theory.   In  fact,  it  will  be  shown  that  by 
Kantrowitz's  equation  (10)  the  curve  X,(v)   passes  through  the  origin 
which  does  not  correspond  at  all  to  the  experimental  results. 

In  order  to  make  the  comparison  between  the  two  theories  more 
accurate,  consider  Kantrowitz's  equation  (10)  and  compare  it  with  equa- 
tions (6)  and  (22)  into  which  the  suggested  theory  fits.   First  of  all, 
equation  (10)  is  transformed  by  taking  the  time   t  as  independent 
variable  in  place  of  the  distance   s,  and  then  the  amplitude   z   com- 
puted in  place  of  the  angular  elongation  9.      On  the  one  hand 

G  __   z        d9_  ldz   d££  _  1  d£z   d_3_£  __   1  d_3_z 
.   "a   ds  ""  a  ds   ds2  "  a  ds2   dg3  "  a  ds3 

on  the  other 

,    dz    1  dz 

s  =  vt   —  = 

ds   v  dt 


d2z 

1  d2z 

d3z 

1  d3z 

ds2 

v2  dt2 

ds3 

v3  dt3 

so  that  equation  (10)  becomes 


Cxv  dt3 


z  =  0 


that  is 


(a  +  e )v  dt3 


(3D 


by  expressing  the  constants  K,   C]_,  C2  by  equations  (28),  (29),  and 
(30).   We  repeat  equation  (6)  of  the  simplified  theory 
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IT 


d35 


dt" 


ID 


(a  +  e)v  ,.3   a  + 


d2z 


dt 


a  dz 
+ +  z 

v  dt 


(6) 


TSa 

and  equation  (22)  of  the  complete  theory  after  multiplying  by  , 

a  +  e 

taking  into  account  the  relationship  Se  =  D 


ITS 


d^z 


IT 


d^z 


R(a  +  e)v2  dt^   (a  +  e)v  dt3 


ID 


a  +  e 


T  +  Sa£ 


R(a  +  e)vf_Jdt 


d2z 


a  dz 

+ +  z 

v  dt 


(22) 

Now  equation  (31)  can  be  compared  with  equations  (6)  and  (22).   Inci- 
dentally, the  turn  coefficient  R   (consequently  the  phenomenon  of 
torsion)  was  implicitly  introduced  in  Kantrowitz's  theory,  although 
its  author  claims  to  have  ignored  the  torsion  of  the  tire.   Kinematic 
shimmy  can  only  be  explained  by  the  effect  of  this  torsion  and  it  is, 
in  fact,  easy  to  prove  that  equation  (6),  which  would  be  a  complete 
theory  if  this  phenomenon  of  torsion  could  be  disregarded,  cannot  give 
an  account  of  kinematic  shimmy.   A  second  preliminary  remark  is  needed. 
When  the  linear  differential  equation  of  the  third  order  with  positive 
coefficients 

Az'"  +  Bz"  +  Cz'  +  D  =  0 

is  considered,  Routh's  rule  is  summed  up  as  follows:   for  the  oscilla- 
tion to  be  convergent  it  is  sufficient  and  necessary  that  BC  >  AD.   In 
equation  (31)  C  =  0;  therefore  the  oscillation  can  never  be  convergent. 
As  a  result  Kantrowitz  could  not  be  led  to  suspect  the  possibility  of 
combating  shimmy  by  a  simple  judicious  combination  of  parameters.   This 
was  unfortunate  since  it  led  to  the  design  of  landing  gears  having  as 
short  a  caster  length  as  possible.   Lastly,  in  comparing  equation  (6) 
of  the  simplified  theory  with  equation  (22)  of  the  complete  theory  it 
is  seen  that  the  only  difference  is  the  presence  in  equation  (22)  of 
two  supplementary  terms: 


ITS 


d^z 


R(a  +  e)v2  dt1* 


and 


T  +  Sa2   d2z 
R(a  +  e)v2  dt2 


Each  one  of  these  terms  carries  the  factor  v2  in  the  denominator.   So, 
when  the  velocity  increases  infinitely,  the  motion  represented  by 
equation  (22)  tends  asymptotically  towards  the  motion  represented  by 
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equation  (6).   In  other  words,  when  the  velocity  is  large  enough  the 
two  equations  (6)  and  (22)  can  be  regarded  as  practically  equivalent. 
So,  for  a  comparison  of  equations  (22)  and  (31)  at  high  velocities, 
it  is  sufficient  to  compare  equations  (6)  and  (31). 

We  shall  now  define  the  relationship  existing  between  velocity  v 
and  divergence  X     in  the  case  of  equation  (31)  (that  is,  in  the 
Kantrowitz  theory),  as  we  did  before  for  equation  (6). 

Employing  the  same  method  as  for  equation  (ll),  the  system 
IT    w,  p   „  os        /ID      T 


XQ,2  _  30,2)  +  _t^_   +  JL_VX2  -  a£)  +  1  =  0 


(a  +  e)v  \a  +  e   Rev2 


(3X2  _  ^)  +  (_^_  +  _LJ\2X  - 


(a  +  e)v 


\a  +  €   Rev2y 


is  obtained.   To  obtain  the  relationship  looked  for  between  v  and  X, 
simply  eliminate  cu2  from  the  two  preceding  equations.   Thus 

8IT    X3  +  s(™-  +  ^-V2  ♦  2(a  +  6)v/lD_  +  _T_\\  =  x   (32) 
(a  +  e)v       Va  +  e   Rgv2/         IT   \a  +  e   Rev2y 


This  equation  assumes  that  time  t   is  the  independent  variable.   Now, 
in  his  calculations  as  well  as  in  the  curves  of  figures  30  and  31, 
Kantrowitz  has  used  the  distance  x  -  vt  as  the  independent  variable. 
The  fundamental  oscillation  of  shimmy 

z  =  e^sin  tut 
is  then 

z  =  e^sin  2  x 
v 


by  putting 


\ 

u  =  - 
v 


To  obtain  the  equation  giving  the  new  expression  of  the  divergence,   A. 
is  simply  replaced  by  vu   in  equation  (32).   We  then  get,  after 
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arrangement  of  the  terras  of  this  equation 


+  e\  V     \Re 


2Iu  £    D  \  k   /  8   b   UD      Up   2  (a  +  e) 

— -±—  (2u  +  -  )  v4  +  (  —  u^  +  — —  u  -  i)v^  +  -^ i-  u  =  0 

a  +  e  V     T/     V^e      TRe     V       IR2e2 


With  respect  to  velocity  v  this  is  an  equation  of  the  fourth  power, 
whose  two  primary  roots,  if  they  are  real,  have  the  same  sign.   In 
order  that  these  roots  are  real,  the  discriminant  must  be  positive. 
Therefore  after  some  simplifications  we  obtain 

l6Tu2  +  8Du  -  Re  ^  0 

which  is  true  when  u  has  a  value  between  U]_  and  up  of  this 

trinomial.   These  roots  being  of  opposite  sign,  it  is  sufficient  that 
u  is  less  than  the  positive  root,  hence 


^D2  +  TRe  -  D 

u  <  1 

kT 

It  is  readily  apparent  that,  when  this  condition  exists,  the  roots 
are,  of  necessity,  positive,  because  the  coefficient  of  v2   is  then 
negative.   In  fact,  this  coefficient  is  a  trinomial  of  the  second  degree 
in  u  allowing  two  roots  of  opposite  sign,  hence  the  trinomial  is 
negative  when  u  is  less  than  the  positive  root.   But  this  root  is 


\w7 


2TRe  -  D 
^1  =  ti 

and  it  is  obviously  greater  than  the  positive  root  of  the  discriminant. 
Thus  the  first  condition  is  sufficient  to  assure  the  second,   u  <  u]_. 
The  result  is  that  the  divergence  u  passes  through  a  maximum 

\/d2  +  TRe 

*- W~ 

for  a  given  velocity  vm  given  by  the  equation 

2      a  +  e 


v, 


m  - 

IRe 


e + f) 
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Thus,    for  the   tire   studied  by  Kantrowitz    (curve  A    in   fig.    36) 

um  =   0.0322     C.G.S.    units 
v  =   36^  cm/sec 
or    in  English  units    of  feet    (30.5  cm) 

um  =  0. 98     English  units 

vm  =   11.9  ft/sec 

Shape  of  curve   X(v)   at  high  velocities 

The  variation  of  A-(v)  when  the  velocity  increases  infinitely, 

will  be  examined  by  putting  A,  =  —  and  trying  to  define  K  in  terms 

of  the  characteristic  coefficients.   Entering  this  value  in  the  preceding 
equation  while  ignoring  the  infinitely  small  terms  gives  the  relationship 


2(a 


+  Qv/ID  \2  K  = 
IT    \a   +  J     v  " 


so  that 

T(a  +  e ) 


K  = 


2D2I 


The  same  calculation  with  equation  (ll)  gives 

(T  -  aD)(a  +  O 


K'  = 


2D2I 


Thus  it  is  seen  that  the  present  theory  produces  lower  values  for  X 
at  high  velocities  than  Kantrowitz 's  theory  and  seems  to  explain  to 
some  extent  the  remark  made  by  him  that  the  values  deduced  by  his 
formula  are  clearly  too  high  for  high  velocities. 

To  complete  the  comparison  of  Kantrowitz 's  theory,  of  the  preceding 
theory  and  of  the  experimental  results,  the  calculations  of  X     and   u 
at  various  velocities  were  made 

(l)  By  Kantrowitz 's  theory  (equations  (32)  and  (32a)) 
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(2)  By  the  present  theory  in  its  approximate  form  (equation  (12) 

(3)  By  the  present  theory  in  its  complete  form  (equation  (22)) 


Kant 

P: 

resent 

theory 

rowitz  's  theory 

Approximate 
form 

Co 

mplete  form 

V 

cm/sec 

X 

C.G.S. 

C.G.S. 

English 
unit 

X 

C.G.S. 

X 
C.G.S. 

C.G.S. 

English 
unit 

0 

0 

0 

100 

.87 

0.0087 

0.265 

112 

0 

0 

0 

200 

4.77 

.0238 

•  730 

10.5 

.0525 

1.6 

294 

Xm  =  24.7 

300 

9-36 

.0312 

•  95 

19.5 

.0650 

1.98 

364 

11.7 

.0322 

•  98 

400 

12.7 

Hm-0318 

Um-97 

20.6 

.0515 

1.57 

500 

14.8 

.0296 

•  90 

20.8 

.0415 

1.27 

6oo 

16.0 

.0267 

.815 

20.8 

.0346 

1.05 

700 

16.6 

.0238 

.725 

8oo 

16.71 

.0209 

•  635 

19.6 

.024 

■  73 

900 

16.65 

.0185 

•  555 

1000 

16.4 

.0164 

•  50 

17.6 

.0176 

.54 

1500 

13.62 

.0091 

.28 

15.3 

.0102 

.31 

2000 

11.48 

. 00575 

.175 

11.3 

.00565 

.17 

2500 

10 

.004 

.12 

10 

.004 

.12 

The  table  is  graphically  represented  in  figure  32  with  u  plotted 
against  the  velocity  v  for 

(1)  Experimental  data  (curve  E) 

(2)  Data  from  Kantrowitz  theory  (curve  K) 

(3)  Data  from  present  theory  in  its  complete  form  (curve  C) 

The  velocity  is  expressed  in  cm/sec  and   u   in  English  units  used  by 
Kantrowitz.   To  obtain  curve  E,  plot  the  15  experimental  points  given 
by  Kantrowitz,  then  draw  the  most  likely  mean  curve  defined  by  those 
15  points. 
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The  plotting  of  curve  K  involves  the  resolution  of  a  linear 
differential  equation  of  the  third  order  which  presents  no  special 
difficulties,  as  the  necessary  calculations  can  be  made  by  any  of 
several  classical  methods;  the  plotting  of  curve  C  reduces  to 
resolving  equation  (22)  which  is  a  linear  differential  equation  of  the 
fourth  order.   It  is  a  more  difficult  problem  and  involves  considerable 
paper  work,  so  a  graphical  method  of  producing  quick  and  excellent 
results  is  desired. 

Graphical  method  for  solving  linear  and  homogeneous  linear  equations 

The  equation  (22)  that  is  to  be  solved,  is  of  the  form 

Az17  +  Bz111  +  Cz11  +  Dz1  +  Ez  =  0 

which,  putting  z  =  Kee   ,  gives  the  characteristic  equation 

f  (x)  =  Ax^-  +  Bx3  +  Cx2  +  Dx  +  E  =  0 

Sought  are  the  real  numbers  or,  complex  x,  solutions  of  this  equation. 
Consider,  therefore,  the  complex  plane  0,  X,  Y,  (fig.  33)-   To  every 
point  m  of  the.  prefix  x  the  preceding  equation  makes  correspond  a 
point  d'  with  prefix  f(x).   This  point  M   is  easily  obtained  from 
point  m  by  simple  graphical  operations  with  rule  and  dividers. 
Point  m  corresponds  to  a  root  x  "hen  M  reaches   0.   Then  consider 
a  point  m1  near  m  to  which  corresponds  a  point  M'  near  M,  and 
which  gives  by  definition 


MM'     ,.  , 
=>  =   f'(x) 


However,  the  value  of  f*(x)   is  not  dependent  on  the  position  of  m'; 
but  solely  on  that  of  m.   Therefore,  if  two  vectors  Pm'*  and  PM  •  * 
equipollent  to  mm1'  and  MM*  are  involved,  the  triangle  Pm'^" 
remains  similar  regardless  of  the  position  of  point  m'  near  m.   Con- 
sequently, in  the  only  condition  where  the  prefix  x  of  m  is  not  one 
of  the  roots  of  the  equation   f '(x)  =  0   (which  would  be  an  exceptional 
case  since  it  is  impossible  to  have  more  than  three  points  in  the  entire 
plane  satisfying  f'(x))  a  point  m*  near  m  is ,  as  a  rule,  easily 
found,  so  that  its  transformed  M'   is  closer  to  0  than  point  M. 
After  a  sufficient  number  of  operations,  the  point  M  can  be  gradually 
brought  as  near  to  point  0  as  desired.   This  is  the  principle  of  a 
graphical  solution  of  algebraic  equations  by  successive  approximations. 


NACA  TM  1337  59 


It  should  be  noted  in  passing  that  the  remark  contains  implicitly  the 
rigorous  proof  of  d'Alembert's  theorem  which  states  that  any  algebraic 
equation  with  real  or  imaginary  coefficients  has  always  at  least  one 
real  or  imaginary  root.   In  many  cases  the  foregoing  operations  can  be 
made  to  converge  very  quickly.   The  method  consists  in  taking  as  point  m 


the  second  point  m]_  so  that 

MO* 


f'(x)  =^3 


mm-i  mm 

M  and  M1  being  constructed  on  the  line  of  m  and  m1,  only  m' 
being  chosen  near  m  on  the  radius  OM.   The  graphical  method  of 

MM' 
defining  the  point  m^  by  utilizing  the  ratio  — y     is  much  more  simple 

mm' 

than  trying  to  determine  m^  directly  from  the  derivative   f'(x). 

After  obtaining  mj_  deduce  M]_,  then  begin  again  with  torque  m^M^, 

the  operations  to  be  made  on  torque  mM.   The  result  is  a  torque  n^Mg. 
As  a  rule,  the  successive  points  M,  M]_,  Mg,  •  •  .   come  closer  and 
closer  to  0  and  the  points  m,  m^_,   nip,  •    •    •   tend  to  a  limiting  point 

whose  prefix  X  +  ico  is  the  root  of  the  equation  f(x)  =  0.   The 
function 

z  =  e^sin  cut 

is  then  the  corresponding  solution  of  the  differential  equation  that 
was  to  be  solved. 

The  unusual  feature  of  this  method,  whose  principle  is  apparent 
from  Newton's  method  for  finding  real  roots  of  algebraic  equations  is 
the  small  number  of  operations  necessary  to  obtain  convergence  as  soon 
as  a  root  is  approached.   This  graphical  method  is  absolutely  general 
and  can  be  applied  to  algebraic  equations  of  any  degree. 

After  this  digression  the  study  of  the  three  curves  in  figure  32 
is  resumed. 

The  comparison  of  the  experimental  curve  E  with  Kantrowitz 's 
theoretical  curve  K  has  shown  th?t  the  maximum  of  the  latter  curve 
is  not  clearly  enough  indicated.   Curve   K  is  much  too  flat  and  drops 
too  slowly  when  the  velocity  increases.   Specifically  curve  K  gives 
slightly  too  small  divergences  near  the  maximum  and  by  way  of  contrast, 
definitely  too  great  divergences  at  high  velocities,  as  the  author 
himself  admits. 
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The  comparison  of  the  experimental  curve   E  with  curve  C   given 
by  the  present  theory  indicates  that  the  latter  gives  value  sensibly 
too  great,  but  that  the  general  shape  of  the  curve  C   (existence  of 
the  same  inversion  velocity  u,  acuteness  of  the  maximum,  rapidity  of 
decrease  at  increasing  velocity)  approaches  curve  E'  more  closely 
than  the  flat  shape  of  curve  K.   The  correspondence  between  curve  C 
and  curve   E   is  almost  a  simple  expansion  of  the  ordinates  in  a  ratio 
approaching  1  as  the  velocity  increases.   So  in  spite  of  the  similarity 
of  forms  between  E   and  C   there  still  is  a  certain  discrepancy 
between  the  experimental  and  theoretical  results  of  the  present  theory. 
It  is  the  first  discrepancy  of  this  kind  that  has  been  found.   A  careful 
examination  to  see  whether  this  can  be  explained  within  the  framework 
of  the  present  theory,  was  therefore  indicated. 

It  seems  preferable  to  have  the  theoretical  discrepancies  on  the 
high  side  rather  than  the  low  (as  in  the  case  of  curve   K  in  the 
neighborhood  of  the  maximum).   The  divergence  of  the  oscillations 
corresponds,  in  fact,  to  a  storage  of  mechanical  energy  in  the  castering 
wheel  and  the  aim  of  the  theory  of  shimmy  is  precisely  to  discover  the 
source  and  the  laws  of  this  energy.   Shimmy  is  accompanied  by  various 
energy  dissipating  phenomena  for  which  the  theory  makes  no  allowance, 
such  as  friction  in  the  castering  axis,  the  slight  slippage  of  the  tire 
close  to  the  periphery  of  the  contact  surface  and  particularly  of  the 
tail  point,  the  reactive  effects  at  the  inside  of  the  tire,  etc.   The 
last  cause  is  far  from  being  negligible:   according  to  the  general 
scheme  of  oscillatory  motions  shimmy  is  accompanied  by  a  periodic  and 
continuous  transformation  of  kinetic  energy  of  the  wheel  into  potential 
stress  energy  of  the  tire  and  vice  versa.   It  is  known  that  a  stressed 
tire  recovers  only  part  of  the  stress  energy.   Drop  tests  on  tire  indi- 
cate, for  example,  the  height  of  rebound  is  always  from  15  to  30  percent 
less  than  the  height  of  drop,  the  loss  being  due  largely  to  the  reaction 
effects  inside  the  tire.   It  is  therefore  to  be  expected  that  in  a 
correct  theory  of  shimmy  the  values  of  divergence  are  always  greater 
than  the  experimentally  observed  values.   Still,  this  explanation 
hardly  justifies  a  difference  exceeding  10  to  20  percent.   The  difference 
between  the  maximum  point  of  the  theoretical  curve  C  and  the  experi- 
mental point  amounts  to  -^-  =  ^0  percent.   There  must  be  some  other 
^  20      ^ 

cause  for  explaining  such  a  marked  discrepancy.   There  is  only  one 
possibility,  namely,  the  numerical  values  used  for  T,  D,  S,  R,  and  e 
are  probably  not  rigorously  exact.   Unfortunately  these  factors  had 
been  defined  indirectly  from  Kantrowitz 's  factors   K,  C]_,  and  C2,  hence 

there  is  a  possibility  that  some  of  the  defined  factors,   T,  D,  S,  R, 
and  e   were  only  approximate. 

This  matter  will  be  discussed  in  a  little  greater  detail:   For- 
mulas (28),  (29),  and  (30)  were  taken  as  a  starting  point.   Formula  (29) 
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resulted  in 

0.1+3  +  e  .  2330  T0.k3 
0.73  +  e    3^10  I1-73 

which  made  it  possible  to  compute   e  =  2.k. 

Obviously,  this  value  for   e   is  only  approximate,  and  for  two 
reasons:   first,  it  was  assumed  that 

I0.1i3  =  I1.73 

although  the  second  of  these  values  is  certainly  greater  than  the  first. 
On  the  other  hand,  formula  (29)  makes  no  allowance  for  inclination  of 
caster.   Rigorously  written  it  should  have  been 


9.1  = 


M    F(a  +  e)  +  p9 


IA        IA 
with  A  =  TF  and  6  =   DF,  that  is 


a  +  g   PD 
u-i  = +  — 

— L      jT       IT 


the  parameter  p   being  defined  later.   However,   T  was  computed  by 
formula  (29)  from  the  derived  value  for  e,  which  is  the  reason  that 
the  value  obtained  for  this  parameter  cannot  be  entirely  rigorous. 
Then  R  was  computed  by  formula  (28).   This  formula,  utilizing  the 
length  of  kinematic  shimmy,  seems  to  be  rigorous.   Unfortunately,  it 
contains   €  and  T  which  limits  the  accuracy  to  be  expected  from  R. 
Lastly,   D  was  computed  by  formula  (30)  and  S  by  the  relation  Se  =  D. 
Our  formula  (30)  is  simply  Kantrowitz 's  formula  (7),  which  is  equivalent 
to  assuming  that  drift  is  instantaneous.   And  this  assumption  is,  as 
already  shown  in  section  III-B,  only  a  simplifying  assumption  frequently 
permitted,  but  which,  in  this  instance,  is  likely  to  taint  the  calcula- 
tion of  CQ  with  an  error. 

A  rigorous  determination  of  the  equation  of  the  curve  studied  by 
Kantrowitz  in  figure  28  will  now  be  attempted.   This  curve  was  used  to 
determine  the  coefficient  Co. 
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Let   y  represent  the  distance  of  the  center  M  of  the  contact 
tread  with  respect  to  the  undeflected  plane  of  the  wheel.   With  our 
usual  notations  (fig.  3^-) 

d2y 

— *-  =   -KM 

dx2 


M  =  -  -  Fe 

S 


ay 

a  =  — 
dx 


1,     ,    y 


F  =  -(z  -  y)  = 

T  T- 


So  the  final  equation  after  eliminating  M,  a  ,  and  F  from  the  pre- 
ceding equations  reads 

R  dx2   S  dx   T 


The  most  general  solution  of  this  differential  equation  is  of  the 
following  form 

1 

S-iX       SoX 

y  =  Ae  x  +  Be  ^ 


A  and  B  are  two  coefficients  depending  on  the  initial  conditions 
and  S]_  and  sp  the  two  roots  of  the  equation  of  the  second  degree 

Mi  +  s  +  1   =  0 
R    S   T 

Therefore,  the  equation  of  curve  L,   which  represents  the  path  of  the 
center  M  of  the  contact  tread  (fig.  3*0  is 

GiX     SoX  -CqX     -  mx 

y  =  Ae  x  +  Be  ^      and  not  y  =  Ce  ^    =  Ce  l 
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as  assumed  by  the  American  author  without  adequate  reason.   As  a  con- 
sequence the  parameter  D  computed  by  wrongfully  assuming  that  it  is 
the  equation  of  curve   L,  may  very  well  be  tainted  with  an  appreciable 
error,  which. by  itself  might  be  sufficient  to  explain  the  discrepancy 
between  curve   E  and  C   in  figure  32. 


V.  EFFECT  OF  RESTORING  TORQUE 
(inclination  of  Castering  Axis) 

EFFECT  OF  FLUID  AND  SOLID  FRICTION 
(Damping  of  Shimmy) 


A.  Effect  of  Restoring  Torque 


In  the  foregoing  it  has  been  assumed  that  the  castering  axis  was 
vertical.   The  next  logical  question  is:   what  is  the  effect  of  an 
inclination  of  the  castering  axis  on  the  stability  of  the  vibrations? 
It  is  readily  seen  that,  provided  it  is  not  too  great,  such  an  inclina- 
tion should  be  equivalent  to  the  existence  of  a  restoring  torque  propor- 
tional to  the  angle   0,  this  torque  being  moreover,  negative  in  the 
sense  of  figure  25,  that  is  tending  to  move  the  wheel  from  its  mean 
position  0=0. 

Now  the  general  problem  of  the  effect  of  a  restoring  torque  Cr  =  p0 
will  be  analyzed.   To  express  the  problem  in  equation  form,  the  same 
line  of  reasoning  used  for  the  system  of  four  equations  (l6),  (17),  (18), 
and  (19)  is  followed.   It  is  apparent  that,  in  these  conditions,  equa- 
tions (l6),  (17),  and  (19)  are  not  changed  by  the  existence  of  a 
restoring  torque.   Equation  (l8)  can  be  replaced  by 

I  — ?  =  -F(a  +  e)  -  M  -  p9 
dt^ 

or 

r<£|=  _F(a  +  e)  _M  _  p  z 
a  dt2  -    a 

since   0  =  — .   To  the  moments   -F(a  +  e)   and   -M  produced  by  the 

reactions  of  the  ground  on  the  contact  tread,  must  be  added  the  moment 
-p0  of  the  restoring  torque,  p  being  positive  when  a  true  restoring 
torque  is  involved  and  negative  when  a  torque  of  opposite  effect  is 
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involved  (as  for  example  in  the  case  of  a  disposition  of  the  caster  as 
illustrated  in  fig.  25). 

The  result  is  the  equation  system 

±  &  =  -m  (16) 

v^  dtd 


z  -  y  =  TF  (17) 


I  Ll=   _F(a  +6)  -  M  -  £  z  (18a) 

a  dt2  "   a 


M  =  -i-  ^  +  —  -  F€  (19) 

Sv  dt   Sa  v 


The  values  of  F  and  M  of  equations  (17)  and  (19)  entered  in  equa- 
tions (l6)  and  (l8a)  give  the  two  equations 


1   d2y  ^ 


.  JL  &  +  1  y  +  (1.   .  l\B  =  0  (20) 

Rv2  dt2   Sv  dt   T     Vsa   T^ 

-L  *Z  -  2:  y  +  I  £5.  +  f £  +  -L  +  al  =  0  (21a) 

Sv  dt   T  J   a  dt2   \a   Sa   T/ 

The  values  of  z  from  equation  (20)  posted  in  equation  (21a)  give  the 
differential  equation 

_I_  dV  +  _i_   d3y  +  /][*  +  _1_  +  _a_  +  _P_\afy  + 
Rav2  dt^   Sav  dt3   \Ta   SEav2   TRv2   Rav2/dt2 

/a_i±  +  _£_W  +  (!L±±  +  iEV  =  0  (33) 

\  TSv    Sav/dt   \  TSa    aT/ 

It  is  easily  seen  that  the  second  remark  made  below  equation  (22)  is 
still  applicable  here.   Therefore,  to  obtain  the  differential  equation 
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of  motion  of  the  wheel  y  Is  simply  replaced  by  z  in  the  preceding 
equation.  The  result  is,  then,  a  linear  differential  equation  of  the 
fourth  order  of  the  form 

ar,z   +  anz    +  aQz   +  a5z   +  ai.z  =  0 


In  connection  with  equation  (22)  it  was  shown  that  Routh's  conditions 
for  stable  oscillation  are  as  follows 


ala2  >  a0a3 


o     p 
a]_aoao  >  aQao^  +  a-^  a^ 


Applied  to  equation  (33)  the  first  of  these  conditions  gives,  after 
reductions 

Ie     T      g 
—  + >  -i— 

a    SRav2   Rv2 


or,  if  the  fundamental  relation  Se  =  D  is  taken  into  account 

p   aD  -  T  ,ol  , 

V  >  "iRD-  '         (3i° 

The  second  of  the  preceding  equations  gives,  after  reductions 


Ta   SRav2   SRa2v2(a  +  e)   TRv2   Sa2   Rav2(a  +  e 


or 


IRv2  -  a )(Sea  -  T)  >  0 

a  +  e/ 
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that  is  if  the  fundamental  relationship  is  taken  into  consideration 

ZlRv2  -  a  -  — — J(aD  -  T)  >  0  (3^a) 

To  complete  the  discussion  of  these  stability  conditions,  two  cases  are 
distinguished  according  to  the  sign  of  the  factor  p. 

(a)  When  a  true  restoring  torque  is  involved,  that  is  when  p  >  0, 
the  discussion  presented  following  equations  (23)  and  (23a)  applies 
here  also.   Consider  the  velocity  u  defined  by  the  equality 

IRu2  =  a  +  -PT 


a  +  e 


First  case:   aD  >  T.   The  condition  (3*+a)  is  proved  when  v  >  u. 

Second  case:   aD  <  T.   The  condition  (3^-a)  is  proved  when  v  <  u. 

It  is  plain  that  in  either  case,  when  equation  (3*t-a)  is  proved 
equation  (3*0  is  also  confirmed.   Therefore,  the  conditions  of  stability 
of  oscillation  are 

First  case:  aD  >  T     with  v  >  u 

Second  case:   aD  <  T    with  v  <  u 

The  problem  therefore  is  like  that  discussed  following  equation  (22) 
and  it  is  clear  that  the  effect  of  a  true  restoring  torque  on  the  sta- 
bility conditions  is  simply  an  increase  in  the  inversion  velocity  u. 

(b)  When  a  negative  torque  is  involved,  that  is  when  p   is  nega- 
tive, which  is  reached  at  an  inclination  of  the  caster  as  shown  in 
figure  25,  it  can  no  longer  be  asserted  that  the  confirmation  of  equa- 
tion (3^-a)  automatically  entails  the  confirmation  of  equation  (3*0. 

It  is  seen  that,  if  aD  <  T  this  condition  (3*0  is  always  proved, 
but  if  aD  >  T,  the  condition  (3**-)  is  proved  only  when  v  >  w,  w 
denoting  the  velocity  defined  by  the  equality 

IRw2  =  a  -  £ 
D 

The  stability  conditions  of  the  oscillation  are  therefore 
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First  case:   aD  >  T  with  v  >  u  and  simultaneously  v  >  w 
Second  case:   aD  <  T  with  v  <  u 

On  trying  to  compare  u  and  w  it  is  readily  apparent  that  in 

a  +  6 
order  that  u  >  w  it  is  necessary  and  sufficient  that  -p  >  — - — . 

From  the  previously  defined  values  for  e   and  D  and  the  value  of  p 
that  will  be  computed  subsequently  (formula  35)  it  follows  that,  except 
for  large  inclination  of  the  axis  and  substantial  wheel  loads,   u  >  w, 
regardless  of  the  caster  length  a. 

The  discussion  is  again  similar  to  that  relevant  to  equation  (22) 
to  the  effect  that  a  negative  restoring  force  on  the  stability  condi- 
tions is  simply  a  decrease  of  the  inversion  velocity  u. 

Let  us  now  refer  to  figure  30  and  compare  the  experimental  curves  A 
and  B  given  by  Kantrowitz.   By  extending  these  curves  up  to  their 
intersection  point  with  the  velocity  axis  the  corresponding  inversion 
velocities  are  obtained.   Obviously  the  inversion  velocity  relative  to 
curve   B   is  substantially  lower  than  that  with  respect  to  curve  A. 
But  curve  B  is  exactly  that  which  has  been  obtained  for  a  20  degree 
inclination  of  the  caster  axis,  corresponding  to  a  negative  restoring 
torque.   This  experimental  check  although  purely  qualitative  (because 
the  author  failed  to  mention  the  load  applied  to  the  wheel,  and  which 
is  needed  for  defining  the  factor  p)  is  an  added  proof  of  the  present 
theory. 

The  decrease  in  u  following  the  inclination  of  the  caster  axis 
is  much  more  evident  in  this  example  since  the  experiments  relative  to 
the  curve  B  of -figure  30  had  been  made  with  a  caster  length  (a  =  1.73  m) 
greater  than  that  relative  to  curve  A   (a  =  0.1+3  m).   So  without  incli- 
nation the  inversion  velocity  would  have  been  greater,  and  it  actually 
would  have  resulted  in 


u  = 


/I.73  X  108         OOK        / 
\ — -1^ =  225  cm/sec 

1.^+0  X  2.38 


for  curve  B,  instead  of  u  =  112  cm/sec  for  curve  A. 

Calculation  of  coefficient  p   for  an  inclination  cp  of  the  spindle 
axis. 

For  the  case  of  an  inclination  cp  of  the  spindle  axis,  the  coef- 
ficient p  of  the  negative  restoring  torque  can  be  defined  in  terms 
of  the  load  P  of  the  wheel  on  the  ground,  the  caster  length  a,  and 
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the  angle  cp.   Consider  figure  35  which  represents  a  tire  touching  the 
ground  at  M  and  assume  that  AP  is  the  spindle  axis.   When  the  wheel 
turns  about  AP   it  can  be  assumed  for  first  approximation  that  M 
describes  a  circle  centered  on  AP,  that  is  a  circle, of  center  H. 
Axis  AP   is  assumed  to  be  stationary  while  the  wheel  turns  through  a 
small  angle  9     around  this  axis.   Now  it  is  readily  seen  that  the 
point  M   is  raised  up  to  a  height 

—  a2  g2 

h  =  MH  —  sin  cp  =  a  —  sin  cp  cos   cp 
2  v  2  ^  v 

The  work  performed  by  the  force  P  is  therefore 

92 
Ph  =  Pa  —  sin  cp  cos  cp 

with  Cr  =  p0  representing  the  negative  restoring  torque  it  is  evident 

that  the  foregoing  work  is  equal  to  that  performed  by  the  torque  Cr, 
that  is 

Cr  d©  =  /   p0  d9  =  p  i 


—T 


0 


2 


Comparison  of  this  expression  with  the  preceding  one,  gives  the  sought 
for  expression  of  the  coefficient  p 

p  =  -Pa  sin  cp  cos  cp  (35) 

So,  if  the  American  author  had  not  forgotten  to  give  us  the  pres- 
sure P  of  the  wheel  on  the  ground,  we  would  have  been  able  to  use  his 
experimental  results  for  a  numerical  check  of  the  present  theory.   In 
the  absence  of  these  data  the  foregoing  considerations  can,  however, 
still  be  utilized  for  a  numerical  calculation  by  adopting  an  inverse 
method  of  looking  for  the  evaluation  of  P  from  an  examination  of  the 
curve  B. 

From  an  examination  of  the  curve  B  in  figure  30  it  is  found  that 

the  inversion  velocity  u  is  reduced  to  a  value  approaching  zero,  by 

the  inclination  of  the  spindle  axis.   Hence  we  assume  u  =  0  It  was 
seen  that 

2    a       PT 

IR   (a  +  e)lR 
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hence 


a   Pa  sin  cp  cos  cpT 
IR      (a  +  e)lR 


consequently 


a  +    e                                  1.73   +  2.k  ,co  w   ,,.5    , 

P   =  =  ^ =    153  X   ICk   dynes 

T  sin  cp  cos   cp        81+  X   10"8  X  0.31+2  X  O.9U 
that   is 

P  #  15. 3  kg  weight 

So,  for  the  first  series  of  experiments  (fig.  30,  curves  A  and  A') 
we  get 

p  =  _io7  x  0.1*3  X  0.087  x  O.996  #  -  3.5  X  105  C.G.S. 
and  for  the  second  (curves  B  and  Bf) 

p  =  -10?  X  1.73  X  0.3^2  X  O.9U  #  -  55  X  105  C.G.S. 

Comparative  study  of  the  curves  of  figure  30. 

Curve  B  presents  a  slight  deviation  with  respect  to  curve  A. 
It  is  logical  to  hunt  for  an  explanation  for  this  different  shape  with 
the  aid  of  the  present  theory.   To  start  from  the  complete  theory, 
that  is  from  equation  (33)  would  involve  considerable  paper  work. 
Recourse  is  therefor^  had  to  the  simplified  theory  which  obviously 
produces  the  equation 

_^_   d3y  +  ^  dfy   a_j_1  +  _p_  dy  +  a_^  +  ep   =  Q       (   } 
Sav  dt3   Ta  di_2    TSv    Sav  dt    TSa    aT 

obtained  from  equation  (33)  by  disregarding  all  the  terms  having  R 

in  the  denominator.   The  same  method  that  produced  equation  (ll)  gives, 

therefore,  when  starting  from  equation  (35),  the  system 
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2T\(X2  +  co2)  +  vD(X.2  +  co2) : =  0 


T(3X2  -  co2)  +  2vD\  +  a(&  +  e)  *  Tp  =  0 
The  second  of  these  equations  gives 


o£  =  3X2   +  2  vm   a(a  +  Q   P  (  g) 

T       IT      I 


which,  suhstituted  in  the  first  equation,  gives 

v(a  +  6  +  Dp) 


b 


(2TX  +  vD) \hX2   +  2  ™  +  a(&+  °  +  P 

T       IT      I 


=  0 


that  is  after  transformations 

2D2I\v2  +  J8DTI  +  (aD  -  T)(a  +  e)]v  +  2TX.  \klTX2  +   a(a  +  e  )  +  TpJ  =  0 

(37) 

Compared  to  equation  (12)  this  equation  differs  only  by  the  presence 
of  coefficient  p   in  the  constant  term.   It  is  seen  that  for  a  real 
value  of  X     there  always  are  real  values  of  v,  provided  only  that  the 
determinant  is  positive,  that  is  that 

16DT2I(X2  +  Dp)  <  (a  +  e)(aD  -  T) 
The  maximum  value  Xm     is  then 

2  =  (a  +  0(aD-T)2  _  Dp  (38) 

i6dt^i 

So,    for  the   curve     A      (a  =  0.1+3   cm,      p   =  -3-5  X  105  C.G.S.    system) 

\a  =  21.1+7 
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and  for  the  curve  B   (a  =  1.73,  P  =  -55  X  105  C.G.S.  system) 

Xm  =  21.38 

Thus  we  obtain  two  practically  equal  values,  which  is  in  good  agreement 
with  the  shape  of  curves  A  and  B. 

We  shall  now  turn  to  the  calculation  of  the  corresponding  veloc- 
ities  vm.   Since  the  discriminant  of  equation  (37)  is  zero  in  this 
cage , 


2  =  -^-[UlTX.2  +  a  (a  +  e)  +  TpJ 


vm'     o 

Tri 

which  gives 

vm  =  255  cm/sec  =  8. K   ft/sec   for  curve  A 

vm  =  31^+  cm/sec  =  10.3  ft/sec   for  curve  B 

These  values  are  almost  identical  with  the  experimental  values  in 
figure  30.   But  above  all  this  calculation  has  proved  that  the  present 
theory  allows  for  the  deviation  of  curve  B  relative  to  curve  A, 
which  consists  essentially  of  a  displacement  of  the  maximum  and  the 
adjacent  parts  toward  the  higher  values  of  the  velocity. 

Equation  (36)  enables  the  values  of  the  angular  frequency  a) 
corresponding  to  the  maximum  X,m  to  be  computed  and  consequently  also 

the  values  of  the  frequency  n  =  — .   Thus 

n  =  10.5  cycles/sec   for  curve  A' 

n  =  12. 4-  cycles/sec   for  curve  B' 

These  values  very  closely  approximate  those  of  the  experimental  curves  A' 
and  B1   of  figure  30.   It  should  be  remembered  above  everything  else 
that  the  present  theory  justifies  a  higher  frequency  at  corresponding 
velocities  for  curve   B'  which  is  plainly  indicated  by  curves  A'   and 
B«. 

The  considerations  are  completed  by  the  calculation  of  the  limit 
of  the  frequency  N  when  the  velocity  increases  infinitely.   Concerning 
equation  (33),  it  is  readily  apparent  that,  when  v   increases  infinitely, 
the  majority  of  its  terms  approach  zero,  so  that  at  the  limit  this 
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equation  is  reduced  to 


Ie  dfy   /a 
Ta  dt2 


i  +  g     gp\ 

TSa   h  aT/* 


This  equation  represents  a  sinusoidal  function  of  the  frequency  go 
hence 


2   1/1   a    \   a  +  e   p 
I\S   D    /    ID     I 


Since  W  =  — 
2n 


N  =  17  cycles/sec   for  curve  A' 
N  =  18. 3  cycles /sec   for  curve  B' 


These  figures  are  in  perfect  agreement  with  the  experimental  curves   A' 

and  B » . 

In  conclusion  it  should  be  borne  in  mind  that  the  inclination  of 
the  spindle  axis,  that  is  the  angle  of  play,  has  only  a  small  effect 
on  the  stability  of  the  oscillations.   Therefore  this  angle  might  just 
as  well  be  determined  from  other  considerations  such  as  structural, 
geometrical,  or  operational  characteristics. 


B.  Effect  of  Viscous  Friction 

Viscous  friction  is  called  a  resistance  proportional  to  the  veloc- 
ity of  motion.   Now  suppose  that  the  castering  of  the  wheel  is  braked 

d© 
by  a  torque  equal  to  f  =  — ,  that  is 


dtJ 


f  dz 
a  dt 


The  same  line  of  reasoning  employed  in  the  foregoing  for  treating 
the  case  of  restoring  torque  will  produce  the  system 


2, 


_L  §~z  -. 


2   2 
v^  dt^ 


=  -RM 


(16) 
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z  -  y  =  TF  (17) 

l£*  =   _F(a  +  e)  -M  -  £  **                                          (18b) 

a  dt2             -  a  dt 


M  =  i  ^  +  ^  -  F6  (19) 

-   Sv  dt   Sa 

The  values  of  F  and  M  from  equations  (17)  and  (19)  entered  in 
equations  (l6)  and  (l8)  yield  the  two  equations 


1     d2y        1    dy        e  /  i        g\ 

— + +-y+ —  z  =  o 

Rv2  dt2       Sv  dt       T  VSa       Tj 


1    dy       a            i   d2z 
_  _  y  + + 

Sv  dt       t  a  dt2 


(20) 


£  *£  +  /I    +  aV   =  0  (21b) 

a  dt        \Sa       T/ 


The   values   of     z      from  equation    (20)   entered   in  equation    (21b)    give 
the   differential  equation 

d  y  .  /  I   .   f  \d3y  .  /ie     1   .   a   .   f  \d2y  j 


i     d  y  +  /  i    +     f  \d^y  +  /le_  +   1   +  a  +  f  \ 

Rav2  dtk  +   VSav  +  Rav2ydt3  +  \Ta  +  SRav2  +  TRy2  +  Siv/ 


dt2 


c 


a-L-i  +  £1)^1  +  ^jL±  y  =  o  (39) 

TSv    aT/dt    TSa 


It  is  readily  seen  that  the  second  remark  made  following  equa- 
tion (22)  is  applicable  here  too.   To  obtain  the  differential  equation 
of  motion  of  the  wheel,  simply  replace  y  by  z   in  the  preceding 
equation.   The  result  is  a  linear  differential  equation  of  the  4th  order 
of  the  form 

aQZ   +  a^z11-1-  +  agz   +  aoz1  +  a^z  =  0 
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In  connection  with  equation  (22),  it  was  seen  that  Routh's  condi- 
tions for  such  a  stable  oscillation  are 

a^a2  >  aga-^ 


2     ? 
ala2a3  >  a0a3  +  al  aU 


The  direct  study  of  these  equations  being  too  laborious,  will  be 
abandoned  in  favor  of  the  simplified  theory.   It  is  evident  that  this 
theory  produces  the  equation 

^_  ^1  +  (IL  +  -L)£l  +  (±J-±  +  £iW  +  *±JL  7  m  o    (Uo) 

Sav  dt3   \Ta   Sav/dt2   ^  TSv    aT/dt    TSa 


obtained  by  disregarding  all  terms  having  R  in  the  denominator  of 
equation  (39)- 

Incidentally,  this  is  the  same  equation  obtained  by  disregarding 
all  the  terms  of  higher  degree  in  l/v,  in  equation  (39) >  that  is  the 
terms  with  l/v2. 

The  motion  defined  by  the  simplified  theory  can,  therefore,  still 
be  regarded  as  the  asymptotic  limit  of  motion  defined  by  the  general 
equation  when  the  velocity  increases  indefinitely.   This  remark  applies 
also  to  equation  (35)  obtained  in  the  theory  of  the  restoring  torque 
from  equation  (33 ) ?  and  to  equation  (6)  with  respect  to  equation  (22). 

Equation  {ko)    is  of  the  form 

aoyI11  +  aiyi:[  +  a2yI  +  a3y  =  ° 


and  all  its  coefficients  are,  of  necessity,  positive.   The  unique  con- 
dition of  Routh's  stability  is  then 

ala2  ^  a0a3 
that  is,  by  taking  the  fundamental  relationship  Se  =  D  into  account, 
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\T    vj       \     v     a 


fD\  >  I(a  +  Q 

av 


which  can  also  be  written 

ID2fv2  +  JDTf2  +  I  (a  +  e)(aD  -  T)]v  +  Ta(a  +  e)f  >0      (*J-l) 

In  relation  to  the  velocity  v  the  left  hand  side  represents  a 
trinomial  of  the  2nd  order  admitting  of  two  real  roots  of  the  same  sign, 
V]_  and  V2,   or  else  imaginary.   Thus  the  condition  of  stability  is 
realized  in  the  following  two  cases: 

(1)  If  the  discriminant  A  of  the  trinomial  (4l)  is  negative 

(2)  If  the  discriminant  A  is  positive  and  if,  in  addition, 

v  <  v-<   or  V2  <  v 

This  argument  assumes  that  v  can  take  any  algebraic  value.   But  v 
is  necessarily  positive,  which  also  restricts  the  scope  of  the  con- 
clusions, as  will  be  seen.   The  discriminant  A  is  written 

A  =  D2T2f^  -  2IDT(a  +  e)(aD  +  T)f2  +  I2(a  +  e)(aD  -  T)2 

With  respect  to  f2,  A  is  a  trinomial  of  the  second  degree  having 
two  positive  roots 


1  ■  iP^(fc  -  tf«) 

If    dt 

f2  =  ^szn^  +  vs)  (u2) 


DT 

Before  proceeding  to  the  discussion,  three  preliminary  remarks  are  noted: 

(l)  If  aD  >  T,  the  left  side  of  the  inequality  (kl)    is  a  sum  of 
positive  terms  and  the  stability  condition  is  confirmed. 
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(2)    If     DTf2  +   I  (a  +  e  )(aD  -  T)  >   0,    that   is    if     f  >   cp     by 
substituting 


V    TYP 


DT 

the  stability  condition  is  proved  again  for  the  same  reason. 

(3)  When  aD  <  T,  then   f  1  <   cp  <  f2. 

The   second  of  this   double    inequality  is   evident.      As   to  the   first, 
it   is   sufficient  to  show  that 

aD   +  T   -   2  VaDT  <  T   -  aD 

0 

that   is 

2aD  <  2l/aDT 

or  else 

aD  <  T 

which  is  true  by  assumption. 

So,  if  aD  <  T  and  f]_  <  f  the  oscillation  is  always  stable 
whatever  the  velocity  may  be,  because 

(a)  either  f]_  <  f  <  fg  and  A  <  0 

(b)  or  f  >  cp 

When  f  <  f^,  then  A  >  0  and  the  oscillation  is  stable  only  when 
v  <  v-j_,  or  Vg  <  v. 

To  sum  up:   stability  exists  no  matter  what  the  velocity  may  be 

(1)  if  aD  >  T 

(2)  if  aD  <  T  and   f  >  f  1 

But,  when  aD  <  T  and  f  <  f-^,  the  oscillations  are  unstable  for  any 
velocities  between  v]_  and  vp. 
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According  to  these  conclusions  when  a  tire  shimmies,  the  phenomenon 
is,  as  a  rule,  observed  only  in  a  certain  speed  range  (v-|_,  vp),  which 
can  become  quite  extended. 

To  overcome  shimmy,  if  aD   is  not  >T,  it  is  sufficient  to  use  a 
hydraulic  damper  having  a  coefficient  of  viscous  friction   f  greater 
than  the  coefficient  of  friction   f]_  defined  by  equation  (^2).   This 
solution  has  been  frequently  used,  especially  by  the  Americans. 
Unfortunately,  this  solution  is  not  without  drawbacks,  because  it  is 
not  simple,  it  increases  the  cost,  requires  a  certain  amount  of  main- 
tenance, and  increases  the  weight  of  the  vehicle,  which  is  particularly 
annoying  when  an  airplane  is  involved.   The  best  solution  is  very 
likely  the  simple  realization  of  the  relation  aD  >  T. 

The  coefficient  of  minimum  viscous  friction  to  overcome  shimmy  is 
therefore   f]_.   For  a  value  of  f  slightly  below  f]_,  shimmy  still 
exists  for  velocities  near  to  vq  given  by  the  equation 

v02  =  Ta(a  +  °  (1*3) 

ID2 


obtained  by  substituting     f  =   f]_      in  equation   (^l).      This   velocity     vfi 

can  therefore  be  regarded  as   the  most  dangerous   velocity  for  shimmy. 
With  respect   to  the  velocity     vm     defined  by  equation   (12a),    that   is 
relative  to  the  maximum   inversion  velocity  in  the   absence   of  the   damper, 
we   get  therefore 


2\/aDT 


(MO 


'm       aD  +  T 


It  follows,  as  is  readily  seen,  that,  as  a  rule 

v0^  vm 


C.  Effect  of  Solid  Friction 

Solid  friction  or,  as  it  is  also  called,  constant  friction,  is  a 
constant  resistance  independent  of  the  speed  of  motion.   Suppose  that 
the  castering  of  the  wheel  is  braked  by  a  constant  torque  C.   The  case 
of  solid  friction  is  more  difficult  to  treat  analytically  than  that  of 
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viscous  friction,  but  if  it  is  only  a  question  of  defining  the  torque 
of  solid  friction  C  sufficient  to  prevent  shimmy,  the  following 
expedient  can  he  resorted  to:   the  solid  or  the  viscous  friction  damper, 
both  act  through  the  dissipation  of  energy  transmitted  to  the  wheel  by 
the  reactions  of  the  ground.   It  therefore  seems  logical  to  assume  that 
the  solid  friction  damper  will  be  sufficient  to  eliminate  shimmy  when 
it  absorbs  at  each  oscillation  the  same  amount  of  energy  as  the  viscous 
friction  damper,  that  is  with  f]_   (formula  (U2))  for  coefficient  of 
viscous  friction. 

The  energy  dissipated  during  a  quarter  cycle  of  amplitude  0m  by 
solid  friction  damper  is  W  =  Cdm.  On  the  other  hand,  the  energy  dis- 
sipated by  a  viscous  friction  damper  under  the  same  conditions  is 

f0m   d0 

w  =  /   f  —  <ae 

Jn  dt 


When  this  damper  is  exactly  able  to  make  this  oscillation  represent 

the  boundary  between  stability  and  instability,  a  sinusoidal  oscillation 

9  =  9m   sin  cut  and  f  =  f i 

is  produced. 

The  preceding  integral  gives 

W  =  J  f^2 
Equating  the  two  values  of  W  produces 

C  =  J  tjcxtin  (k$) 


In  this  formula  all  parameters,  save  the  angular  frequency  <d,  are 
known.   It  was  shown  that  cd  is  an  increasing  function  with  the  veloc- 
ity but  varying  little  with  this  velocity.   So,  when  v  increases 
indefinitely  the  limiting  value  of  go  should  be  taken.   This  value  is, 
as  seen. 


{* 


a  +   e 
oj  =  J 

ID 
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and  formula  (h^>)   becomes  thus 


it  a  +  e    >/T  -  VaD 
k       D  /f 


c  =  *  i±i  yi  -JaiJ  em  (ii6) 


Hence  the  expression  of  the  solid  friction  torque  is  sufficient  and 
even  a  little  more  than  sufficient  to  stifle  shimmy.   This  torque  is 
proportional  to  the  maximum  amplitude  0m  while  the  viscous  friction 
coefficient   f]_  necessary  to  achieve  the  same  purpose  is  independent 
of  0m.   This  remark  is  extremely  important,  because  it  explains  the 
fact  pointed  out  frequently  that,  when  solid  friction  is  involved, 
shimmy  is  not  produced  in  general  in  the  absence  of  an  initial  dis- 
placement of  sufficient  amplitude.   In  other  words,  when  solid  friction 
is  involved  -  and  there  always  is  more  or  less  friction  in  the  spindle 
axis,  even  in  the  absence  of  a  special  damper  -  shimmy,  in  order  to  be 
produced,  has  to  be  initiated  or  stimulated.        ( 

To  avoid  such  initiation  at  landing,  it  might  be  useful  to  employ 
a  wheel  castering  mechanism  to  keep  the  angle  small  at  the  moment  of 
contact  with  the  ground.   So  in  formula  (k6)    for  the  practical  calcu- 
lation of  torque  C,  one  should  use  for  %  the  maximum  amplitude 
likely  to  be  accidently  produced.   The  determination  of  this  angle  9m 
raises  some  difficulties.   It  has  already  been  proposed  to  take  for  9m 
the  angle  at  which  the  skidding  of  the  tire  on  the  ground  occurs.   This 
is  justified  if  one  assumes  that  for  the  greater  angles  sufficient 
braking  action  is  assured  by  the  skidding  of  the  tire  on  the  ground. 
Under  these  conditions  formula  (k6)   can  be  given  a  different  form. 

Consider  a  sinusoidal  oscillation 

9  =   0m  sin  out 

and  assume,  as  in  the  approximate  theory,  that  the  oscillation  of  the 
wheel  is  largely  due  to  the  lateral  force  F  exerted  at  a  distance   e 
behind  the  point  of  contact.   In  this  case 

z  =  a0m  sin  cut  — -  =  -aSnjO^sin  cot 

dt2 

o 
-  £-£■  =   -F(a  +   e)      hence       F(a  +  e)   =   ISmO^sin  cut 
a  dt2 
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and  consequently 


3     =  F 
7m    m 


a  +  e 


laf 


Fm  being  the  lateral  force  which  produces  the  skidding.   On  assuming 
that  is?      is  still  approximately  equal  to 

a  +  e 


ID 


formula    (h-6)   becomes 


it  (a  +    g)    |/t  -   \ZaD 


'm 


(V7) 


The  limitation  of  the  amplitude  of  shimmy  as  a  result  of  the  skidding 
of  the  tire  when  the  lateral  force  F  reaches  the  limit  of  adhesion  Fm, 
explains  in  particular  why  the  terrain  most  favorable  to  intense  shimmy 
is  macadam  and  dry  and  rough  concrete,  and  why  on  the  other  hand,  gravel 
and  light  soil  have  a  certain  damping  effect  on  the  oscillations  of 
shimmy  and  thus  limit  the  amplitude.   It  also  explains  why  the  load  is 
a  factor  having  considerable  effect  on  the  intensity  of  shimmy.   All 
these  facts  have  been  checked  many  times  on  automobile  wheels  as  will 
be  shown  later. 

Many  designers  have  adopted  the  solid  friction  damper  as  a  means 
to  stop  shimmy.   This  solution  has  the  advantage  of  being  simpler  and 
cheaper  than  the  hydraulic  damper,  but  it  has  the  drawback  of  "stiffening" 
the  castering  of  the  wheel  and  so  making  it  harder  to  maneuver  the  air- 
plane on  the  ground.   This  drawback  makes  this  solution  to  the  problem 
extremely  annoying  in  the  case  of  airplanes  of  large  weight,  and  prac- 
tically limits  its  use  to  light  aircraft.   Furthermore,  it  is  very 
difficult  to  ensure  a  well  defined  and  perfectly  constant  value  for 
solid  friction.   There  always  is  the  risk  of  too  much  or  not  enough 
tightness.   All  these  drawbacks  show  that  solid  friction  is  not  the 
ideal  remedy  against  shimmy. 


Tests  made  during  the  war  in  Germany  by  the  Motor  Institute,  a 
branch  of  the  Institute  of  Technology  at  Stuttgart,  and  by  the 
Consolidated  Vultee  Aircraft  Corporation  in  the  U.  S.  A.  (See  Aero 
Digest  of  March  15,  19^,  pp.  13^  to  137;  Conquest  of  Nose  Wheel  Shimmy 
by  C.  B.  Livers  and  J.  B.  Hurt)  have  shown  that  nose  wheel  shimmy  of 
the  tricycle  landing  gear  could  be  eliminated  by  the  use  of  two  wheels 
side  by  side  on  one  axle.   This  remedy  is  similar  to  the  preceding  one, 
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because  the  two  wheels  are  clamped  on  the  same  axle  with  the  possibility 
of  differential  rotation.   The  castering  of  the  unit  is,  of  necessity, 
accompanied  by  a  certain  skidding  of  the  wheels  on  the  ground,  skidding 
equivalent  to  a  solid  friction.   This  solution  presents  the  same  draw- 
backs as  the  foregoing,  it  makes  handling  on  the  ground  difficult, 
especially  when  a  heavy  airplane  is  involved. 

Considerations  of  practical  hydraulic  dampers. 

In  the  aforementioned  study  of  viscous  friction  it  had  been 
assumed  that  the  resistance  was  rigorously  proportional  to  the  rate 
of  motion.   In  practice  this  law  of  damping  is  never  realized  for 
actual  hydraulic  dampers,  because  the  true  law  of  resistance  in  terms 
of  velocity  is  the  resultant  of  three  effects: 

(1)  An  effect  of  the  viscosity  of  the  liquid  corresponding  some- 
what to  a  resistance  proportional  to  the  velocity. 

(2)  A  kinetic  effect  or  so-called  hydraulic  friction  or  turbulence, 
corresponding  to  a  resistance  proportional  to  the  square  of  the  velocity. 
This  is  a  resistance  of  purely  kinetic  origin  that  is  produced  in  every 
device  forcing  a  liquid  through  a  narrow  orifice,  even  if  a  liquid 
without  appreciable  viscosity  is  involved. 

(3)  A  valve  effect,  the  preceding  orifices  being  often  duplicates 
of  spring  valves  that  remain  closed  as  long  as  the  pressure  remains 
below  a  certain  value,  depending  on  the  force  of  the  spring.   On  opening 
they  increase  the  area  of  the  passage  and  prevent  the  torque  from 
increasing  with  the  square  of  the  speed.   All  sorts  of  dampers  with 
varied  characteristics  are  obtained  by  appropriate  control  of  the  two 
compensating  effects. 

As  a  rule,  a  damper  is  characterized  by  the  curve  representing 
the  resistance  it  offers  to  motion  in  terms  of  the  rate  of  the  motion. 
The  ordinary  dampers  have  generally  characteristic  curves  similar  to 
that  shown  in  figure  36-   The  dashed  lines  represent  the  characteristic 
curve  of  a  solid  friction  damper  and  the  theoretical  curves  of  pure 
viscous  and  pure  hydraulic  friction.   The  viscosity  of  a  liquid  decreases 
when  the  temperature  increases,  so  that  the  typical  curve  of  an  actual 
damper  is  always  deformed  through  diminution  of  all  its  ordinates. 
With  such  a  graphical  representation  it  can  be  immediately  determined 
whether  or  not  shimmy  disappears  if  one  assumes  that  two  dampers  are 
equivalent  from  the  point  of  view  of  oscillation  stability  when  they 
absorb  the  same  amount  of  energy  during  one  period.   In  that  case,  the 
viscous  damping  indicated  by  the  theory  (formula(^2) )  is  simply  trans- 
ferred on  to  the  characteristic  diagram  of  the  damper  which  gives  a 
straight  line  passing  through  the  origin.   Considering  the  shape  of 
the  curve  given  for  an  actual  damper  it  is  seen  that  shimmy  is  likely 
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to  occur  at  low  angular  velocities  and  hence  at  low  amplitudes.   If 
the  curve  of  the. viscous  damping  required  is  tangent  to  the  damper 
characteristic  at  zero  velocity,  the  damper  is  exactly  strong  enough 
to  prevent  divergent  oscillations.   If  the  slope  is  less  than  that  of 
the  damper  characteristic  at  zero  velocity,  it  is  more'  than  sufficient 
to  eliminate  shimmy.   Conversely,  if  the  slope  is  greater,  the  torque 
at  low  velocities  should  be  greater  than  that  that  can  he  given  hy  the 
damper,  and  shimmy  results.   Lastly,  if  curve  of  the  damping  necessary 
intersects  the  curve  of  the  damper  at  a  comparatively  high  velocity 
above  which  the  damper  gives  a  more  than  ample  torque,  shimmy  amplitude 
is  limited  to  that  at  which  the  energy  received  by  the  wheel  at  each 
cycle  is  equal  to  the  energy  dissipated  by  the  damper.   If  the  curve  of 
the  resistance  is  definitely  above  the  characteristic  curve  of  the 
damper,  the  latter  is  powerless  to  limit  the  amplitude  of  shimmy.   As 
the  smallest  amplitudes  possible  are  already  negligible,  the  only  solu- 
tion is  to  stifle  shimmy  completely.   On  the  other  hand,  the  theory  of 
shimmy  being  applicable  only  to  low  amplitudes,  if  the  damping  provided 
is  insufficient  for  stifling  shimmy  completely,  the  amplitude  attained 
may  be  high  enough  so  that  the  equations  are  no  longer  applicable,  and 
the  motion  continues  then  to  increase  in  intensity  instead  of  remaining 
in  a  certain  state  of  equilibrium.   This  is  the  reason  why,  in  certain 
cases,  it  was  deemed  preferable  to  adopt  a  greater  damping  than  strictly 
necessary  in  order  to  assure  convergence;  any  oscillatory  effects  due 
to  deviations  from  linearity  that  could  occur  are  thus  damped  out. 

Equation  (h-'j)   presents  an  unusual  feature.   Taking  the  coefficients 
which  we  computed  for  the  Schlippe-Dietrich  tire,  that  is 

6  =  h   cm      T  =  157  X  lO"10  (C.G.S.       D  =  15-7  X  10-10  (C.G.S. 

system)  system) 

and  computing  by  this  equation  the  variations  of  the  friction  coeffi- 
cient C  necessary  to  stifle  shimmy  for  various  caster  lengths  a, 
the  following  results  are  obtained: 


0  cm  C  =  3.  Ill  Fm  C.G.S.  system  a' =  5  C  =  2.09  Fm  C.G.S.  system 
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The  data  are  represented  graphically  in  figure  37-   Although  this  curve 

T 
decreases  consistently  when  a  varies  from  0  to  —  =  10,  it  is  seen 
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that  C  still  is  comparatively  great  for  a  values  of  the  order  of 
magnitude  of  6  cm,  the  decrease  of  C   being  especially  rapid  when  a 
approaches  T/D.   But  since  C  measures  to  some  extent,  the  eventual 
intensity  or  the  danger  of  shimmy,  it  follows  that  the  usual  caster 
lengths  are  displaced  precisely  in  the  zone  particularly  dangerous  for 
shimmy.   The  matter  becomes  plainer  if  one  notes  that  on  certain  tires 
the  curve  C   passes  through  a  maximum  for  a  value  of  a  varying 
between  0  and  T/D.   In  fact  if  it  is  assumed  that  e  =  2   instead 
of  e  =  h     for  the  preceding  tire,  a  figure  still  likely  for  a  tire  of 
that  size,  equation  (1+7)  then  permits  the  calculation  of  the  following 
values: 

C.G.S.  system  a  =  5  cm  C  =  1.63  Fm  C.G.S.  system 


0  cm     C 

=  1.57  Fm 

0.5 

1-52 

1 

1.61 

2 

l.fk 

3 

1.77 

k 

1.7^ 

5   cm 

C     : 

=  1-63  Fm 

6 

1.1* 

7 

1.18 

8 

.865 

9 

.VT5 

10 

0 

The  results  are  represented  by  the  curve  of  figure  38.   The  curve  C 
passes  through  a  maximum  near  a  =  3.   This  remark  can  be  illustrated 
by  an  analytical  calculation:   computing  the  derivative  of  C   from 
equation  (h'j)    gives 


!^(3»'-#^) 


dC 
3a" 


when  a   is  infinitely  small  the  foregoing  quantity  assumes  a  negative 

infinite  value  which  makes  it  possible  to  define  the  shape  of  the 

T 
curve  C   for  very  small  values  of  a.   When  —  >  3? ,   the  bracketed 

term  (a  trinomial  of  the  second  degree  with  respect  to  \fa.)   cancels 
out  for  the  two  roots 


/a  =  — 


3e 


It  is  readily  apparent  that  the  smallest  of  these  values  of  a 
corresponds  to  a  minimum  of  C  and  the  largest  to  a  maximum.   In  the 
present  case  this  calculation  gives 
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a^_  =  0. 15  cm 


l2 


=  2.96 


cm 


_?min 


rmax 


1.14-8  Fm  C.G.S.  system 
1.77  Fm  C.G.S.  system 


The  condition  necessary  and  sufficient  for  curve  C   to  assume  a 
maximum  is 

This  condition  is  realized  on  certain  tires.   The  usual  caster  lengths 
are  then  precisely  the  most  dangerous  for  shimmy.   This  remark  stresses 
the  importance  of  increasing  the  caster  length  to  combat  shimmy  at  its 
very  source.   Equation  (12)  makes  it  possible  to  establish  that,  when 
the  velocity  v  varies,  the  maximum  divergence   X,m   is  given  by  the 
equation 

lSm2!^2  =    (a  +  e)(aD  -  T)2 

m 

But  a  calculation  readily  proves  that  when  =-  >  2e  (a  condition  realizec 
for  a  large  number  of  tires)  and  a  is  made  to  vary,  this  divergence  X,m 
passes  through  a  maximum 


A 


T  +  De  T  +  De 


6DT 


31 


for  a  caster  length 


T  -  2De 
3D 


The  coefficient  C   seems  better  than  the  divergence  X  for 
appraising  the  danger  of  shimmy. 

Doctor  Langguth's  Experiments  (19*4-1). 

It  is  timely  to  compare  those  theoretical  conclusions  with  the 
experimental  results  obtained  in  19*+1  by  Doctor  Langguth  and  reported 
by  M.  P.  Mercier.   He  dealt  with  experiments  on  260  X  85  tires,  mounted 
on  a  rolling  mat,  at  a  speed  of  50  km/h.   The  author  recorded  the 
maximum  angular  swivel  or  deflections  of  the  wheel  in  terms  of  caster 


lengths  for  castering  angles  of  -3  ,  0  ,  and  +3  .   Figure  39  summarizes 
the  results.   The  four  essential  conclusions  of  this  work  are  the 
following: 
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(1)  The  caster  angle  is  of  little  importance 

(2)  The  angular  deflections  pass  through  a  maximum  at  a  =  7  cm 
approximate ly 

(3)  They  cancel  out  at  around  a  =  12  cm 

(h)    Increasing  the  load  increases  the  angular  deflections  prac- 
tically in  the  same  proportions 

Before  interpreting  these  results,  it  is  advisable  to  understand 
the  meaning  of  angular  deflection  (swivel)  studied  by  the  German  author, 
since  this  quantity  does  not  figure  directly  in  the  present  theory. 
Nevertheless  the  theory,  considered  superficially,  seems  to  indicate 
that,  every  time  there  is  an  oscillatory  divergence,  the  angular  deflec- 
tion should  increase  indefinitely.   This  apparent  contradiction  is  due 
to  the  fact  that,  in  order  to  treat  the  problem,  the  data  must  be 
linearized,  which  is,  perfectly  legitimate  so  long  as  the  angular 
deflections  attain  no  unduly  great  values.   Furthermore,  the  present 
theory  assumes  the  adhesion  of  the  tire  to  the  ground  to  be  infinite, 
while  in  reality  this  adherence  is  limited.   The  limit  of  adherence 
explains  in  particular  why  the  maximum  angular  deflections  were  approxi- 
mately proportional  to  the  load  on  the  tire. 

Undoubtedly  the  present  theory  does  not  permit  the  ordinates  of 
the  curves  of  Langguth  to  be  tied  quantitatively  to  the  characteristic 
parameters  of  the  tire,  but  it  does  seem  justified  to  assume  that, 
everything  else  being  equal,  the  maximum  angle  must  vary  as  C,  and  the 
existence  of  the  maximum  discovered  by  Langguth  seems  to  be  in  good 
agreement  with  the  foregoing  theoretical  conclusions.   But  from  the 
present  writer's  point  of  view  the  major  importance  of  these  experiments 
is  that  they  establish  the  existence  of  a  caster  length  of  about  12  cm 
below  which  (for  a  <  12)  the  oscillations  are  unstable  and  above  which 
(for  a  >  12)  the  oscillations  remain  stable.   In  other  words,  Langguth *s 
experiments  confirm  quite  well  the  essential  conclusion  of  the  present 
theory,  that  is  the  stability  of  the  oscillations  with  a  sufficiently 
large  caster  length.   And  this  conclusion  is  quantitative  as  well  as 
qualitative:   Langguth *s  260  X  85  tire  was  of  the  same  size  as  that 
used  by  Schlippe  and  Dietrich.   Therefore,  it  is  legitimate  to  assume 
that  these  tires  have  practically  the  same  characteristics  T  and  D. 

m 

But,  it  has  been  shown  that  the  caster  length  a  =  —  of  the  Schlippe  - 

Dietrich  tire,  necessary  for  stable  oscillations  was  about  10  cm.   The 
quasi-agreement  of  this  figure  with  Langguth's  experimental  value  is 
extremely  remarkable,  in  as  much  as  the  slight  residual  discrepancy 
between  the  two  values  can  be  readily  attributed  to  a  difference  in 
tire  inflation  in  the  two  test  series. 
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VI.  MATHEMATICAL  THEORY  OF  COMPLEX  SHIMMY 

The  term  complex  shimmy  describes  shimmy  with  two  degrees  of 
freedom,  the  second  degree  of  freedom  resulting  from  a'  certain  lateral 
elasticity  of  the  pivoting  axis. 

As  in  the  study  of  elementary  shimmy,  we  shall  start  with  the 
simplified  theory,  disregarding  the  torsional  moment  of  the  bottom 
part  of  the  tire,  that  is  assume  that  the  drift  directly  accompanies 
the  lateral  force  F  which,  in  more  exact  terms,  amounts  to  assuming 
the  coefficient  of  turn  R  of  the  tire  to  be  infinite. 

To  put  the  problem  into  mathematical  form  we  return  to  figure  8 
and  to  the  line  of  reasoning  that  produced  the  system  of  equations  (l), 
(2),  and  (3).   But,  first,  in  order  to  account  for  the  lateral  elasticity 
of  the  pivoting  axis,  it  is  advisable  to  replace  figure  8  by  figure  ^0: 

With   Z  as  the  lateral  displacement  of  the  pivot  under  the  action 
of  the  lateral  force  Fq, 

I   =  TqFo 

where  Tq  is  called  the  coefficient  of  the  lateral  elasticity  of  the 
pivot. 

Equation  (l)  always  gives  the  relation 

5   =   9  +  i|r    =  DF 
but  now 

e  =  z  -  z 

a 
so  that  equation    (l)   must  be   replaced  by 


1  djr  +  £  _   l_  =  DF 

v  dt        a       a 


Equation  (2)  is  obviously  not  modified. 
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As  to  equation  (3)>  it  is  transformed  as  follows:   If  Iq   is  the 
inertia  of  the  wheel  about  the  vertical  axis  passing  through  the  center 
of  gravity  of  the  wheel,  then  I  =  Iq  +  ma^ f   where  m  =  mass  of  wheel. 
The  angular  acceleration  of  the  wheel  about  this  axis  ,1s  a  function  of 
the  sum  of  the  moments  of  the  forces  F  and  Fq,  hence 


&26 
d^ 


Fna  -  Fe 


now 


Z 
and  Fn  =  =— 
0   T0 


Therefore 


The  introduction  of  a  new  variable   Z   calls  for  a  fourth  equation, 
which  is  obtained  by  noting  that  the  acceleration  of  the  center  of 
gravity  0   is  due  to  the  sum  of  the  forces   F  and  Fq.   Thus 

m  —  =  -F  -  F0  =  -F  -  =L 
dt2       °      To 


Hence  the  following  system 


I  &  +  £  _  I  =  DF  m 

v  dt   a   a 


z  -  y  =  TF  (k9) 


(50) 
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m 


d2z 
dt2 


(5D 


To  obtain  the  differential  equation  in  z  defining  the  motion  of 
the  wheel,  simply  eliminate  F,  Z,  and  y  from  these  four  equations, 
as  follows: 

The  four  equations  are  linear  and  homogeneous.   So,  if 

y  =  Aeat 

is  a  particular  solution  of  the  system,  three  particular  solutions 
for  z,  Z,  and  F  of  the  form 


z  =  Be 


st 


I   =   Ce 


st 


F  =  De 


st 


are  obtained.   Thus  the  equation  in  s  of  the  foregoing  system  or  the 
characteristic  equation  of  the  system  can  be  obtained  directly  by 
replacing  the  unknowns  y,  z,  Z,  and  F  by  the  values  above  and  then 
eliminating  A,  B,  C,  and  D  from  the  four  equations  of  the  system  by 
Cramer's  method  of  determinants.   Hence  we  obtain  an  equation  in  s 
in  the  form 


(y) 


(z) 


(i) 


(F) 


1 

—  s 
v 


1 
a 


10  s2 


ms 


1 
a 


l0  g2 
a  s 


a 
Tn 


-D 
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The  development  of  the  determinant  gives  the  characteristic  equation 


mlnE 


E*o  s5  +  ^  sk  +   i  &  +  Zo 


In  +  ma^ 


T 


v  V  T 


s3  +  I 

T 


D(l0  +  ma?] 


me 


S2  + 


a(a  +  e )  s  +  a  +  e  _  q 

n-0 


TT0v 


The  differential  equation  sought  is  thus 


mI0TT0  d5z    _  nm   d^z   IOTO  +  IT  d3z 

+  itiIqDTq  + +  (ID  +  meTo) 4 


,d2z 


dt5 


dt 


dt 


3 


dt' 


a(a  +  € )  dz   /     v    _ 
— - —  +  (a  +  e)z  =  0 


dt 


(52) 


It  should  be  noted  that,  when  Tq  =  0,  that  is  when  the  pivot  has  no 
lateral  elasticity,  the  foregoing  equation  reverts  to  equation  (6). 
Two  particular  cases  with  important  conclusions  will  be  analyzed. 

1.  Wheel  Absolutely  Rigid 


Assume  that  simultaneously 

T  =  0  D  =  0 


0 


that  is  a  wheel  without  any  lateral  elasticity,  hence  without  drift. 
Equation  (52)  is  simplified  and  becomes 


!oT0  d3z   a(a  +  e )  dz    ,      N 

-  + +  (a  +  e  )z  =  0 


dt- 


dt 


It  is  readily  seen  that  the  condition  of  Routh's  convergence  can 
never  be  realized  with  such  an  equation.   In  such  a  case,  if  Tq  is 
not  zero,  that  is  if  the  pivot  has  a  certain  lateral  elasticity,  there 
always  will  be  a  divergent  oscillation,  that  is  shimmy. 
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The  lateral  elasticity  of  the  pivot  can  he  the  sole  generator  of 
shimmy  without  it  heing  necessary  for  the  mechanical  properties  of  the 
tire  themselves  being  involved.   For  certain  particular  combinations 
of  tire  factors  and  speed,  the  lateral  elasticity  of  the  pivot  can 
therefore  become  the  factor  facilitating  shimmy. 

2.  Pivot  With  Great  Lateral  Freedom 

Centering  attention  on  parameter  Tq  it  is  assumed  that  To  =  °°, 
that  is  that  the  pivot  is  absolutely  free  to  shift  laterally  under  the 
action  of  a  force  no  matter  how  small.   Equation  (52)  then  becomes 

raI0T  d5z    T  _  d'+z   Io  d^z      d2z   _ 

v  dt5      dt*   v  dc3    dt2 


d2z 
Putting  =  u,  the  preceding  equation,  which  is  independent 

dt2 

of  a  is  transformed  into  a  differential  equation  of  the  3rd  order 

for  which  Routh's  stability  condition  becomes,  after  simplifications, 

I0  >  m  |  e 

It  is  readily  apparent  that  this  condition  is  realized  generally 
with  ordinary  tires.   The  ratio  T/D  is  a  length  generally  very  close 
to  the  radius  r  of  the  tire;  and  considering  that 

In  =  -  mr2 
u   K 

K  is  a  coefficient  generally  near  2.   The  preceding  inequality  thus 
becomes  r  >  2e,  which  is  amply  confirmed  for  all  customary  tires. 

Since  the  pivot  is  completely  free  to  shift  laterally,  the  oscil- 
lations are,  naturally,  convergent,  that  is  there  can  be  no  shimmy. 
Hence  the  possibility  of  a  new  remedy  for  removing  shimmy  by  giving 
the  pivot  complete  lateral  freedom. 

Based  upon  more  elementary  reasoning  of  qualitative  rather  than 
quantitative  nature,  various  American  technicians  had  been  led  to 
assume  that  the  possibility  of  a  wheel  to  shift  laterally  without 
pivoting  should  be  capable  of  suppressing  shimmy.   To  realize  this  aim, 
two  types  of  mechanical  arrangements  were  recommended  and  tested,  either 
on  the  airplane  or  on  the  belt-machine  apparatus. 
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The  first  method  consists  in  letting  the  free  wheel  slide  on  its 
axis,  by  allowing  the  necessary  play  at  the  end  of  the  hub.   The  axis 
was  slightly  curved  so  that  the  wheel  had  a  tendency  to  remain  in  the  ' 
center.   This  method  was  tested  by  the  NACA  on  the  Weicjc  W-l  airplane 
and  revealed  itself  as  being  effective  as  long  as  the  pivoting  remained 
below  13°  in  one  direction  or  the  other.   But  other  tests  have  given 
no  satisfaction.   In  the  light  of  the  present  theory  it  is  easy  to 
understand  the  reason  for  this  failure.   Whatever  the  perfection  of 
lubrication,  the  friction  of  the  wheel  on  its  axle,  in  the  motion  of 
the  lateral  displacement,  cannot  be  considered  as  negligible.   This 
friction  reestablishes,  in  a  certain  measure,  a  lateral  restraint.   It 
follows  that,  if  the  stability  condition 

J0  >  m  5  € 

is  not  greatly  exceeded,  the  phenomenon  of  shimmy  may  reappear. 

The  second  method  consists  in  utilizing  two  pivots  coupled  as 
illustrated  in  figure  kl.      The  NACA  tested  this  system  in  the  laboratory 
and  found  that  it  eliminated  shimmy.   It  was  also  studied  by  E.  R.  Warner 
(ref.  h) .      The  operation  is  obvious.   This  arrangement  realizes  complete 
lateral  freedom  for  the  wheel  pivot.   It  seems  quite  superior  to  the 
first,  although  structural  complications  are  involved  and  landing  gear 
retraction  is  made  more  difficult.   The  American  authors,  E.  S.  Jenkins 
and  A.  F.  Donovan  (ref.  5)  believe,  however,  that  it  could  be  made  as 
light  as  a  damped  system  and  that  the  elimination  of  the  detrimental 
effects  of  damping  friction  on  ground  handling,  as  well  as  the  elimina- 
tion of  the  upkeep  of  a  damper,  would  have  appreciable  advantages.   The 
additional  complications  introduced  by  this  system  of  double  pivot  are 
compensated  for  in  part  by  dispensing  with  the  scissors  which  transmit 
the  shimmy  torque  to  the  damper  in  the  ordinary  arrangements. 

Clearly,  all  the  advantages  presented  by  this  arrangement  are 
found  in  that  advocated  in  the  present  report  (simple  lengthening  of 
caster  length)  and  which  has  the  added  advantage  of  greater  simplicity. 

The  two  foregoing  examples  prove  that,  in  the  fight  against  shimmy, 
the  existence  of  a  certain  amount  of  lateral  elasticity  may  sometimes 
turn  out  beneficial,  sometimes  harmful,  depending  on  the  magnitude  of 
this  elasticity  and  the  values  of  the  characteristic  coefficients  and 
the  velocity.   A  clarification  of  this  point  calls  for  a  complete  study 
of  the  stability  condition  of  equation  (52).   Unfortunately,  too  much 
paper  work  is  involved  to  be  undertaken  by  us.   Furthermore  there 
does  not  seem  to  be  much  practical  interest  in  such  a  discussion;  in 
fact,  it  has  been  shown  that  it  is  comparatively  easy  to  stop  shimmy 
i.n  the  case  of  a  rigid  pivot;  and  the  realization  of  a  lateral  elasticity, 


92  NACA  TM  1337 


properly  speaking,  in  the  pivot  will  present  difficulties  of  construc- 
tion and  drawbacks  of  practical  nature. 
t 

Effect  of  a  restoring  torque 

Assume  a  restoring  torque  of  the  form  Cr  =  p0.   Obviously,  equa- 
tions (k8) ,    (^9)>  and  (51)  remain  unaltered.   Equation  (50)  is  replaced 

by 

Wd2z   d2z\  =  al  _ 


a\dt2   dt2/   T0        a 

the  same  way  as  equation  (l8)  was  replaced  by  equation  (l8a)  to  obtain 
equation  (33). 

The  result  is  a  system  of  four  equations  which  is  resolved  by 
means  of  the  method  used  for  equations  (h8) ,  (U-9) ,  (50),  and  (51). 
Consequently 

mI0TTo  d5z    T     d^z   IT  +  JoTo  +  mTToP  d3z 

"  U  ^4  +  mloDTo  ^f  +  — +  (ID  +  meT0  + 

V    dtP         dt^  v        dt3 


d2z   a(a  +  e)  +  (T  +  T0)p  dz    ,        _  . 
mDToP) + +  (a  +  e  +  Dp)z  =  0 

dt2  v  dt 


If  p  =  0,  equation  (52)  is  obtained  again. 


Complete  Theory  of  Complex  Shimmy 

If  the  torsional  elasticity  and  the  turning  effect  of  the  tire 
are  involved,  the  foregoing  line  of  reasoning  results  in  a  system  of 
five  equations. 


1  d^y 
v2  dt2 


2. 

RM 


z  -  y  =  TF 
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I0/d2z   d2Z\    _    aZ 


a 


_Fe  +  —  -  M 
,dt2   dt2/         T0 


1  dy     z  -  l 

M  =  —  — -  +  — Fe 

—   sv  dt    Sa 


m 


d2_z 
dt2 


{'♦*) 


The  differential  equation  sought  is  obtained  by  eliminating  F,  M, 
Z,  and  y  from  the  five  equations.   This  system  can  also  be  studied  and 
resolved  more  rapidly  by  the  method  employed  for  equations  (U8),  (h^) , 
(50),  and  (51).   Even  then  the  calculations  are  long  and  therefore  not 
attempted. 


VII.  SHIMMY  OF  AUTOMOBILE  WHEELS 

Ever  since  the  first  automobiles  started  to  appear  in  public 
certain  cars  revealed  the  existence  of  a  very  peculiar  oscillatory 
phenomenon  of  the  front  wheels  consisting  of  a  periodic  lateral 
wavering.   This  phenomenon  generally  occurred  when  the  car  .reached  or 
exceeded  certain  critical  speeds.   At  first,  it  was  called  waddling  or 
tacking.   In  the  beginning,  it  seemed  to  present  no  particularly  serious 
problem  and  car  designers  paid  little  attention  to  it.   Moreover,  at 
the  slow  speeds  of  cars  in  that  era,  it  was  successfully  overcome  by 
all  sorts  of  empirical  means,  among  which  might  be  mentioned  at  random: 

Toein  of  front  wheels 

Increasing  the  angle  of  camber 

Use  of  tight  shock  absorbers 

Stiff  front  springs 

Stiff  and  irreversible  steering 

Higher  tire  pressure,  etc. 
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All  these  attempts  were  no  more  than  palliatives,  the  simple  effect 
of  which  was  to  put  the  critical  speed  for  the  appearance  of  shimmy 
above  the  speed  of  utilization  of  the  car.   Shifting  the  critical  speed 
was  comparatively  easy,  except  for  racing  cars  where  more  energetic 
remedies  were  called  for.   The  only  truly  effective  remedy  for  this 
type  of  car  was  the  use  of  very  stiff  suspensions  and  hard  shock 
absorbers  which  meant  the  complete  loss  of  comfort. 

But  this  state  of  affairs  became  suddenly  worse  when  the  Americans 
tried  to  introduce  the  balloon  tire.   The  much  greater  elasticity  of 
this  tire,  so  pleasant  for  greater  comfort,  actually,  lowered  the 
critical  speed  of  oscillation  of  the  front  wheels  considerably  and. this 
was  aggravated  as  the  suspension  became  more  flexible  and  front  wheel 
brakes  came  into  wider  use.   On  some  cars,  these  phenomena  were  so 
violent  that  it  became  impossible  to  speed  up  even  moderately.   The 
oscillation  was  transmitted  integrally  to  the  axle  and,  because  of 
gyroscopic  interactions,  was  not  just  limited  to  a  simple  motion  in  a 
horizontal  plane  but  also  in  the  vertical  plane,  which  accentuated  and 
compounded  the  manifestations  to  an  unusual  degree.   In  some  cases  the 
axle  seemed  to  do  a  fervent  acrobatic  dance,  and  this  is  why  the 
Americans  got  to  designate  the  whole  of  these  phenomena  by  the  name 
of  an  exotic  dance  then  in  vogue  in  the  U.  S.  A.  and  the  name  shimmy 
has  remained  classical  ever  since. 

For  some  years  shimmy  was  regarded  as  a  particularly  grave  ailment 
of  balloon  tires  and  even  as  a  redhibitory  vice  to  their  commerical 
distribution.   Nevertheless,  the  greater  riding  comfort  on  balloon 
tires  tempted  the  engineers  of  every  country  to  all  sorts  of  ways  to 
combat  this  parasitic  and  annoying  phenomenon.   But  the  technicians 
themselves  did  not  agree  on  the  origin  and  nature  of  shimmy,  and  even 
less  on  a  rational  method  of  avoiding  its  effects.   The  question 
reached  such  importance  that  a  whole  session  of  the  U.  C.  Congress  of 
Automobile  Engineers  in  1925?  was  devoted  to  it.   Everybody  stuck 
loyally  to  the  results  of  his  own  observations  and  to  his  personal 
opinion.   But  the  diversity  of  opinion  of  assumptions  and  remedies 
suggested,  the  utter  lack  of  a  general  viewpoint  on  the  nature  of  the 
subject,  far  from  bringing  any  enlightenment,  actually  left  the  impres- 
sion of  profound  inconsistency.   It  is  only  necessary  to  consult  "La 
technique  automobile  et  aerienne,"  vol.  16,  no.  130,  3rd  quarter  1925) 
which  published  the  principal  reports  read  at  that  session.   It  is 
particularly  surprising  to  find  that  not  one  referred  to  the  true  origin 
of  shimmy.   The  underlying  cause  of  shimmy  lies  in  the  peculiar  mechani- 
cal properties  of  the  tire.   It  is  a  complex  shimmy,  the  sizeable 
elasticity  of  the  longitudinal  suspension  springs  plays  the  part  of  the 
lateral  elasticity  of  the  pivot. 

However,  the  complete  mathematical  theory  of  this  shimmy  of  auto- 
mobile wheels  is  more  complicated  than  that  of  complex  shimmy  described 
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in  the  foregoing  and  its  external  manifestations  themselves  are  a  little 
more  dissimilar  by  reason  of  the  following  two  factors: 

(1)  The  reactions  of  the  linkage  or  of  the  steering  gear 

(2)  The  gyroscopic  reactions 

The  first  idea  of  seeking  the  cause  of  shimmy  in  the  mechanical  pro- 
perties of  the  tire  is  much  earlier  than  the  recent  reports  of.  American 
experimenters  on  airplane  tires.   On  May  22,  1925  a  French  engineer, 
known  for  his  many  publications  on  the  dynamics  of  the  automobile  and 
for  his  frequently  original  views,  A.  Broulhiet,  gave  a  very  interesting 
speech  before  the  Societe  des  Ingenieus  Civile  de  France  (see  the 
report  of  this  society  (1925,  PP-  5^0-55^))  in  which  he  stated  that 
the  drift  of  the  tire  was  the  true  cause  of  shimmy.   Perhaps  Broulhiet 's 
idea  was  not  convincing  enough,  perhaps  it  was  not  accompanied  by  a 
sufficiently  correct  mathematical  theory  nor  developed  enough  to  have 
attention  called  to  it  as  the  only  correct  explanation  of  shimmy.   In 
any  case,  it  was  not  taken  up  by  other  theorists  and  did  not  receive 
the  publicity  it  should  have  had.   Incidentally,  the  development  of 
automobile  engineering  has  a  rather  unusual  history  and  which  proves 
that  even  a  wrong  theory  can  sometimes  be  very  practical. 

At  this  meeting  of  the  Society  of  American  Automobile  Engineers 
in  1925,  A.  Healy,  of  the  Dunlop  Tire  Company,  presented  a  report  in 
which  he  charged  the  gyroscopic  effect  of  the  wheels  as  being  the 
cause  of  shimmy.   The  idea  was  taken  up  in  France  by  an  engineer  of 
Russian  extraction,  Dimitri  Sensaud  De  Lavaud,  a  specialist  in  automobile 
dynamics,  who  undertook  to  publish  a  complete  mathematical  theory  of 
this  gyroscopic  effect  in  two  installments  (refs.  6  and  7). 

It  is  true  that  such  a  study  is  of  technical  interest,  because 
the  gyroscopic  reactions  have  an  undesirable  influence  on  the  shimmying 
motions  if  the  two  front  wheels  are  mounted  on  the  same  axle.   But 
they  can  in  no  case  be  the  true  cause  of  shimmy  because  they  are  con- 
servers  of  energy.   Consequently  the  gyroscopic  effect  absolutely  can- 
not explain  the  fundamental  characteristic  of  shimmy  which  is  the 
divergence  or  instability  of  oscillations  resulting  from  a  continuous 
transformation  of  kinetic  energy  of  translation  of  the  car  into  oscil- 
lating energy  of  the  wheel. 

This  matter  did  not  escape  De  Lavaud 's  attention,  who  explicitly 
remarked  on  it.   To  explain  the  divergence  of  the  oscillations  he 
attempted  to  bring  in  the  rocking  motion  which  naturally  accompanies 
shimmying.   To  be  sure,  there  is  a  real  possibility  of  divergence  but 
this  scarcely  perceptible  rocking  motion  is  so  minute  that  its  amount 
of  energy  per  cycle  is  very  small  and  surely  too  small  to. explain  the 
intensity  of  shimmy  especially  when  considering  the  many  sources  of 
damping  which  restrain  the  turning  motion  of  the  wheel. 
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In  De  Lavaud's  theory  this  imperceptible  rocking  motion  of  the 
car  is  the  sole  possible  cause  of  the  divergence  of  oscillation.   But 
it  is  plainly  inadequate.   The  author  was  therefore  compelled  to  regard 
shimmy  not  as  a  phenomenon  of  instability,  that  is  divergence  of  oscil- 
lations, but  as  a  simple  resonance  phenomenon  arising  from  defective 
balancing  or  centering  of  the  wheels.   Today  it  is  known  as  an  absolute 
fact  that  shimmy  is  not  a  resonance  phenomenon  but  due  solely  to  the 
divergence  of  oscillations,  which  is  entirely  different.   Even  if  it 
is  true  that  the  unbalance  of  the  front  wheels  may  facilitate  and,  in 
certain  cases,  actually  cause  artificial  shimmy,  experience  proves 
that  defective  balancing  of  wheels  or  tires,  whether  automobile  or 
airplane,  is  not  the  true  cause  of  shimmy.   No  resonance  phenomenon 
could  ever  explain  the  intensity  so  often  observed.   Furthermore, 
shimmy  has  been  observed  on  perfectly  balanced  wheels  rolling  on  level 
ground.   In  addition  the  phenomenon  designated  by  the  Americans  with 
the  name  of  kinematic  shimmy,  indicates  clearly,  that  it  can  have  an 
oscillatory  motion  without  the  least  effect  of  inertia,  which,  however, 
is  necessary  in  order  to  have  resonance. 

To  sum  up,  the  gyroscopic  actions  which  accompany  automobile  wheel 
shimmy  and  which  complicate  the  phenomenon  can  never  be  entirely  the 
cause.   This  is  readily  apparent  when  considering  airplane  wheel  shimmy 
studied  under  the  name  of  elementary  shimmy.   In  this  case  the  gyro- 
scopic actions  are  actually  without  any  effect  on  the  motion  of  the 
wheel,  since  their  moment  is  then  perpendicular  to  the  pivoting  axis 
which  itself  is  absolutely  rigid. 

In  short,  De  Lavaud  mathematically  defined  the  concept  of  the 
gyroscopic  coupling  of  two  oscillations,  very  important,  it  is  true, 
but  which  did  not  evolve  the  true  source  of  energy  brought  into  play 
by  shimmy.   As  a  logical  consequence  of  his  theory,  De  Lavaud  was 
forced  to  advocate  the  elimination  of  the  gyroscopic  effect  by  an 
appropriate  mode  of  linkage  of  the  front  wheels  to  the  chassis  per- 
mitting the  wheels  only  a  vertical  displacement  in  their  own  plane. 
Through  this  De  Lavaud  became  the  champion  and  forerunner  of  this  type 
of  suspension,  which  under  the  name  of  independent  front  wheel  suspen- 
sion achieved  a  certain  amount  of  success  for  two  reasons: 

First,  the  independent  front  wheel  eliminated,  as  foreseen,  the 
gyroscopic  reactions  and  their  vitiating  effects,  not  only  in  the  case 
of  shimmy  but  also  for  the  passage  of  the  wheel  over  any  obstacle,  and 
this  advantage  alone  was  enough  to  justify  its  use.   But,  and  this  is 
a  curious  historical  fact,  the  independent  front  wheels  also  completely 
eliminated  shimmy.   In  the  face  of  these  advantages  one  after  the  other 
automobile  builder  adopted  this  mode  of  suspension.   So,  with  the  pro- 
blem of  shimmy  practically  solved,  the  automobile  engineers  dropped  its 
study. 
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Without  any  other  justification  it  was  accepted  as  an  accomplished 
fact  that  the  gyroscopic  effects  were  the  true  cause  of  this  phenomena, 
and  up  to  this  day  the  theory  of  De  Lavaud  is  regarded  as  complete 
explanation.   But  now  it  has  been  proved  that  the  gyroscopic  actions 
are  simply  a  phenomenon  accompanying  but  not  a  primary  cause  of  shimmy. 
How  then  can  the  fact  that  the  independent  wheels  eliminate  shimmy  be 
explained? 

Actually,  all  independent  front  wheel  suspension  systems  have  in 
common  this  peculiarity  of  totally  suppressing  the  lateral  elasticity 
of  the  front  wheels,  a  lateral  elasticity  which  was  formerly  consider- 
able with  the  longitudinal  springs  generally  used  for  suspending  the 
chassis  from  the  axle  common  to  the  two  wheels.   Owing  to  this  fact, 
automobile  wheels  can  develop  only  elementary  shimmy  with  a  restoring 
torque  due  to  the  elasticity  of  the  steering  gear.   The  law  of  motion 
then  follows  equation  (33).   To  establish  the  stability  conditions, 
reference  is  made  to  the  discussion  following  this  equation.   In  auto- 
mobile wheels  with  a  given  caster  length 

aD  <  T 
It  follows  that  the  unique  condition  of  stability  is  v  <  u,  that  is 

o  ^  a      PT 

IR   (a  +  e)lR 

To  find  the  conditions  of  stability  at  any  speed  the  car  may  reach 
simply  consider  a  direction  in  which  the  steering  is  rigid  or  in  other 
words,  large  enough  p.   This  is  precisely  the  combination  realized 
on  all  modern  cars. 

VIII.  CONCLUSIONS 


The  present  report,  instead  of  starting  with  long  series  of  experi- 
ments and  involving  costly  materials,  begins  with  the  discussion  and 
exploitation  of  all  previous  experiments  in  enough  detail  to  be  prac- 
ticable.  Life  is  short  and  it  serves  no  useful  purpose  to  lose  time 
in  repeating  the  experiments  already  made  with  care  and  of  which  the 
results  have  been  published.   This  method  of  work  certainly  saves  time 
and  money,  because,  taking  into  consideration  that  the  comparison  of 
the  present  theory  with  the  many  results  obtained  earlier  by  other 
experimenters,  could  leave  no  doubt  about  the  physical  value  of  this 
theory  and  it  is  this  physical  value  which  is  the  principal  achievement 
of  this  work. 
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However,  some  experiments  were  carried  out  by  the  O.N.E.R.A.  at 
Toulouse  in  order  to  support  this  study  which  involved  an  experimental 
check  of  the  fundamental  stability  equation  (23a). 

(iRv2  -  a)(aD  -  T)  >  0 

This  formula  can  be  checked  by  ascertaining  that  the  conditions  of 
stability  and  those  of  shimmy  in  terms  of  caster  length  a  and  velo- 
city v  are  those  represented  in  figure  k-2.      The  experiments  themselves 
merely  involve  checking  the  pressure  or  absence  of  shimmy,  hence  com- 
paratively easy  and  inexpensive  measurements,  which  are,  in  any  case, 
less  difficult  and  shorter  than  Kantrowitz's  divergence  measurements. 

In  all  cases,  the  practical  use  of  this  theory  requires  the 
preliminary  and  systematic  measurement  of  the  characteristic  factors  T, 
D,  S,  R,  and  e  of  the  different  tires  used  in  aeronautics.   It  also  is 
desirable  to  study  the  variations  of  these  factors  in  terms  of  tire 
inflation  and  load.   Such  a  study  would  make  it  possible  to  decide  the 
form,  size,  and  inflation  pressure  required  to  obtain  the  minimum  for 
the  quotient  T/D,  so  as  to  be  able  to  combat  shimmy  without  having  to 
resort  to  long  caster  lengths. 

To  conclude,  it  is  a  theory  which  seems  to  have  already  been 
supported  by  many  experimental  proofs  and  which  is  herewith  presented 
to  tire  technicians  and  experimenters;  I  now  let  them  speak. 


Translated  by  J.  Vanier 
National  Advisory  Committee 
for  Aeronautics 
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Figure  25.-    Definition  of  forms.     Illustrated  directions  positive. 
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Wheel   central    plane 
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Figure  26.-    Geometrical  relations  involved  in  kinematic  shimmy,  illustrating 

that    dX  =  -eds. 
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Figure  40 
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